1. Multiple Integration 


1. 
. Double Integrals over Rectangular Regions 
. Double Integrals over General Regions 

. Double Integrals in Polar Coordinates 

. Triple Integrals 


CON MU BW N 


Introduction 


. Calculating Centers of Mass and Moments of Inertia 
. Change of Variables in Multiple Integrals 


Introduction 
class="introduction" 


The City of 
Arts and 
Sciences in 
Valencia, 
Spain, has a 
unique 
structure 
along an 
axis of just 
two 
kilometers 
that was 
formerly the 
bed of the 
River Turia. 
The 
1’ Hemisfeéric 
has an 
IMAX 
cinema with 
three 
systems of 
modem 
digital 
projections 
onto a 
concave 
screen of 
900 square 
meters. An 
oval roof 
over 100 
meters long 
has been 


made to 
look like a 
huge human 
eye that 
comes alive 
and opens 
up to the 
world as the 
“Eye of 
Wisdom.” 
(credit: 
modificatio 
n of work 
by Javier 
Yaya Tur, 
Wikimedia 
Commons) 


In this chapter we extend the concept of a definite integral of a single 
variable to double and triple integrals of functions of two and three 
variables, respectively. We examine applications involving integration to 
compute volumes, masses, and centroids of more general regions. We will 
also see how the use of other coordinate systems (such as polar, cylindrical, 
and spherical coordinates) makes it simpler to compute multiple integrals 
over some types of regions and functions. As an example, we will use polar 
coordinates to find the volume of structures such as |’ Hemisféric. (See 
[link].) 


In the preceding chapter, we discussed differential calculus with multiple 
independent variables. Now we examine integral calculus in multiple 
dimensions. Just as a partial derivative allows us to differentiate a function 
with respect to one variable while holding the other variables constant, we 
will see that an iterated integral allows us to integrate a function with 
respect to one variable while holding the other variables constant. 


Double Integrals over Rectangular Regions 


e Recognize when a function of two variables is integrable over a rectangular region. 

e Recognize and use some of the properties of double integrals. 

e Evaluate a double integral over a rectangular region by writing it as an iterated integral. 

e Use a double integral to calculate the area of a region, volume under a surface, or average value of a function 
over a plane region. 


In this section we investigate double integrals and show how we can use them to find the volume of a solid over a 
rectangular region in the zy-plane. Many of the properties of double integrals are similar to those we have already 
discussed for single integrals. 


Volumes and Double Integrals 


We begin by considering the space above a rectangular region R. Consider a continuous function f (x, y) > 0 of 
two variables defined on the closed rectangle R: 
Equation: 


R=[a,0] x [c,d] = {(2,y) € R’la<a<bc<y<d} 


Here [a,b] x [c, d] denotes the Cartesian product of the two closed intervals [a, b] and |c, d]. It consists of 
rectangular pairs (x, y) such thata < « < bandc < y < d. The graph of f represents a surface above the ry- 
plane with equation z = f (x, y) where z is the height of the surface at the point (, y). Let S' be the solid that lies 
above R and under the graph of f ((link]). The base of the solid is the rectangle R in the xy-plane. We want to find 
the volume V of the solid S. 


Zh 


The graph of f (x, y) over the rectangle R in the ry- 
plane is a curved surface. 


We divide the region R into small rectangles R,;, each with area AA and with sides Az and Ay ([link]). We do 
this by dividing the interval [a, b] into m subintervals and dividing the interval [c, d] into n subintervals. Hence 


= ArAy. 


and AA 


n? 


d—c 


Rectangle R is divided into small rectangles R,;, each with area AA. 
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The volume of a thin rectangular box above R 
each R,; as shown in the following figure. 


Zh 


A thin rectangular box above R;; with height 


x x 


f (x55, Yi): 


Using the same idea for all the subrectangles, we obtain an approximate volume of the solid S as 


mm n * * 
Ve S S f (ij. yi;) 4A. This sum is known as a double Riemann sum and can be used to approximate the 
i=1 j=l 
value of the volume of the solid. Here the double sum means that for each subrectangle we evaluate the function at 
the chosen point, multiply by the area of each rectangle, and then add all the results. 


As we have seen in the single-variable case, we obtain a better approximation to the actual volume if m andn 
become larger. 


Equation: 
mn 5. " min 3 ‘ 
= ii Vig AA =i ye )AA. 
Y= ites 2y dy Mew vgAAorV= lim | > > fey dy) 


Note that the sum approaches a limit in either case and the limit is the volume of the solid with the base R. Now we 
are ready to define the double integral. 


Note: 

Definition 

The double integral of the function f (x, y) over the rectangular region R in the xy-plane is defined as 
Equation: 


| / f(2,y)dA = lim Pty yj) AA. 


. i=1 j=l 


If f (x, y) = 0, then the volume V of the solid S, which lies above R in the xy-plane and under the graph of f, is 
the double integral of the function f (a, y) over the rectangle R. If the function is ever negative, then the double 
integral can be considered a “signed” volume in a manner similar to the way we defined net signed area in The 
Definite Integral. 


Example: 
Exercise: 


Problem: 
Setting up a Double Integral and Approximating It by Double Sums 


Consider the function z = f(x, y) = 3x? — y over the rectangular region R = (0, 2] x (0, 2] ({link]). 


a. Set up a double integral for finding the value of the signed volume of the solid S that lies above R and 
“under” the graph of f. 

b. Divide R into four squares with m = n = 2, and choose the sample point as the upper right corner 
point of each square (1,1), (2, 1), (1, 2), and (2, 2) ([link]) to approximate the signed volume of the 
solid S that lies above R and “under” the graph of f. 

c. Divide R into four squares with m = n = 2, and choose the sample point as the midpoint of each 
square: (1/2,1/2), (3/2, 1/2), (1/2, 3/2), and (3/2, 3/2) to approximate the signed volume. 

z=f(x,y) = 3 -y 


The function z = f(z, y) graphed over the 
rectangular region R = [0,2] x [0,2]. 


Solution: 


a. As we can see, the function z = f(z, y) = 3x? — y is above the plane. To find the signed volume of S, 
we need to divide the region R into small rectangles R,;, each with area A.A and with sides Ax and 
Ay, and choose (2;;,4;;) as sample points in each R;;. Hence, a double integral is set up as 
Equation: 


V= // (30? —y)dA= lim y s [3(«:,) = u AA. 
R 


i=1 j=l 


b. Approximating the signed volume using a Riemann sum with m = n = 2 we have 
AA = AvAy =1 x 1 = 1. Also, the sample points are (1, 1), (2, 1), (1, 2), and (2, 2) as shown in the 
following figure. 


Choose upper right 
corners of squares 


y ——n 


(1, 2) (2, 2) 


(2, 1) 


Each square has area AA = 1 


Subrectangles for the rectangular region 
A — (0,2) 0) 2 1 


Hence, 
Equation: 


2 * 
PPC 2i7¥,,)AA 


j=l 


IM» iM 


(F(ai1, Yai) + F(a, Yin) AA 


s. 
ll 
i 


= Fle, Yi JAA + f(91, Yn JAA + f(212, Yi2JAA + F(222, 422) AA 
= f(1,1)(1) + F(2,1)) + F, 2)(1) + F(2,2)(1) 

= (3 — 1)(1) + (12 — 1)(1) + (8 — 2)(1) + (12 — 2)(1) 
=24+114+1+410= 24. 


c. Approximating the signed volume using a Riemann sum with m = n = 2, we have 
AA = AzvAy=1 x 1 = 1. In this case the sample points are (1/2, 1/2), (3/2, 1/2), (1/2, 3/2), 
and (3/2, 3/2). 
Hence 
Equation: 


= f(@1,y,)AA + f(z a ¥aAA + f(ai2, ¥i2)AA + f(@29, Yor) AA 
= f(1/2,1/2)(1) + f(3/2,1/2)(1) + f(1/2, 3/2)(1) + (3/2, 3/2)(1) 
(=) a (4 Ne ee 2 (eae 2) (0) 


Analysis 


Notice that the approximate answers differ due to the choices of the sample points. In either case, we are 
introducing some error because we are using only a few sample points. Thus, we need to investigate how we can 
achieve an accurate answer. 


Note: 
Exercise: 


Problem: Use the same function z = f(z, y) = 3x? — y over the rectangular region R = [0,2] x [0,2]. 


Divide R into the same four squares with m = n = 2, and choose the sample points as the upper left corner 
point of each square (0,1), (1,1), (0, 2), and (1, 2) ([link]) to approximate the signed volume of the solid S 
that lies above # and “under” the graph of f. 


Solution: 
2 ®@ ae 
V= Se SS F(2ij,Yij)AA = 0 
i=l j=l 
Hint 


Follow the steps of the previous example. 


Note that we developed the concept of double integral using a rectangular region R. This concept can be extended 
to any general region. However, when a region is not rectangular, the subrectangles may not all fit perfectly into R, 
particularly if the base area is curved. We examine this situation in more detail in the next section, where we study 
regions that are not always rectangular and subrectangles may not fit perfectly in the region R. Also, the heights 
may not be exact if the surface z = f («, y) is curved. However, the errors on the sides and the height where the 
pieces may not fit perfectly within the solid S approach 0 as m and n approach infinity. Also, the double integral of 
the function z = f (2, y) exists provided that the function f is not too discontinuous. If the function is bounded 
and continuous over R except on a finite number of smooth curves, then the double integral exists and we say that 
f is integrable over R. 


Since AA = AzAy = AyAz, we can express dA as dx dy or dy dz. This means that, when we are using 


rectangular coordinates, the double integral over a region R denoted by | / f(x, y)dA can be written as 


|] tear ay or f f(x, y)dy dx. 
R R 


R 


Now let’s list some of the properties that can be helpful to compute double integrals. 


Properties of Double Integrals 


The properties of double integrals are very helpful when computing them or otherwise working with them. We list 
here six properties of double integrals. Properties 1 and 2 are referred to as the linearity of the integral, property 3 
is the additivity of the integral, property 4 is the monotonicity of the integral, and property 5 is used to find the 
bounds of the integral. Property 6 is used if f(z, y) is a product of two functions g(a) and h(y). 


Note: 

Properties of Double Integrals 

Assume that the functions f(a, y) and g(«, y) are integrable over the rectangular region R; S and T are subregions 
of R; and assume that m and M are real numbers. 


i. The sum f(z, y) + g(2, y) is integrable and 


Equation: 
k If Fe. 9) + ole, ania = h If sean | ff eterna. 


ii. If c is a constant, then cf (a, y) is integrable and 


Equation: 
|[etewaa = c ff #(e,yaa. 
R R 


iii. If R = SUT and ST = 0 except an overlap on the boundaries, then 


Equation: 
|] tenaa= ff tena ff teyaa. 
R S ae 


iv. If f(z, y) > g(a, y) for (x, y) in R, then 


Equation: 
|] tena 2 |[sevaa. 
R R 


v. lfm < f(x,y) < M, then 
Equation: 


m x A(R) < |[ #enaas M x A(R). 
R 


E. 


vi. In the case where f(z, y) can be factored as a product of a function g(a) of x only and a function h(y) of y 


only, then over the region R = {(x, y)|a <  < b,c < y < d}, the double integral can be written as 


Equation: 
b d 
i s(e,s)aa= (f° aayae) ( [ ne)a). 


These properties are used in the evaluation of double integrals, as we will see later. We will become skilled in 
using these properties once we become familiar with the computational tools of double integrals. So let’s get to 
that now. 


Iterated Integrals 


So far, we have seen how to set up a double integral and how to obtain an approximate value for it. We can also 
imagine that evaluating double integrals by using the definition can be a very lengthy process if we choose larger 
values for m and n. Therefore, we need a practical and convenient technique for computing double integrals. In 
other words, we need to learn how to compute double integrals without employing the definition that uses limits 
and double sums. 


The basic idea is that the evaluation becomes easier if we can break a double integral into single integrals by 
integrating first with respect to one variable and then with respect to the other. The key tool we need is called an 
iterated integral. 


Note: 

Definition 

Assume a, b, c, and d are real numbers. We define an iterated integral for a function f(z, y) over the rectangular 
region R = [a,b] x [c,d] as 


a. 
Equation: 
b d bf d 
| [ teadyac= [ [teva dx 
b. 
Equation: 
d b d[— b 
[ [ tenacay= [| [ r(0,s)de} ay 


b d 


The notation / / f(x, y)dy | dx means that we integrate f(x, y) with respect to y while holding x constant. 
a c 
df b 
Similarly, the notation di / f(x, y)dx } dy means that we integrate f(x, y) with respect to x while holding y 
c a 


constant. The fact that double integrals can be split into iterated integrals is expressed in Fubini’s theorem. Think 
of this theorem as an essential tool for evaluating double integrals. 


Note: 

Fubini’s Theorem 

Suppose that f(z, y) is a function of two variables that is continuous over a rectangular region 

ie = AG y) €Rla<2<bc<y< Ue Then we see from [link] that the double integral of f over the region 
equals an iterated integral, 


Equation: 
[foeoma= [f semtear= J [romavae= ff eae 


More generally, Fubini’s theorem is true if f is bounded on R and f is discontinuous only on a finite number of 
continuous curves. In other words, f has to be integrable over R. 


~ 


wee we wm em em em em we ew ee ee ee ee eS 
~ 


seem ew eee eee eee ee ee ee = S 


Height 
f(x, y) 


Height 
f(x, y) 


Area A(x) 


(a) (b) 


(a) Integrating first with respect to y and then with respect to z to find the area A(x) and then the volume V; 
(b) integrating first with respect to x and then with respect to y to find the area A(y) and then the volume V. 


Example: 
Exercise: 


Problem: 
Using Fubini’s Theorem 


Use Fubini’s theorem to compute the double integral | / f(x, y)dA where f(x, y) = x and 
R 


#—|0,2| x (0,1): 


Solution: 


Fubini’s theorem offers an easier way to evaluate the double integral by the use of an iterated integral. Note 
how the boundary values of the region R become the upper and lower limits of integration. 
Equation: 


y=1 me 
= / 2dy = 2y|)—0 = 2, 
y 


The double integration in this example is simple enough to use Fubini’s theorem directly, allowing us to convert a 
double integral into an iterated integral. Consequently, we are now ready to convert all double integrals to iterated 
integrals and demonstrate how the properties listed earlier can help us evaluate double integrals when the function 
f(x,y) is more complex. Note that the order of integration can be changed (see [link]). 


Example: 
Exercise: 


Problem: 
Illustrating Properties i and ii 


Evaluate the double integral // (cy — 3ay”)dA where R = {(2,y)|0< a < 2,1 <y <2}. 
R 


Solution: 


This function has two pieces: one piece is zy and the other is 3ay”. Also, the second piece has a constant 3. 
Notice how we use properties i and ii to help evaluate the double integral. 
Equation: 


il i (Bea ee) 


R 


= // cydA + i (—3ay”)dA Property i: Integral of a sum is the sum of the int 
R R 


y=2 a=2 y=2 a2=2 
= i / cy dz dy — / 3xy"dax dy Convert double integrals to iterated integrals. 
y=1 z=0 y=1 z=0 


y=2 y=2 = 
= i (2y)|"=ody — 3 | (27) oy Integrate with respect to x, holding y constant. 
y=1 y=1 
y=2 y=2 
= i 2y dy — / 6y*dy Property ii: Placing the constant before the inte 
y=1 y=1 
2 2 
= y dy — 6 | y’dy Integrate with respect to y. 
1 1 
2 2 
= 25 |, om 64 1 
a2) 2 3 |2 
=y"|; — 29°]; 


= (d=1) = 2(8— 1) 
=3— 37) = 3-14 = 11. 


Example: 
Exercise: 


Problem: 
Illustrating Property v. 


Over the region R = {(z, y)|1 <a < 3,1 <y< 2}, we have 2 < x? + y? < 13. Finda lower and an 


upper bound for the integral // (ae 4p y’) dA. 
R 


Solution: 


For a lower bound, integrate the constant function 2 over the region R. For an upper bound, integrate the 
constant function 13 over the region R. 


Equation: 
2 3 2 2 
i i Qde dy = / [zelf]ay =f 2(2)dy = 4y|? = 4(2-1) =4 
1 1 1 1 


2 3 2 2 
/ / 13dz dy = | [1320/3] dy - i: 13(2)dy = 26y|? = 26(2 — 1) = 26. 
1 1 1 1 


Hence, we obtain 4 < // (ae -- y)dA SHG. 
R 


Example: 
Exercise: 


Problem: 
Illustrating Property vi 


Evaluate the integral // e’cos x dA over the region R = { (x,y) |0 SS ES ine 
R 


Solution: 


This is a great example for property vi because the function f(z, y) is clearly the product of two single- 
variable functions e¥ and cos x. Thus we can split the integral into two parts and then integrate each one as a 
single-variable integration problem. 


Equation: 
1 1/2 
// e’cosrdA = i i e’cos x dx dy 
- 0 JO 
1 1/2 
= (/ a) i cos x dx 
0 0 
= (cla) (sin 2|9””) 
=e—I. 
Note: 
Exercise: 
Problem: 


a. Use the properties of the double integral and Fubini’s theorem to evaluate the integral 


Equation: 
1 3 
a i (3 — «+ 4y)dy da. 
0 J-1 
1 


b. Show that 0 < [fs mz cosTy dA < ag where R = (0, =) = =). 
R 


Solution: 


a. 26 b. Answers may vary. 


Hint 


Use properties i. and ii. and evaluate the iterated integral, and then use property v. 


As we mentioned before, when we are using rectangular coordinates, the double integral over a region R denoted 


by // f(x, y)dA can be written as // f(x, y)dx dy or // f(x, y)dy dx. The next example shows that the 
R R R 


results are the same regardless of which order of integration we choose. 


Example: 
Exercise: 


Problem: 
Evaluating an Iterated Integral in Two Ways 


Let’s return to the function f(z, y) = 3x? — y from [link], this time over the rectangular region 


R= [0,2] x [0,3]. Use Fubini’s theorem to evaluate i / f(x, y)dA in two different ways: 
R 


a. First integrate with respect to y and then with respect to x; 
b. First integrate with respect to x and then with respect to y. 
Solution: 
[link] shows how the calculation works in two different ways. 


a. First integrate with respect to y and then integrate with respect to x: 
Equation: 


|[ tenes =f . / "(be" — y)dy dx 
R 


y = 
2=2 = 
— / (9a? 2\dx = oy —t2\" 5 = 15 
z=0 
b. First integrate with respect to x and then integrate with respect to y: 
Equation: 
y=3 2=2 
ff senaa =f foe? pave ay 
y=0 a=0 
R 
y=3 2=2 y=3 a 
= / (/ ie = ide dy = / Es a yl | dy 
y=0 «=0 y=0 7 
y=3 =, 
=| (8 — 2y)dy = 8y — |" = 15. 
y=0 : 
Analysis 


With either order of integration, the double integral gives us an answer of 15. We might wish to interpret this 
answer as a volume in cubic units of the solid S below the function f(z, y) = 3x7 — y over the region 


R= [0,2] x [0,3]. However, remember that the interpretation of a double integral as a (non-signed) volume 
works only when the integrand f is a nonnegative function over the base region R. 


Note: 
Exercise: 


iy) a=5 

Problem: Evaluate / / (2 Poa y”) dex dy. 
a=—s) des 

Solution: 


__ 1340 
3 


Hint 


Use Fubini’s theorem. 


In the next example we see that it can actually be beneficial to switch the order of integration to make the 
computation easier. We will come back to this idea several times in this chapter. 


Example: 
Exercise: 


Problem: 
Switching the Order of Integration 


Consider the double integral le sin(xy)dA over the region R = {(x, y)|0 < x < 3,0 < y < 2} ([link)). 


a. Express the double integral in two different ways. 
b. Analyze whether evaluating the double integral in one way is easier than the other and why. 
c. Evaluate the integral. 


The function z = f(x, y) = x sin(ay) over the rectangular region 
R—|0;x| x |l,2h 


Solution: 
a. We can express | / zx sin(zy)dA in the following two ways: first by integrating with respect to y and 


R 
then with respect to x; second by integrating with respect to x and then with respect to y. 
Equation: 


| i ren(aydA 


x sin(xy)dy dx Integrate first with respect to y. 


x sin(xy)dx dy Integrate first with respect to z. 


b. If we want to integrate with respect to y first and then integrate with respect to x, we see that we can use 
the substitution u = xy, which gives du = x dy. Hence the inner integral is simply / sin u du and we 


can change the limits to be functions of x, 
Equation: 


[fesiacena a- ff 


R z 


il x sin(xy)dy dz = i fi sin(u)du } dx. 
=1 


y z=0 U=z 


However, integrating with respect to x first and then integrating with respect to y requires integration by 
parts for the inner integral, with uw = « and du = sin(ry)dz. 


Then du = dz andv = — teu) sO 
Equation: 
te l—1 coeds 
x cos(xy) 1 
x sin(xy)d xz sin(xy)dz dy = cos(xy)dz } dy. 
R 1 y 2=0 y ; 
We iy 


Since the evaluation is getting complicated, we will only do the computation that is easier to do, which 
is clearly the first method. 
c. Evaluate the double integral using the easier way. 


Equation: 
aan y=2 
[[esincenaa = i, i} x sin(xy)dy dx 
R z=0 y=l 
van [ u=2e von can 
= / / sin(u)du| dx = i |—cos u aa ike = i (—cos 2x + cos x)dx 
z=0 |u=« 2=0 z=0 
= —1sin 22 + sin z|F_5 = 0. 
Note: 
Exercise: 


Problem: Evaluate the integral H l ze“4dA where R = [0,1] x [0,1In 5). 


Solution: 


4—In5 
Ind 


Hint 


Integrate with respect to y first. 


Applications of Double Integrals 


Double integrals are very useful for finding the area of a region bounded by curves of functions. We describe this 
situation in more detail in the next section. However, if the region is a rectangular shape, we can find its area by 


integrating the constant function f(x, y) = 1 over the region R. 


Note: 

Definition 

The area of the region R is given by A(R) = | i 1dA. 
R 


This definition makes sense because using f(x, y) = 1 and evaluating the integral make it a product of length and 
width. Let’s check this formula with an example and see how this works. 


Example: 
Exercise: 


Problem: 
Finding Area Using a Double Integral 


Find the area of the region R = {(x,y)|0 < « < 3,0 < y < 2} by using a double integral, that is, by 
integrating 1 over the region R. 


Solution: 


The region is rectangular with length 3 and width 2, so we know that the area is 6. We get the same answer 
when we use a double integral: 


Equation: 
H 3 2 2 
A(R) = ff idedy= f {eli)av= [ say=3 
0 


0 0 0 


dy = 3y|3 = 3(2) =6. 


— 


We have already seen how double integrals can be used to find the volume of a solid bounded above by a function 
f(x, y) over a region R provided f(x, y) > 0 for all (x, y) in R. Here is another example to illustrate this concept. 


Example: 
Exercise: 


Problem: 
Volume of an Elliptic Paraboloid 


Find the volume V of the solid S that is bounded by the elliptic paraboloid 2x? + y? + z = 27, the planes 
x = 3 and y = 3, and the three coordinate planes. 


Solution: 
First notice the graph of the surface z = 27 — 2x” — y? in [link](a) and above the square region 


R, = [-3,3] x [—3,3]. However, we need the volume of the solid bounded by the elliptic paraboloid 
2a? + y? + z = 27, the planes x = 3 and y = 3, and the three coordinate planes. 


(a) The surface z = 27 — 2x? — y* above the square region R; = [—3,3] x [—3,3]. (b) The solid S 
lies under the surface z = 27 — 2x” — y* above the square region R2 = [0,3] x [0,3]. 


Now let’s look at the graph of the surface in [link](b). We determine the volume V by evaluating the double 
integral over Ro: 


Equation: 
Vv = |[-4a- ff er-2 —y’)dA 
R R 
y=3 2=3 
= / i (27 aU = y’) dx dy Convert to iterated integral. 
g=0 2=0 
v=3 
3 27 |2=3 : 
= [27x =a” =D) z| ae Integrate with respect to x. 
y=0 
y=3 
= / (64 — 3y”)dy = 63y — ae = 162. 
y=0 
Note: 
Exercise: 
Problem: 


Find the volume of the solid bounded above by the graph of f(x, y) = xy sin(x?y) and below by the zy- 
plane on the rectangular region R = [0,1] x [0,z]. 


Solution: 


[3 


Hint 


Graph the function, set up the integral, and use an iterated integral. 


Recall that we defined the average value of a function of one variable on an interval [a, b] as 
Equation: 


b 
1 
fave => i [ teae. 


Similarly, we can define the average value of a function of two variables over a region R. The main difference is 
that we divide by an area instead of the width of an interval. 


Note: 

Definition 

The average value of a function of two variables over a region R is 
Equation: 


fave = sag If flewaa. 
R 


In the next example we find the average value of a function over a rectangular region. This is a good example of 
obtaining useful information for an integration by making individual measurements over a grid, instead of trying to 
find an algebraic expression for a function. 


Example: 
Exercise: 


Problem: 
Calculating Average Storm Rainfall 


The weather map in [link] shows an unusually moist storm system associated with the remnants of Hurricane 
Karl, which dumped 4-8 inches (100-200 mm) of rain in some parts of the Midwest on September 22-23, 
2010. The area of rainfall measured 300 miles east to west and 250 miles north to south. Estimate the 
average rainfall over the entire area in those two days. 


Effects of Hurricane Karl, which dumped 4-8 inches (100-200 mm) of rain in some parts of southwest 
Wisconsin, southern Minnesota, and southeast South Dakota over a span of 300 miles east to west and 
250 miles north to south. 


Solution: 

Place the origin at the southwest corner of the map so that all the values can be considered as being in the 
first quadrant and hence all are positive. Now divide the entire map into six rectangles (m = 2 and n = 3), 
as shown in [link]. Assume f(z, y) denotes the storm rainfall in inches at a point approximately x miles to 
the east of the origin and y miles to the north of the origin. Let R represent the entire area of 


250 x 300 = 75000 square miles. Then the area of each subrectangle is 
Equation: 


1 
AA = = (75000) = 12500. 


* * 
Assume (x55, Yi;) are approximately the midpoints of each subrectangle R;;. Note the color-coded region at 
each of these points, and estimate the rainfall. The rainfall at each of these points can be estimated as: 


At (211, yi1) the rainfall is 0.08. 
At (212, yiz) the rainfall is 0.08. 
At (213, yi3) the rainfall is 0.01. 
At (£21, ya1) the rainfall is 1.70. 
At (292, Yo2) the rainfall is 1.74. 


At (#23, y23) the rainfall is 3.00. 


os aye 


ars (os 9 A RE ee 


Storm rainfall with rectangular axes and showing the midpoints of each subrectangle. 


According to our definition, the average storm rainfall in the entire area during those two days was 
Equation: 


foe = sake ff fe, wde dy = asian ff fe, v)de dy 
R R 


* * 


f (255; y;)AA 


IR 
a 
iM 
Me 


IIe 


75000 ies Yi JAA +t f(212, 412) AA 


+f(213) ¥13)AA + f(01, Ya JAA + F (2295 Yo2) AA + f (293, Yo3) NA] 
T5000 [0.08 + 0.08 + 0.01 + 1.70 + 1.74 + 3.00)AA 
mn [0.08 + 0.08 + 0.01 + 1.70 + 1.74 + 3.00]12500 


[0.08 + 0.08 + 0.01 + 1.70 + 1.74 + 3.00] 
10. 


~| 


me ge 


During September 22—23, 2010 this area had an average storm rainfall of approximately 1.10 inches. 


Note: 
Exercise: 


Problem: A contour map is shown for a function f(«, y) on the rectangle R = [—3,6] x [—1,4]. 


a. Use the midpoint rule with m = 3 and n = 2 to estimate the value of | it f(x, y)dA. 
R 


b. Estimate the average value of the function f(z, y). 


Solution: 


Answers to both parts a. and b. may vary. 
Hint 


Divide the region into six rectangles, and use the contour lines to estimate the values for f(z, y). 


Key Concepts 


e We can use a double Riemann sum to approximate the volume of a solid bounded above by a function of two 
variables over a rectangular region. By taking the limit, this becomes a double integral representing the 
volume of the solid. 

e Properties of double integral are useful to simplify computation and find bounds on their values. 

We can use Fubini’s theorem to write and evaluate a double integral as an iterated integral. 

Double integrals are used to calculate the area of a region, the volume under a surface, and the average value 

of a function of two variables over a rectangular region. 


Key Equations 


¢ Double integral 


n 


[f sea = tim SOY Heinui)aA 
R 2=1 


j=l 
¢ Iterated integral 


[ [ tenaeay= [ [[ rewau a 
hh ie f(z, y)dx dy = [ [fs i)ae| dy 


e Average value of a function of two variables 


foe = scr ff f2v)de dy 
R 


In the following exercises, use the midpoint rule with m = 4 and n = 2 to estimate the volume of the solid 
bounded by the surface z = f(z, y), the vertical planes ¢ = 1, x = 2, y = 1, and y = 2, and the horizontal plane 


2= 0; 
Exercise: 


Problem: f(z, y) = 4x + 2y + 8xy 
Solution: 
27. 

Exercise: 


Problem: f(x,y) = 16x” + 4 


In the following exercises, estimate the volume of the solid under the surface z = f(a, y) and above the 
rectangular region R by using a Riemann sum with m = n = 2 and the sample points to be the lower left corners 
of the subrectangles of the partition. 

Exercise: 


Problem: f(x,y) = sinz — cosy, R= [0,7] x [0,7] 


Solution: 


0. 


Exercise: 


Problem: f(x,y) = cos z+ cosy, R = [0,7] x [0, =| 
Exercise: 


Problem: 
Use the midpoint rule with m = n = 2 to estimate | / f(a, y)dA, where the values of the function f on 
R 


R= [8,10] x [9,11] are given in the following table. 


y 
X 9 9.5 10 10.5 11 
8 9.8 5 6.7 5 5.6 
8.5 9.4 4.5 8 5.4 3.4 
9 8.7 4.6 6 5.5 3.4 
9.5 6.7 6 4.5 5.4 6.7 


10 6.8 6.4 5.5 5.7 6.8 


Solution: 


2153; 
Exercise: 


Problem: 


The values of the function f on the rectangle R = [0,2] x [7,9] are given in the following table. Estimate 


the double integral | / f(x, y)dA by using a Riemann sum with m = n = 2. Select the sample points to be 


R 
the upper right corners of the subsquares of R. 


yo =7 yi =8 y2=9 
rp = 0 10.22 10.21 9.85 
rj=1 6.73 9.75 9.63 
L2=2 5.62 7.83 8.21 
Exercise: 
Problem: 


The depth of a children’s 4-ft by 4-ft swimming pool, measured at 1-ft intervals, is given in the following 
table. 


a. Estimate the volume of water in the swimming pool by using a Riemann sum with m = n = 2. Select 
the sample points using the midpoint rule on R = [0,4] x [0, 4]. 
b. Find the average depth of the swimming pool. 


DA 
x 0 1 2 3 + 
0 1 185 2 2.5 3 
1 1 1.5 2 2.5 3 
2 1 1.5 1.5 2.5 3 
3 1 1 1.5 2 2.5 


Solution: 


a. 28 ft® b. 1.75 ft. 
Exercise: 
Problem: 
The depth of a 3-ft by 3-ft hole in the ground, measured at 1-ft intervals, is given in the following table. 
a. Estimate the volume of the hole by using a Riemann sum with m = n = 3 and the sample points to be 


the upper left corners of the subsquares of R. 
b. Find the average depth of the hole. 


y 
x 0 1 2 3 
0 6 6.5 6.4 6 
1 6.5 7 vies) 6.5 
2 6.5 6.7 6.5 6 
3 6 6.5 5 5.6 
Exercise: 
Problem: 


The level curves f(z, y) = & of the function f are given in the following graph, where k is a constant. 


a. Apply the midpoint rule with m = n = 2 to estimate the double integral | / f(x, y)dA, where 
R 


R= (0.2,1] x [0,0.8]. 
b. Estimate the average value of the function f on R. 


1 1 

on: de 

== k=0 k=; 
1 1 

—k=> —k=5 —k=1 


Solution: 


a.0.112b. fave © 0.175; here f(0.4, 0.2) ~ 0.1, f(0.2,0.6) ~ —0.2, (0.8, 0.2) ~ 0.6, and 
(0.8, 0.6) ~ 0.2. 


Exercise: 


Problem: 


The level curves f(z, y) = & of the function f are given in the following graph, where k is a constant. 


a. Apply the midpoint rule with m = n = 2 to estimate the double integral | i f(a, y)dA, where 
R 


R= (0.1,0.5] x [0.1, 0.5]. 
b. Estimate the average value of the function f on R. 


0 01 02 03 04 05 06 0.708 09 1% 


1 a: 1 1 
~k=y Tks k=3 k=% 
1 _ 3 pe 
=e k=_ k=1 
Exercise: 
Problem: 


The solid lying under the surface z = \/4 — y? and above the rectangular region R = [0,2] x [0,2] is 


illustrated in the following graph. Evaluate the double integral | fl f(x, y)dA, where f(z, y) = /4 — y?, by 


R 
finding the volume of the corresponding solid. 


Solution: 


Qn. 


Exercise: 


Problem: 


The solid lying under the plane z = y + 4 and above the rectangular region R = [0,2] x [0,4] is illustrated 
in the following graph. Evaluate the double integral | f(x, y)dA, where f(x,y) = y+ 4, by finding the 


volume of the corresponding solid. 


In the following exercises, calculate the integrals by interchanging the order of integration. 
Exercise: 


1 2 
Problem i / (2a + 3y+ 5)dx | dy 
=) \=2 
Solution: 
40. 
Exercise: 


1 
Problem: / fe + 2e¥ — 3)dx | dy 
0 


Exercise: 


2 
Problem: i J (Wa + 3/y)dy | dx 
i 


Solution: 


+ 39V2. 


Exercise: 


16 


8 
Problem: / [ (e+ 2vaey dz 
1 


1 


Exercise: 
In3 1 
Problem: / / e* Ydy | dx 
In2 0 
Solution: 
e—l. 
Exercise: 
2 1 
Problem: / / 3° dy | dx 
0 0 
Exercise: 
6 95 
Problem: / / V4 dy dx 
x 
1 2 
Solution: 
10V2 
15 — ar aoe 
Exercise: 


ey a a 
Problem: / i Mo dz 
y 
1 \4 


In the following exercises, evaluate the iterated integrals by choosing the order of integration. 
Exercise: 
a 1/2 
Problem: / / sin(2x)cos(3y)dx dy 
0 0 


Solution: 


0. 


Exercise: 


1/8 1/3 
Problem: / / [cot « + tan(2y)|dzx dy 
m/12 4/4 


Exercise: 


e e 1 1 
problem: | f/| © sin(in x) + —cos(In y) | dx dy 
x y 
id 


Solution: 


(e —1)(1 + sin1 — cos 1). 


Exercise: 


e€ e€ 7 1 1 
Problem: / / ene eeny) dx dy 
ry 
1 1 


Exercise: 
2 2 
probiem: | [ (2! + scr )avee 
Tig 2 
Solution: 


31n (3) +21n?2 —In2. 
Exercise: 
e€ 


2 
Problem: [ f 2°In(x)dy dz 
1 1 


Exercise: 


V3 


2 
1 
Problem: / 7 y arctan (=) dy dz 
x 
1 


1 


Solution: 


1|(2v3-3)n+6In2]. 


Exercise: 


1 1/2 
Problem: / 7 (arcsin x + arcsin y)dy dx 
0 0 


Exercise: 


1 2 
Problem: / / xe? Ydy dx 
i 


0 


Solution: 
+e (e4 — ay 
Exercise: 
2 1 
Problem: / / ze” Ydydzx 
1 0 
Exercise: 
In tn Ing 
Problem: dy dz 
va Ve 
Solution: 
A(e — 1) (2- ve). 
Exercise: 
xin Yo 7 Ing 
Problem: dy dz 
JY 
Exercise: 
1 2 
Problem: / / (ata a) ty dx 
0 1 
Solution: 


= + In (4) sin 2 + arctan 2. 


Exercise: 


2 


1 
Problem: / : Z dy dx 
r+y? 
0 1 


In the following exercises, find the average value of the function over the given rectangles. 
Exercise: 


Problem: f(x,y) = —x + 2y, R = [0,1] x [0,1] 


Solution: 


i 

= 
Exercise: 

Problem: f(x,y) = z* + 2y’, R= [1,2] x [2,3] 
Exercise: 

Problem: f(x,y) = sinh x + sinh y, R= [0,1] x [0,2] 

Solution: 

+(2cosh 1 + cosh 2 — 3). 


Exercise: 


Problem: f(x,y) = arctan(zy), R = [0,1] x [0,1] 
Exercise: 


Problem: 


Let f and g be two continuous functions such that 0 < m, < f(x) < M; for any z € [a, bj and 
0 < m2 < g(y) < Mp for any y € [c, d]. Show that the following inequality is true: 


b d 
m m2(b — a)(c — d) < [| te)owdydr < MiMalo—a)(e-«), 


In the following exercises, use property v. of double integrals and the answer from the preceding exercise to show 
that the following inequalities are true. 
Exercise: 


Problem: 4 < [[-rra <1, where R = [0,1] x [0,1] 
R 
Exercise: 


Problem: x < [[smzcosyda < = where R = [4,4] x [ 


R 


|3 
wa 
ta 


Exercise: 

Problem: 0 < [[eveose aa < 4, where R = [0, +] x [0, +] 
Exercise: 

Problem: 0 < // (In z) (Iny)dA < (e — 1)”, where R = [1,e] x [1,e] 


Exercise: 


Problem: 


Let f and g be two continuous functions such that 0 < m, < f(x) < M; for any x € [a, 6] and 
0 < m2 < g(y) < Mp for any y € [c, d]. Show that the following inequality is true: 


(my t Mz3)(b o<| [ves y)|dy dx < < (M, + M))(b — a)(c— d). 


In the following exercises, use property v. of double integrals and the answer from the preceding exercise to show 
that the following inequalities are true. 
Exercise: 


Problem: 2 < // (e* + ev )dA < 2, where R = [0,1] x [0,1] 


Exercise: 
Problem: x < // (sinz+cosy)dA < v3 where R= [z, z] x [Z, q] 
Exercise: 


Problem: £e~7/? < // (cos z +e ¥)dA < 7, where R = [0,2] x [0,4] 


R 


Exercise: 


Problem: 4+ < ie —Inz)dA < 2, where R= [0,1] x [0,1] 


In the following exercises, the function fis given in terms of double integrals. 


a. Determine the explicit form of the function f. 

b. Find the volume of the solid under the surface z = f(a, y) and above the region R. 

c. Find the average value of the function f on R. 

d. Use a computer algebra system (CAS) to plot z = f(x, y) and z = faye in the same system of coordinates. 


Exercise: 
yY « 
Problem: [T] f(z, y) = i} [oes + yt)ds dt, where (x,y) € R= [0,1] x [0,1] 
0 0 


Solution: 


a 1 
a. f(z, y) = Fay (2? + y?) b.V = | [ tear ay = be toe Sa 
0 


0 


f(x, Y) = fies + yt)dsadt 


Exercise: 


zy 
Problem: [T] f(x, y) = / [tc0s(s) + cos (t)|dt ds, where (x, y) € R = [0,3] x [0,3] 
0 0 


Exercise: 


Problem: Show that if f and g are continuous on [a, 6] and [c, d], respectively, then 


/ fue) + g(y)|dy da = (d — ce) J rt) 


bod d d 0b 
+ / / gy) dy dx = (b — a) / oy) dy+ / / f(x)da dy. 


Exercise: 
Problem: 


b 


b d d 
Show tat [vt + xg(y)dydx = 5 (@? — ) | tae + 5 (0? — a’) [say 


a 


Exercise: 
Problem: [T] Consider the function f(x,y) = e~* ~¥, where (x, y) € R = [—1,1] x [—1,1]. 
a. Use the midpoint rule with m = n = 2, 4,..., 10 to estimate the double integral J = // eV dA. 


R 
Round your answers to the nearest hundredths. 


b. Form = n = 2, find the average value of f over the region R. Round your answer to the nearest 
hundredths. 


c. Use a CAS to graph in the same coordinate system the solid whose volume is given by | u eV dA 


and the plane z = fave- 


Solution: 


aForm=w=2,f=4e 22430. fc=e 20.61; 


c. 
f(x,y) =e*-¥ 
0.92 
0.72 
4 
0.52 
0.32 
0.12 * 
—-1.0 
Exercise: 


Problem: [T] Consider the function f(x,y) = sin (x”)cos (y”), where (x,y) € R = [-1,1] x [-1,1]. 
a. Use the midpoint rule with m = n = 2, 4,..., 10 to estimate the double integral 


L= i / sin (x”) cos (y’)dA. Round your answers to the nearest hundredths. 


b. Form = n = 2, find the average value of f over the region R. Round your answer to the nearest 
hundredths. 
c. Use a CAS to graph in the same coordinate system the solid whose volume is given by 


fl sin (x”)cos (y”)dA and the plane z = faye. 


In the following exercises, the functions f,, are given, where n > 1 is a natural number. 


a. Find the volume of the solids S,, under the surfaces z = f,(x, y) and above the region R. 
b. Determine the limit of the volumes of the solids S,, as n increases without bound. 


Exercise: 


Problem: f(x,y) = a" +y"+ zy, (x,y) € R= (0,1) x [0,1] 


Solution: 
2 1 1 
da. n+l + a b. 7 


Exercise: 


Problem: f(x,y) = <> + wr (,¥) € R= [1,2] x [1,2] 
Exercise: 


Problem: 
Show that the average value of a function f on a rectangular region R = [a,b] x [c,d] is 


me 2 * * * * 
Fave a > S> (xi, vi) , where (xij, Yi;) are the sample points of the partition of R, where 


i=1 j=l 
l<it<mandl<j<n. 
Exercise: 
Problem: 


Use the midpoint rule with m = n to show that the average value of a function fon a rectangular region 
R= [a,b] x [c,d] is approximated by 


Equation: 
12 1 1 
fave >of (Fen + @i), 5 (yja t u)). 
nm jj=l 
Exercise: 
Problem: 


An isotherm map is a chart connecting points having the same temperature at a given time for a given period 
of time. Use the preceding exercise and apply the midpoint rule with m = n = 2 to find the average 
temperature over the region given in the following figure. 


40 50 60 70 80 


Solution: 


56.5° F; here f(1,y;) = 71, f(e2,41) = 72, f(w2,41) = 40, f(x, 42) = 43, where «; and y; are the 
midpoints of the subintervals of the partitions of [a, b] and [c, d], respectively. 


Glossary 


double integral 
of the function f(, y) over the region R in the zy-plane is defined as the limit of a double Riemann sum, 


| / fle,y)dA = lim SOY Fein yisAA. 
R 


t=1 j=l 


double Riemann sum 
n 


f (x;,, y,;)AA where R is divided into smaller 
j=l 


Mas 


of the function f(, y) over a rectangular region R is 


7 


ll 
e 


* * 
subrectangles R;; and (x;;,y;;) is an arbitrary point in Rj; 
Fubini’s theorem 


if f(z, y) is a function of two variables that is continuous over a rectangular region 
R= { (2, y) € Rla< a <b,c< y< d}, then the double integral of f over the region equals an iterated 


integral, // f(x, y)dy dz = i [ f(x, y)dx dy = a iE f(x, y)dx dy 
R 


iterated integral 
for a function f(z, y) over the region R is 


af i: eae i , i “Hess de, 
b. f : [F Hless)ae dy = [ ; / “Sle, v)ae] dy, 


where a, b, c, and d are any real numbers and R = [a,b] x [c,d] 


Double Integrals over General Regions 


e Recognize when a function of two variables is integrable over a general region. 

e Evaluate a double integral by computing an iterated integral over a region bounded by two vertical lines and 
two functions of x, or two horizontal lines and two functions of y. 

e Simplify the calculation of an iterated integral by changing the order of integration. 

e Use double integrals to calculate the volume of a region between two surfaces or the area of a plane region. 

¢ Solve problems involving double improper integrals. 


In Double Integrals over Rectangular Regions, we studied the concept of double integrals and examined the tools 
needed to compute them. We learned techniques and properties to integrate functions of two variables over 
rectangular regions. We also discussed several applications, such as finding the volume bounded above by a 
function over a rectangular region, finding area by integration, and calculating the average value of a function of 
two variables. 


In this section we consider double integrals of functions defined over a general bounded region D on the plane. 
Most of the previous results hold in this situation as well, but some techniques need to be extended to cover this 
more general case. 


General Regions of Integration 
An example of a general bounded region D on a plane is shown in [link]. Since D is bounded on the plane, there 


must exist a rectangular region R on the same plane that encloses the region D, that is, a rectangular region R 
exists such that D is a subset of R(D C R). 


g(x,y) = fy) 


g(x, y) = 0 


For a region D that is a subset of R, we can 
define a function g(x, y) to equal f(z, y) at 
every point in D and 0 at every point of R not 
in D. 


Suppose z = f (z, y) is defined on a general planar bounded region D as in [link]. In order to develop double 
integrals of f over D, we extend the definition of the function to include all points on the rectangular region R and 
then use the concepts and tools from the preceding section. But how do we extend the definition of f to include all 
the points on R? We do this by defining a new function g (x, y) on Ras follows: 

Equation: 


f(z,y) if (z,y)isinD 
0 if (x, y) is in R but not in D 


(ew) ={ 


Note that we might have some technical difficulties if the boundary of D is complicated. So we assume the 
boundary to be a piecewise smooth and continuous simple closed curve. Also, since all the results developed in 
Double Integrals over Rectangular Regions used an integrable function f (a, y), we must be careful about g (2, y) 
and verify that g (x, y) is an integrable function over the rectangular region R. This happens as long as the region 
D is bounded by simple closed curves. For now we will concentrate on the descriptions of the regions rather than 
the function and extend our theory appropriately for integration. 


We consider two types of planar bounded regions. 


Note: 

Definition 

A region D in the (a, y)-plane is of Type I if it lies between two vertical lines and the graphs of two continuous 
functions g; (a) and go (a). That is ({link]), 

Equation: 


D={(z,y)la<z <b, (x) Sy < g (2)}. 
A region D in the zy plane is of Type II if it lies between two horizontal lines and the graphs of two continuous 
functions hy (y) and hg (y). That is ([link]), 
Equation: 


D={(z,y)\e< y Sd, hi (y) < & < ho (y)}- 


y = 9,(x) 


Y = go(x) 


y = 9,(x) eran 


o}---------- 


Crleeaaane 


1 
1 
1 
' 
1 
' 
1 

a 


A Type I region lies between two vertical lines and the graphs of two functions of z. 


x = h{y) 


A Type II region lies between two horizontal lines and the graphs of two 
functions of y. 


Example: 
Exercise: 


Problem: 
Describing a Region as Type I and Also as Type II 


Consider the region in the first quadrant between the functions y = ./z and y = z° ([link]). Describe the 
region first as Type I and then as Type II. 


0.2 0.4 0.6 0.8 1.0 1.2% 


Region D can be described as Type I or as Type II. 


Solution: 


When describing a region as Type I, we need to identify the function that lies above the region and the 
function that lies below the region. Here, region D is bounded above by y = ,/z and below by y = ee in the 
interval for x in [0, 1]. Hence, as Type I, D is described as the set {(x,y)|0 <2 <1,a?<y< yz}. 


However, when describing a region as Type II, we need to identify the function that lies on the left of the 
region and the function that lies on the right of the region. Here, the region D is bounded on the left by 


x = y* and on the right by x = </y in the interval for y in [0, 1]. Hence, as Type II, D is described as the set 
(ie pilav heya ae 35). 


Note: 
Exercise: 


Problem: 


Consider the region in the first quadrant between the functions y = 2x and y = x”. Describe the region first 
as Type I and then as Type II. 


Solution: 


Type I and Type II are expressed as Ae, yl0<2<2,2? <y< 2x} and 
{(z,y)|0<y<4, Fy <a < py}, respectively. 
Hint 


Graph the functions, and draw vertical and horizontal lines. 


Double Integrals over Nonrectangular Regions 


To develop the concept and tools for evaluation of a double integral over a general, nonrectangular region, we need 
to first understand the region and be able to express it as Type I or Type II or a combination of both. Without 
understanding the regions, we will not be able to decide the limits of integrations in double integrals. As a first 
step, let us look at the following theorem. 


Note: 

Double Integrals over Nonrectangular Regions 

Suppose g (x, y) is the extension to the rectangle R of the function f (x, y) defined on the regions D and R as 
shown in [link] inside R. Then g (2, y) is integrable and we define the double integral of f (2, y) over D by 


Equation: 
|[seuaa = |[oc y)dA. 
D R 


The right-hand side of this equation is what we have seen before, so this theorem is reasonable because FR is a 


rectangle and 7 / g (a, y)dA has been discussed in the preceding section. Also, the equality works because the 


values of g(x, y) are 0 for any point (2, y) that lies outside D, and hence these points do not add anything to the 
integral. However, it is important that the rectangle R contains the region D. 


As a matter of fact, if the region D is bounded by smooth curves on a plane and we are able to describe it as Type I 
or Type II or a mix of both, then we can use the following theorem and not have to find a rectangle R containing 
the region. 


Note: 
Fubini’s Theorem (Strong Form) 
For a function f (x, y) that is continuous on a region D of Type I, we have 


Equation: 
b g2(x 


|f e.naa= ff teaver = [ foo de, 
D D ) 


@ gi(% 


Similarly, for a function f (x, y) that is continuous on a region D of Type II, we have 


Equation: 
d| haf 
i] f(2,y)d4 = | f(z, y)de dy = | ; 


y 


) 
f(z, y)dz | dy. 


The integral in each of these expressions is an iterated integral, similar to those we have seen before. Notice that, 
in the inner integral in the first expression, we integrate f (x, y) with x being held constant and the limits of 
integration being g; (x) and gp («). In the inner integral in the second expression, we integrate f (x, y) with y 
being held constant and the limits of integration are hy (x) and hg (2). 


Example: 
Exercise: 


Problem: 
Evaluating an Iterated Integral over a Type I Region 


Evaluate the integral | / z?e7¥dA where D is shown in [link]. 


Solution: 


First construct the region D as a Type I region ([link]). Here D = Hee y)\0<2<2,g2<y< in Then 
we have 
Equation: 


We can express region D 


as a Type I region and 
integrate from y = $2 to 
y = 1, between the lines 
o = Oana = 2, 


Therefore, we have 


Equation: 
a gel n=2 y=1 
| [> *e%dydx = / / xe™dy| dx Iterated integral for a Type I region. 
z=0 yoig z=0 | y=1/22 
r=2 = Integrate with respect to y using 
= 2? = ' dx u-substitution with u = zy where z is held 
x ee constant. 


Integrate with respect to x using 


u-substitution with u = +2. 


((link]). 


This is a Type II region and the integral would then look like 
Equation: 


me xe%dax dy. 
=0 


However, if we integrate first with respect to x, this integral is lengthy to compute because we have to use 
integration by parts twice. 


Example: 
Exercise: 


Problem: 
Evaluating an Iterated Integral over a Type II Region 


Evaluate the integral // (3a? + y*)dA where = {(z,y)| -2<y<3,y2-3<a<y+3}. 
D 


Solution: 


Notice that D can be seen as either a Type I or a Type II region, as shown in [link]. However, in this case 
describing D as Type I is more complicated than describing it as Type II. Therefore, we use D as a Type II 
region for the integration. 
Type | region: Integrate Type II region: Integrate 
Yh first with respect to y Ya first with respect to x 


(a) (b) 


The region D in this example can be either (a) Type I or (b) Type II. 


Choosing this order of integration, we have 
Equation: 


y=3 z=y+3 
// (3a? + y)dA = i] Hl (az + y’)dx dy Iterated integral, 1] 
D 


y=—2 z=y?-3 


y=3 yt3 
= i (a? + ay’) dy Integrate with res} 
v= y2—-3 
y=3 
3 
= / ((y 3) ver dy (3) 3)y") dy 
ae 
3 
= i (54 + 27y — 12y? + 2y* + 8y* — y®)dy Integrate with resy 


[54y + ue Ay? vt aur | 


— 2375 
=“ 


Note: 
Exercise: 


Problem: 

Sketch the region D and evaluate the iterated integral | / xy dy dx where D is the region bounded by the 
D 

curves y = cos x and y = sin z in the interval |[—37/4, 7/4]. 


Solution: 


1/4 
Hint 


Express D as a Type I region, and integrate with respect to y first. 


Recall from Double Integrals over Rectangular Regions the properties of double integrals. As we have seen from 
the examples here, all these properties are also valid for a function defined on a nonrectangular bounded region on 
a plane. In particular, property 3 states: 


If R= SUT and SOT = 0 except at their boundaries, then 


Equation: 
[f1naa- [fteneas ff senee 


Similarly, we have the following property of double integrals over a nonrectangular bounded region on a plane. 


Note: 
Decomposing Regions into Smaller Regions 
Suppose the region D can be expressed as D = D; U D2 where D; and Dy do not overlap except at their 


boundaries. Then 
I ft e.w)aa = i / f (0,y)dA + i / f(o,y)aa. 


Equation: 


This theorem is particularly useful for nonrectangular regions because it allows us to split a region into a union of 
regions of Type I and Type II. Then we can compute the double integral on each piece in a convenient way, as in 
the next example. 


Example: 
Exercise: 


Problem: 
Decomposing Regions 


Express the region D shown in [link] as a union of regions of Type I or Type II, and evaluate the integral 
Equation: 


i / (22 + 5y)dA. 


This region can be decomposed into 
a union of three regions of Type I or 
Type Il. 


Solution: 


The region D is not easy to decompose into any one type; it is actually a combination of different types. So 
we can write it as a union of three regions D;, D2, and D3 where, 


Dy = {(z,»)| =23<n7<0,0<7< (x +2)?}, D = eae NOES SEs (y — iy*)}. 


These regions are illustrated more clearly in [link]. 


Breaking the region into three 


subregions makes it easier to set up 
the integration. 


Here Dj is Type I and D2 and D3 are both of Type II. Hence, 


Equation: 
[fc +5y)dA = fees 5y)dA + ffee+ 5y)dA + [fees 5y)dA 
D D, D2 D3 
2=0 y=(2+2)” y=4 x=y—(1/16)y° y=0 v=y—(1/16)y? 
= i / (2¢ + 5y)dy dx + i (2 + 5y)da dy + 
z=-2 y=0 y=0 z=0 y=—4 z=—-2 
z=0 y=4 
= —(2+2)?(20 + 24a + 5x?)| 4 / z y® if y' + 6y? 
256 16 
xz=—2 y=0 
apt 7 
6 4 2 
+ 6y~ + 10 4 
a / aaa fe oe | 
y=—4 
= 4 ag = 


Now we could redo this example using a union of two Type II regions (see the Checkpoint). 


Note: 
Exercise: 


Problem: 


Consider the region bounded by the curves y = In z and y = e” in the interval [1, 2]. Decompose the region 
into smaller regions of Type II. 


Solution: 


ei g =a re MUNG ae ae NUN ile =agice mye oe) 


Hint 


Sketch the region, and split it into three regions to set it up. 


Note: 
Exercise: 


Problem: Redo [link] using a union of two Type II regions. 


Solution: 


Same as in the example shown. 


Hint 


(eps ys tot yg = es (9 > vr) PU ey) Fs a0 2 eee GY} 


Changing the Order of Integration 


As we have already seen when we evaluate an iterated integral, sometimes one order of integration leads to a 
computation that is significantly simpler than the other order of integration. Sometimes the order of integration 
does not matter, but it is important to learn to recognize when a change in order will simplify our work. 


Example: 
Exercise: 


Problem: 
Changing the Order of Integration 


Reverse the order of integration in the iterated integral / ze* dy dx. Then evaluate the new 
=() 


iterated integral. 


Solution: 


The region as presented is of Type I. To reverse the order of integration, we must first express the region as 
Type II. Refer to [link]. 


y 
(0, 2) (0, 2) 
y=2-x2 x= \2-y 
(0, 0) (2,0) = * — (0,0) (2,0) * 


Type | Type Il 


Converting a region from Type I to Type I. 


We can see from the limits of integration that the region is bounded above by y = 2 — x? and below by 
y = 0, where z is in the interval lo, v3] . By reversing the order, we have the region bounded on the left by 


az = 0 and on the right by 2 = ,/2 — y where y is in the interval [0, 2]. We solved y = 2 — x? in terms of x 
to obtain x = /2 =p 


Hence 
Equation: 


V2 2-2? V2-y Reverse the orc 


2 
/ i ze* dydz = / / ae® dx dy integration the: 
0 


0 substitution. 


2 
L = i 
Feat 2 j dy = [ye — 1)dy = —1(e?-¥ +y) |, 


II 
SI= See 


Example: 
Exercise: 


Problem: 
Evaluating an Iterated Integral by Reversing the Order of Integration 


Consider the iterated integral | / f (x, y)dx dy where z = f (x, y) = x — 2y over a triangular region R that 
R 


has sides on z = 0, y = 0, and the line x + y = 1. Sketch the region, and then evaluate the iterated integral 
by 


a. integrating first with respect to y and then 
b. integrating first with respect to x. 
Solution: 


A sketch of the region appears in [link]. 


A triangular region R for 
integrating in two ways. 


We can complete this integration in two different ways. 


a. One way to look at it is by first integrating y from y = 0 to y = 1 — z vertically and then integrating x 
ional 4p = (O)inaae = Ils 
Equation: 


(x — 2y)dy dx = fi ry — 2y?| [ace 
2=0 


[[seaaedy = ff 
R =(1) 


= lr eo i [-14 30 -20"]de = |-2+ Se -= 
=0 


b. The other way to do this problem is by first integrating x from x = 0 tox = 1 — y horizontally and 
then integrating y from y = Oto y = 1: 


Equation: 
yee — y— 1 pees 
i j@wecy = / / (a — 2y)dz dy = / Ea — 20) dy 
R y=0 =O y=0 ee 
yl v— 


Note: 
Exercise: 


Problem: 
Evaluate the iterated integral | i) (a 4F y)dA over the region D in the first quadrant between the functions 


y = 2x and y = x”. Evaluate the iterated integral by integrating first with respect to y and then integrating 
first with resect to x. 
Solution: 


216 
35 


Hint 


Sketch the region and follow [link]. 


Calculating Volumes, Areas, and Average Values 


We can use double integrals over general regions to compute volumes, areas, and average values. The methods are 
the same as those in Double Integrals over Rectangular Regions, but without the restriction to a rectangular region, 


we can now solve a wider variety of problems. 


Example: 
Exercise: 


Problem: 
Finding the Volume of a Tetrahedron 


Find the volume of the solid bounded by the planes x = 0, y = 0, z = 0, and 2x + 3y+ z=6. 
Solution: 


The solid is a tetrahedron with the base on the zy-plane and a height z = 6 — 2” — 3y. The base is the 
region D bounded by the lines, « = 0, y = 0 and 2x + 3y = 6 where z = 0 ([link]). Note that we can 


consider the region D as Type I or as Type II, and we can integrate in both ways. 
z=6-2x—- 3y 


x 


A tetrahedron consisting of the three coordinate planes and the plane z = 6 — 2a — 3y, with the 
base bound by z = 0, y = 0, and 2% + 3y = 6. 


First, consider D as a Type I region, and hence D = {(z,y)|0 <a <3,0<y<2- a}. 


Therefore, the volume is 


Equation: 
z=3 y=2—(2x/3) x=3 3 y=2—(22/3) 
Va— (6 — 2x — 3y)dy dr = / (sy — 2ry — 3”) dx 
nt y=0 z=0 - 
2=3 9 
z=0 


Now consider D as a Type II region, so D = {(x,y)|0 < y < 2,0 < x <3 — 3y}. In this calculation, the 
volume is 
Equation: 


I 
| 
— ~~ 
< 
| 
i) 
na 
bo 
| cr | 
Q 
K— 
I 
@ 


y=0 


Therefore, the volume is 6 cubic units. 


Note: 
Exercise: 


Problem: 


Find the volume of the solid bounded above by f (x, y) = 10 — 2a + y over the region enclosed by the 
curves y = 0 and y = e”, where z is in the interval [0, 1]. 


Solution: 


< + 10e — 2 cubic units 
Hint 


Sketch the region, and describe it as Type I. 


Finding the area of a rectangular region is easy, but finding the area of a nonrectangular region is not so easy. As 
we have seen, we can use double integrals to find a rectangular area. As a matter of fact, this comes in very handy 
for finding the area of a general nonrectangular region, as stated in the next definition. 


Note: 

Definition 

The area of a plane-bounded region D is defined as the double integral | / 1dA. 
D 


We have already seen how to find areas in terms of single integration. Here we are seeing another way of finding 
areas by using double integrals, which can be very useful, as we will see in the later sections of this chapter. 


Example: 
Exercise: 


Problem: 
Finding the Area of a Region 


Find the area of the region bounded below by the curve y = x” and above by the line y = 2z in the first 
quadrant ([link]). 


The region bounded by 
y = x? and y = 22. 


Solution: 


We just have to integrate the constant function f (2, y) = 1 over the region. Thus, the area A of the bounded 


a2 y=2e eal Ba 
region is i i dy dz or / / dx dy: 
C=) YE y=0  2=y/2 
Equation: 
e=2 y=2e mth DES a 
A [[azay- / ldydz = / [vit Jae = / (2a = x”) de oe fy] =. 
. aioe a #=0 2=0 0 
Note: 
Exercise: 
Problem: 


Find the area of a region bounded above by the curve y = x? and below by y = 0 over the interval (0, 3]. 


Solution: 


81 P 
7 Square units 
Hint 


Sketch the region. 


We can also use a double integral to find the average value of a function over a general region. The definition is a 
direct extension of the earlier formula. 


Note: 

Definition 

If f (x, y) is integrable over a plane-bounded region D with positive area A (D), then the average value of the 
function is 

Equation: 


fave = ap If fenaa 


Note that the area is A (D) = // 1dA. 
D 


Example: 
Exercise: 


Problem: 
Finding an Average Value 


Find the average value of the function f (x, y) = 7xy? on the region bounded by the line x = y and the 


curve x = VY ([link]). 
y (1, 1) 


0.75 
0.5 


0.25 


The region bounded by x = y and 


z= /Y. 


Solution: 


First find the area A (D) where the region D is given by the figure. We have 
Equation: 


y=l t=V/y y= y=1 F 


A(D)= ff1a= ff ravay= f [eyav= f Wa-way= 50°” ler 


y=0 w=y y=0 y=0 


Then the average value of the given function over this region is 


Equation: 
y=l tay y=1 
if = aty ff Ftc )dA = : | teptae dy = : male 
ave —_ A(D) Y RK A(D y Y= 1 6 ) y Bet Yy 
D y=0 «r=y y=0 
y=1 y=1 4 
7 7 2 (4A 5 
= [Ev (u-¥")]av=6 EW -v) v= 2 (4-4)|- 8-3 
y=0 y=0 
Note: 
Exercise: 
Problem: 


Find the average value of the function f (x, y) = xy over the triangle with vertices (0,0), (1,0) and (1, 3). 


Solution: 


3 


4 
Hint 


Express the line joining (0,0) and (1,3) asa function y = g(z). 


Improper Double Integrals 


An improper double integral is an integral | / f dA where either D is an unbounded region or f is an 
D 


unbounded function. For example, D = {(«, y)| |z — y| > 2} is an unbounded region, and the function 
f (x,y) =1/ (1 — x? — 2y?) over the ellipse x? + 3y? < 1 is an unbounded function. Hence, both of the 


following integrals are improper integrals: 


i, [[2y4A where D = {(a, y)| |x — y| = 2}; 


D 
1 
iff Tog? cage if where D = {(2,u)la? + 39° < 1}. 


In this section we would like to deal with improper integrals of functions over rectangles or simple regions such 
that f has only finitely many discontinuities. Not all such improper integrals can be evaluated; however, a form of 
Fubini’s theorem does apply for some types of improper integrals. 


Note: 

Fubini’s Theorem for Improper Integrals 

If D is a bounded rectangle or simple region in the plane defined by {(z, y):a < « < b, g(x) < y < h(x)} and 
also by {(z, y):¢ < y< d,j(y) < x < k(y)} and f is a nonnegative function on D with finitely many 
discontinuities in the interior of D, then 


Equation: 


z=b y=h(x) y=d x=k(y) 
[fsa = f (2, y)dy dx = / / f (x, y)dz dy. 
D T=A y—g(z) y=c  2=j(y) 


It is very important to note that we required that the function be nonnegative on D for the theorem to work. We 
consider only the case where the function has finitely many discontinuities inside D. 


Example: 
Exercise: 


Problem: 
Evaluating a Double Improper Integral 


Consider the function f (x, y) = a over the region D = {(z,y):0 <a<le<y< vz}. 


Notice that the function is nonnegative and continuous at all points on D except (0,0). Use Fubini’s theorem 
to evaluate the improper integral. 


Solution: 


First we plot the region D ([link]); then we express it in another way. 


The function f is continuous at all points of the 
region D except (0, 0). 


The other way to express the same region D is 
Equation: 


DSN Gy) eee 


Thus we can use Fubini’s theorem for improper integrals and evaluate the integral as 
Equation: 


y 
‘i i © dx dy 
y 
y=0 a=y? 
Therefore, we have 
Equation: 
gl B=7 ; vi ; y=1 1 
e€ e = ey 
fl / yo if yl = i ra (y — y?)dy = fe — ye" dy = e — 2. 
y—0 a=y? y=0 y=0 0 


As mentioned before, we also have an improper integral if the region of integration is unbounded. Suppose now 
that the function f is continuous in an unbounded rectangle R. 


Note: 
Improper Integrals on an Unbounded Region 


If R is an unbounded rectangle such as R = {(x,y):a < @ < co,c < y < oo}, then when the limit exists, we 
d 


b d b 
have [fsenes = rola) | [tendy dx = salt. | [tendy dy. 


@ 


The following example shows how this theorem can be used in certain cases of improper integrals. 


Example: 
Exercise: 


Problem: 
Evaluating a Double Improper Integral 


Evaluate the integral | / aye” dA where R is the first quadrant of the plane. 


Solution: 


The region R is the first quadrant of the plane, which is unbounded. So 
Equation: 


az=b 
Sy = : Seay, = : =afiaap 
// rte ees Hl gee Dacian axe / / ES ed ae 
R 0 


= pee (1 = et) (1 i e*) a 


Thus, | I aye * ¥ dA is convergent and the value is + 


R 


Note: 
Exercise: 


Problem: 


Evaluate the improper integral // Fo where D = {(a,y)z >0,y > 0,27 +y? < 1}. 
1—2?—-y? 
D 


Solution: 


Tr 


4 
Hint 


Notice that the integral is nonnegative and discontinuous on x? + y? = 1. Express the region D as 
D= {(2, y:0<2<10<y<vV1-— a} and integrate using the method of substitution. 


In some situations in probability theory, we can gain insight into a problem when we are able to use double 
integrals over general regions. Before we go over an example with a double integral, we need to set a few 
definitions and become familiar with some important properties. 


Note: 

Definition 

Consider a pair of continuous random variables X and Y, such as the birthdays of two people or the number of 
sunny and rainy days in a month. The joint density function f of X and Y satisfies the probability that (X,Y) lies 
in a certain region D: 

Equation: 


P((x,¥)€D)= ff fa,yjaa. 
D 


Since the probabilities can never be negative and must lie between 0 and 1, the joint density function satisfies the 
following inequality and equation: 
Equation: 


f (x,y) 2 Oand |[ #@waa- 1. 
R2 


Note: 

Definition 

The variables X and Y are said to be independent random variables if their joint density function is the product of 
their individual density functions: 

Equation: 


f (x,y) = fi (@) fe (y). 


Example: 
Exercise: 


Problem: 
Application to Probability 


At Sydney’s Restaurant, customers must wait an average of 15 minutes for a table. From the time they are 
seated until they have finished their meal requires an additional 40 minutes, on average. What is the 
probability that a customer spends less than an hour and a half at the diner, assuming that waiting for a table 
and completing the meal are independent events? 


Solution: 


Waiting times are mathematically modeled by exponential density functions, with m being the average 
waiting time, as 
Equation: 


; 0 ift <0, 
FO) 1 e-tim ift > 0. 


If X and Y are random variables for ‘waiting for a table’ and ‘completing the meal,’ then the probability 
density functions are, respectively, 
Equation: 


0 f2<0, 0 ify < 0, 
fi(z) = de-a/i5 fe 0 and fo(y) = qe 4/40 ify > 0. 


Clearly, the events are independent and hence the joint density function is the product of the individual 
functions 
Equation: 


0 ifx < Oor 0, 
fees) = niente) =| <UVory< 


1_,-2/15,-y/60 ; 
s00© /lbe-v/60 if, y > 0. 


We want to find the probability that the combined time X + Y is less than 90 minutes. In terms of geometry, 
it means that the region D is in the first quadrant bounded by the line + y = 90 ([link]). 


y 
100 


The region of integration for a joint probability 
density function. 


Hence, the probability that (X, Y) is in the region D is 
Equation: 


P(X+Y < 90) = P((X,Y)€D)= |[#enaa = If aoote waa 
D D 


Since z + y = 90 is the same as y = 90 — a, we have a region of Type I, so 


Equation: 
D = {(a,y)|0 <2 <90,0< y< 90-3}, 
z=90 y=90—-z z=90 y=90—x 
P(X+Y<90) = & / / e 7/6 -¥/M de dy = sae / / e 7/5 ——9/M da dy 

z=0 y=0 a=0 y=0 

2=90 y=90—x 
= / / eW (2/15 +4/4) dee dy = 0.8328. 
z=—0 y=0 


Thus, there is an 83.2% chance that a customer spends less than an hour and a half at the restaurant. 


Another important application in probability that can involve improper double integrals is the calculation of 
expected values. First we define this concept and then show an example of a calculation. 


Note: 

Definition 

In probability theory, we denote the expected values E (X) and E (Y), respectively, as the most likely outcomes 
of the events. The expected values EF (X) and E (Y) are given by 


Equation: 
E(x) = [fe (z,y)dA and E (Y )= fos (x, y)dA, 


where S is the sample space of the random variables X and Y. 


Example: 
Exercise: 


Problem: 
Finding Expected Value 


Find the expected time for the events ‘waiting for a table’ and ‘completing the meal’ in [link]. 
Solution: 


Using the first quadrant of the rectangular coordinate plane as the sample space, we have improper integrals 
for E (X) and E'(Y). The expected time for a table is 
Equation: 


xL=CO 


4 q=es 
_ 7 -2/15,,—y/40 merle e 2/15 —y/40 

|[- 600° e€ dA= sf fe e€ dA 

S x=0 =0 

= BOO Rete ni f ze */e—¥/Wdz dy 


I 


E(X) 


=a y=b 
2 : —y/40 
= lim re */ da lim y/40 q 
sot a= 5-00 is 
z=0 y=0 


y=b 
y=0 


= sh ( (tim, (-15e-*"* (e +15))) 7”) ( (2 (-40e-¥) 


= gop ( lim (—15e74/9 (@ + 15) + 225)) ( jim (—40¢~°/4° +40) 
el 

Se 0) 

=, 


A similar calculation shows that E (Y) = 40. This means that the expected values of the two random events 
are the average waiting time and the average dining time, respectively. 


Note: 
Exercise: 


Problem: The joint density function for two random variables X and Y is given by 
Equation: 


an (2? +y") if0< 2 <15,0<y<10 


f(z,y) = 2 


0 otherwise 


Find the probability that X is at most 10 and Y is at least 5. 
Solution: 


2 = 0.7638 
Hint 


10 y=10 


1 
Compute the probability P(X <10,Y >5)= 6000 (x? + y*)dy dz. 


Z=—o0o y=5 


Key Concepts 


¢ A general bounded region D on the plane is a region that can be enclosed inside a rectangular region. We can 
use this idea to define a double integral over a general bounded region. 

e To evaluate an iterated integral of a function over a general nonrectangular region, we sketch the region and 
express it as a Type I or as a Type II region or as a union of several Type I or Type II regions that overlap only 
on their boundaries. 

e We can use double integrals to find volumes, areas, and average values of a function over general regions, 
similarly to calculations over rectangular regions. 

e We can use Fubini’s theorem for improper integrals to evaluate some types of improper integrals. 


Key Equations 
e Iterated integral over a Type I region 
b 92(2) 
// f(z, y)dA = // f(x, y)dy dx = / / F(a, y)dy | dx 
D D a | g(x) 
e Iterated integral over a Type II region 
d | ha(y) 
|[tenaa= ff tenaeay= [| f tewae| ay 
D D ¢c hi(y) 


In the following exercises, specify whether the region is of Type I or Type II. 
Exercise: 


Problem: The region D bounded by y = x°, y = x? +1, x = 0, and x = 1 as given in the following figure. 


-02 % 02 04 06 08 1.0 1.2% 


Exercise: 


Problem: 
Find the average value of the function f (x, y) = 3zy on the region graphed in the previous exercise. 


Solution: 


27 
20 


Exercise: 


Problem: Find the area of the region D given in the previous exercise. 
Exercise: 


Problem: 


The region D bounded by y = sinz,y = 1+sinz,x =0, and z = as given in the following figure. 


28 


y=1+sin(x) 


y = sin(x) 


Solution: 


Type I but not Type I 
Exercise: 


Problem: 


Find the average value of the function f (x, y) = cos x on the region graphed in the previous exercise 
Exercise: 


Problem: Find the area of the region D given in the previous exercise. 


Solution: 


T 


2 


Exercise: 


Problem: 


The region D bounded by x = y? — 1 and x = “, 1 — y? as given in the following figure. 


Exercise: 


Problem: 


Find the volume of the solid under the graph of the function f (2, y) = xy + 1 and above the region in the 
figure in the previous exercise. 


Solution: 


= (8 + 37) 


Exercise: 


Problem: The region D bounded by y = 0,2 = —10+ y, and x = 10 — yas given in the following figure. 


Exercise: 


Problem: 


Find the volume of the solid under the graph of the function f (x, y) = x + y and above the region in the 
figure from the previous exercise. 


Solution: 


1000 
3 


Exercise: 


Problem: The region D bounded by y = 0,z = y— 1, x = 4 as given in the following figure. 


Exercise: 


Problem: The region D bounded by y = 0 and y = x? — 1 as given in the following figure. 


Solution: 


Type I and Type II 
Exercise: 


Problem: 


Let D be the region bounded by the curves of equations y = x, y = —a, and y = 2 — x”. Explain why D is 
neither of Type I nor I. 


Exercise: 


Problem: 


Let D be the region bounded by the curves of equations y = cos x and y = 4 — z? and the z-axis. Explain 
why D is neither of Type I nor II. 


Solution: 


The region D is not of Type I: it does not lie between two vertical lines and the graphs of two continuous 
functions gi(x) and g(x). The region D is not of Type II: it does not lie between two horizontal lines and 
the graphs of two continuous functions hi(y) and h2(y). 


In the following exercises, evaluate the double integral | / f (x, y)dA over the region D. 
D 
Exercise: 


Problem: f (x,y) = 2 + 5yand D = {(2,y)|(0<2<1,a><y<a'+1} 


Exercise: 
Problem: f (x,y) = land D= {(z,y)|0 <a < 4,sma <y<1+sinz} 


Solution: 


TT 


2 
Exercise: 


Problem: f (x,y) = 2 and D= {(z,y)|0< y<1,y—1< <2 <arccos y} 


Exercise: 
Problem: f (x,y) = zy and D = {(z,y)| -l<y<ly-l<a¢< VJ1 =} 


Solution: 
0 


Exercise: 


Problem: f (x,y) = sin y and D is the triangular region with vertices (0,0), (0,3), and (3, 0) 


Exercise: 
Problem: f (x,y) = —x + 1 and Dis the triangular region with vertices (0,0), (0,2), and (2, 2) 
Solution: 
2 
3: 


Evaluate the iterated integrals. 
Exercise: 


1 32 
Problem: / i; (a + y’)dy dx 
0 22x 


Exercise: 


1 27x41 
Problem: / i (xy + 1)dy dx 
0 2/z 


Solution: 


AL 
20 


Exercise: 


e? 2 


Problem: / forts u)dv du 


e Inu 


Exercise: 


2 <1 
Problem: / / (8uv)dv du 
Ll -u?-1 


Solution: 


—63 


Exercise: 


Vi-¥ 
Problem: / LL (2a + 4a*) dx dy 


Exercise: 
1/2 /1—4y? 
Problem: / / Adz dy 
0 —~/1—4yP? 
Solution: 
T 
Exercise: 


Problem: Let D be the region bounded by y = 1 — x?, y = 4 — 2?, and the z- and y-axes. 


1 4-2? 2 4-2? 
a. Show that | / zdA = / / x dydx + / / x dy dx by dividing the region D into two regions 
0 1-2? 0 
of Type I. 
b. Evaluate the integral | / xdA. 
D 
Exercise: 


Problem: Let D be the region bounded by y = 1, y = a, y = Ing, and the z-axis. 


1 2 e 
a. Show that |[v dA= / iE dy dx + ' iE dy dz by dividing D into two regions of Type I. 
D 0 0 1 1 


b. Evaluate the integral | fl y dA. 


Solution: 


Be 2 
a. Answers may vary; b. + 


Exercise: 
Problem: 
0 2-2? 1 2-2? 
a. Show that i IE 21h = i / y’dy dx + / / y’ dy dz by dividing the region D into two regions 
-1 -2 


of Type I, mee {(x,y)|y > 2, roe oe ae “y. 


b. Evaluate the integral il / yA. 


Exercise: 


Problem: Let D be the region bounded by y = x”, y = x + 2, andy = —2. 


1 vy 2 Vv¥ 
a. Show that | / zdA= / / x dx dy + / i x dx dy by dividing the region D into two regions of 
D 0 -y 1 y-2 


Type I, where D = {(x,y)|y >a ,y>-a,y<a2t+ 2}. 


b. Evaluate the integral | / adA. 


D 


Solution: 


a. Answers may vary; b. ~ 
Exercise: 
Problem: 


The region D bounded by z = 0, y = x° +1, and y = 3 — 2? is shown in the following figure. Find the area 
A (D) of the region D. 


yi Y=3-2% 


-10 0 02 04 06 08 10 1.2x 


Exercise: 


Problem: 


The region D bounded by y = cos a, y = 4cosa, andz = +4 is shown in the following figure. Find the 
area A (D) of the region D. 


y=4+.cosx 


-15-1.0-0.5 9 
=i 


Solution: 


8r 
3 


Exercise: 


Problem: Find the area A(D) of the region D = {(z,y)|y>1—2?,y<4-—2?,y>0,a > 0}. 
Exercise: 


Problem: 


Let D be the region bounded by y = 1, y = a, y = Ing, and the z-axis. Find the area A (D) of the region D 


Solution: 


e— 


bo|oo 


Exercise: 


Problem: 


Find the average value of the function f (x, y) = sin y on the triangular region with vertices (0, 0), (0, 3), 
and (3, 0). 

Exercise: 
Problem: 


Find the average value of the function f (x, y) = 
and (2,2). 


—a + 1 on the triangular region with vertices (0,0), (0, 2), 
Solution: 


2. 
3 


In the following exercises, change the order of integration and evaluate the integral. 
Exercise: 


m/2 atl 
Problem: / sin x dy dx 
-1 0 
Exercise: 


~~ 

5 

i=} 

i=" 

o 

Z 
ai 


Solution: 

1 1-2 0 ytl 1 ly 
/ [eaae= ff cavays [ x dxdy = + 
0 «1 -1 0 0 0 

Exercise: 
0 Vytl 
Problem: a / y’dx dy 
“1 _ gi 
Exercise: 
1/2 VyP+i 
Problem: / / y dx dy 
-1/2 _/Pqi 
Solution: 

1/2 f+ 2 vei 

i / yr dy= [ | ydydz = 0 
1/2 fy Lo -Vz?=1 

Exercise: 
Problem: 


The region D is shown in the following figure. Evaluate the double integral | / (x? + y) dA by using the 


easier order of integration. 


Exercise: 


Problem: 


The region D is given in the following figure. Evaluate the double integral | / (a? a y’)dA by using the 
D 


easier order of integration. 


Solution: 
1 1-y! 
// (x? —y*)dA= / (x? — y?)\ dx dy = ayia 
D —1 yf-1 
Exercise: 
Problem: 


Find the volume of the solid under the surface z = 2x + y? and above the region bounded by y = x° and 
y= a. 

Exercise: 
Problem: 


Find the volume of the solid under the plane z = 3a + y and above the region determined by y = x” and 
y= a. 


Solution: 

4 

5 
Exercise: 


Problem: 


Find the volume of the solid under the plane z = x — y and above the region bounded by 
z=tany,x = —tany, andz = 1. 

Exercise: 
Problem: 


Find the volume of the solid under the surface z = x? and above the plane region bounded by 
z=siny,xz = —siny,andz = 1. 


Solution: 


Sn 
32 


Exercise: 
Problem: 


Let g be a positive, increasing, and differentiable function on the interval [a, b]. Show that the volume of the 
solid under the surface z = g/(x) and above the region bounded by y = 0, y= g(x), x = a, and x = bis 


given by + (9 (b) — g? (a)). 
Exercise: 
Problem: 
Let g be a positive, increasing, and differentiable function on the interval a, b], and let k be a positive real 


number. Show that the volume of the solid under the surface z = g/(x) and above the region bounded by 
y= 9(x),y=9(x) +k, x =a, and x = bis given by k(g(b) — g(a). 


Exercise: 
Problem: 
Find the volume of the solid situated in the first octant and determined by the planes z = 2, 
z=0,2+y=1,2 =0, andy=0. 
Exercise: 
Problem: 


Find the volume of the solid situated in the first octant and bounded by the planes x + 2y = 1, 
z=0,y=0,z=4, andz=0. 


Solution: 


1 
Exercise: 


Problem: 


Find the volume of the solid bounded by the planes zx + y= 1,2 —-y=1,x =0,z=0,andz=10. 
Exercise: 
Problem: 


Find the volume of the solid bounded by the planesza@+y=1,2—y=1l1,r+y=-—l, 
x—-y=-—l,z=landz=0. 


Solution: 


2 
Exercise: 
Problem: 


Let ; and S¥ be the solids situated in the first octant under the planes x + y+ z= 1landx+y+2z=1, 
respectively, and let S be the solid situated between S;, 52,2 = 0, andy = 0. 


a. Find the volume of the solid S$}. 
b. Find the volume of the solid So. 
c. Find the volume of the solid S' by subtracting the volumes of the solids S; and S9. 


Exercise: 
Problem: 


Let S and S} be the solids situated in the first octant under the planes 2x + 2y+ z2=2andxr+y+z=1, 
respectively, and let S be the solid situated between S;,.S2,2 = 0, andy =0. 


a. Find the volume of the solid S$}. 
b. Find the volume of the solid $5. 
c. Find the volume of the solid S by subtracting the volumes of the solids S; and $3. 


Solution: 


Exercise: 


Problem: 


Let S; and S¥ be the solids situated in the first octant under the plane x + y+ z = 2 and under the sphere 
a? +y? + 2* =4, respectively. If the volume of the solid S'y is =, determine the volume of the solid S 
situated between S$; and S»2 by subtracting the volumes of these solids. 


Exercise: 
Problem: 


Let S; and S‘ be the solids situated in the first octant under the plane x + y+ z = 2 and bounded by the 
cylinder x? + y? = 4, respectively. 


a. Find the volume of the solid S$}. 
b. Find the volume of the solid S». 
c. Find the volume of the solid S' situated between S; and S>2 by subtracting the volumes of the solids S; 


and S. 
Solution: 
a. 4; b. 273 c. am 
Exercise: 
Problem: 


[T] The following figure shows the region D bounded by the curves y = sin z, z = 0, and y = x*. Usea 
graphing calculator or CAS to find the x-coordinates of the intersection points of the curves and to determine 
the area of the region D. Round your answers to six decimal places. 


Exercise: 
Problem: 
[T] The region D bounded by the curves y = cos x, z = 0, and y = 2? is shown in the following figure. 


Use a graphing calculator or CAS to find the x-coordinates of the intersection points of the curves and to 
determine the area of the region D. Round your answers to six decimal places. 


Solution: 


0 and 0.865474; A (D) = 0.621135 
Exercise: 
Problem: 
Suppose that (X, Y) is the outcome of an experiment that must occur in a particular region S in the zy-plane. 


In this context, the region S is called the sample space of the experiment and X and Y are random variables. 
If D is a region included in S, then the probability of (X, Y) being in D is defined as 


P(X,Y)eD = | / p(x, y)dx dy, where p(z, y) is the joint probability density of the experiment. Here, 


p(, y) is a nonnegative function for which | i p(x, y)dx dy = 1. Assume that a point (X,Y) is chosen 
Ss 


arbitrarily in the square [0,3] x [0,3] with the probability density 


p(x, y) = . (x,y) € [0,3] x [0,3], 


0 otherwise. 


Find the probability that the point (X, Y) is inside the unit square and interpret the result. 
Exercise: 


Problem: 


Consider X and Y two random variables of probability densities p, (x) and p2(z), respectively. The random 
variables X and Y are said to be independent if their joint density function is given by 

p(«, y) = pi(x)p2(y). At a drive-thru restaurant, customers spend, on average, 3 minutes placing their orders 
and an additional 5 minutes paying for and picking up their meals. Assume that placing the order and paying 
for/picking up the meal are two independent events X and Y. If the waiting times are modeled by the 
exponential probability densities 


1 e-2/3 >0 1Le-y/5 >0 
3 e @&7 Vv, a 5 e y = 0, 
L).= and = 
pi(z) i: otherwise, aly) . otherwise, 


respectively, the probability that a customer will spend less than 6 minutes in the drive-thru line is given by 
P[IX+Y<6)/= [[e(eu)ae dy, where D = {(z,y)}|a >0,y>0,2+y <6}. Find P[X+Y <6] 


and interpret the result. 
Solution: 


P[X+Y<6)/=1+ ot - _ = 0.45; there is a 45% chance that a customer will spend 6 minutes in the 
drive-thru line. 


Exercise: 
Problem: 


[T] The Reuleaux triangle consists of an equilateral triangle and three regions, each of them bounded by a 
side of the triangle and an arc of a circle of radius s centered at the opposite vertex of the triangle. Show that 


the area of the Reuleaux triangle in the following figure of side length s is = (x - v3) : 


Exercise: 
Problem: 
[T] Show that the area of the lunes of Alhazen, the two blue lunes in the following figure, is the same as the 


area of the right triangle ABC. The outer boundaries of the lunes are semicircles of diameters AB and AC, 
respectively, and the inner boundaries are formed by the circumcircle of the triangle ABC’. 


Glossary 


improper double integral 
a double integral over an unbounded region or of an unbounded function 


Type I 
a region D in the xy-plane is Type I if it lies between two vertical lines and the graphs of two continuous 
functions g; (x) and go (x) 


Type I 
a region D in the xy-plane is Type II if it lies between two horizontal lines and the graphs of two continuous 
functions hy (y) and hg (y) 


Double Integrals in Polar Coordinates 


e Recognize the format of a double integral over a polar rectangular region. 

e Evaluate a double integral in polar coordinates by using an iterated integral. 
e Recognize the format of a double integral over a general polar region. 

e Use double integrals in polar coordinates to calculate areas and volumes. 


Double integrals are sometimes much easier to evaluate if we change rectangular coordinates to polar 
coordinates. However, before we describe how to make this change, we need to establish the concept of a 
double integral in a polar rectangular region. 


Polar Rectangular Regions of Integration 


When we defined the double integral for a continuous function in rectangular coordinates—say, g over a 
region R in the zy-plane—we divided R into subrectangles with sides parallel to the coordinate axes. 
These sides have either constant z-values and/or constant y-values. In polar coordinates, the shape we work 
with is a polar rectangle, whose sides have constant r-values and/or constant 6-values. This means we can 
describe a polar rectangle as in [link](a), with R = {(r, 0)|a<r<ba<6< ph. 


In this section, we are looking to integrate over polar rectangles. Consider a function f (r, @) over a polar 
rectangle R. We divide the interval {a, b] into m subintervals [r;_1,7;] of length Ar = (b — a)/m and 
divide the interval [a, 6] into n subintervals [0;-1, 6;| of width A@ = (8 — a)/n. This means that the 
circles r = r; and rays 0 = 0; for 1 <i < mand1 < j < n divide the polar rectangle R into smaller 
polar subrectangles R;; ([link](b)). 


(rj, 0%) 


ty 
ple 
hin 
A@ Wy M7, os 
Wy ae 


Wao 
‘mee? 


(a) (b) (c) 


(a) A polar rectangle R (b) divided into subrectangles R;;. (c) Close-up of a subrectangle. 


As before, we need to find the area AA of the polar subrectangle R;; and the “polar” volume of the thin 
box above R;;. Recall that, in a circle of radius r, the length s of an arc subtended by a central angle of 0 
radians is s = r@. Notice that the polar rectangle R;; looks a lot like a trapezoid with parallel sides r;_, A 
and r;Aé@ and with a width Ar. Hence the area of the polar subrectangle R;; is 

Equation: 


AA= 5A (r;-1A0 +71A8). 


Simplifying and letting ry = $(ri-1 +7r;),wehave AA = r,,ArAd. Therefore, the polar volume of the 


thin box above R;; ([link]) is 
Equation: 


* * * * 


Finding the volume of the 
thin box above polar 
rectangle R,;;. 


Using the same idea for all the subrectangles and summing the volumes of the rectangular boxes, we obtain 
a double Riemann sum as 
Equation: 


3 


i * * * 
SoS Flrijs Oj)rigAr 8. 


As we have seen before, we obtain a better approximation to the polar volume of the solid above the region 
R when we let m and n become larger. Hence, we define the polar volume as the limit of the double 
Riemann sum, 

Equation: 


* 


m : : 
Li mam bs Frys 65) Piz ATA. 


This becomes the expression for the double integral. 


Note: 

Definition 

The double integral of the function f (r, #) over the polar rectangular region R in the r6-plane is defined 
as 

Equation: 


Again, just as in Double Integrals over Rectangular Regions, the double integral over a polar rectangular 
region can be expressed as an iterated integral in polar coordinates. Hence, 


Equation: 
[[ 200.944 = ff t(0,8)rar a8 - /  t.0yedran 
R R 0=a r=a 


Notice that the expression for dA is replaced by r dr dO when working in polar coordinates. Another way 
to look at the polar double integral is to change the double integral in rectangular coordinates by 
substitution. When the function f is given in terms of x and y, using 

x =rcosé,y=rsin9, anddA = r dr dé changes it to 


Equation: 
i f(z,y)dA = // f (rcos 6,r sin 6)r dr dé. 
R R 


Note that all the properties listed in Double Integrals over Rectangular Regions for the double integral in 
rectangular coordinates hold true for the double integral in polar coordinates as well, so we can use them 
without hesitation. 


Example: 
Exercise: 


Problem: 
Sketching a Polar Rectangular Region 


Sketch the polar rectangular region R = {(r,@)|1 <r < 3,0 <6< 7a}. 


Solution: 


As we can see from [link], 7 = 1 andr = 3 are circles of radius 1 and 3 and 0 < 0 < 7 covers the 
entire top half of the plane. Hence the region R looks like a semicircular band. 


The polar region R lies 
between two semicircles. 


Now that we have sketched a polar rectangular region, let us demonstrate how to evaluate a double integral 
over this region by using polar coordinates. 


Example: 
Exercise: 


Problem: 
Evaluating a Double Integral over a Polar Rectangular Region 


Evaluate the integral // 3a dA over the region R = {(r, 6)|l <r<2,0<0< 7}. 
R 


Solution: 


First we sketch a figure similar to [link] but with outer radius 2. From the figure we can see that we 
have 


Equation: 
Te Use an iterated integral with correct limits 
32dA = J [: 3r cos Or dr dO ; ; s 
of integration. 
R 6=0 r=1 
= cos 6 [r° ie dé Integrate first with respect tor. 
6=0 
0=1 


= 7 cos 6d0 = 7 sin 6|5_. = 0. 


Note: 
Exercise: 


Problem: Sketch the region R = iG. ONL ie ae SS a}, and evaluate |[« dA. 


Solution: 


14 
3 
Hint 


Follow the steps in [link]. 


Example: 
Exercise: 


Problem: 
Evaluating a Double Integral by Converting from Rectangular Coordinates 


Evaluate the integral | / (l 77 = y’)dA where R is the unit circle on the xy-plane. 


Solution: 


The region R is a unit circle, so we can describe itas R = {(r,0)|0 <r < 1,0 <@< 2r}. 


Using the conversion x = rcos 0, y = rsin 6, and dA = r dr dO, we have 


Equation: 
2x 1 
ffa-2-v)aa =f fo- jarae= f flr )dr dO 
R 0 0 
20 
2 471 
=| leva pa f done 
2 4} 4 2 
0 0 
Example: 
Exercise: 
Problem: 


Evaluating a Double Integral by Converting from Rectangular Coordinates 


Evaluate the integral // (a + y)dA where R = {(2z,y)|1 <a? +y? <4,2 <0}. 


Solution: 


We can see that A is an annular region that can be converted to polar coordinates and described as 
Re) 2 Sn} (see the following graph). 


Ya 


The annular region of 
integration R. 


Hence, using the conversion x = rcos 0, y = rsin 0, and dA = r dr dO, we have 
Equation: 


0=37r/2 r=2 
[fternaa = / / (r cos 8 + r sin 6)r dr dO 
R O=n/2  r=1 
r=2 37/2 
= r’dr / (cos @ + sin 0)d0 
r=1 m/2 
eo iFepes 3/2 
= B [sin 6 — cos Al /2 
— 14 
=-%. 


Note: 
Exercise: 


Problem: 


Evaluate the integral i / (4 ie y’) dA where R is the circle of radius 2 on the zy-plane. 
R 


Solution: 


87 


Hint 


Follow the steps in the previous example. 


General Polar Regions of Integration 


To evaluate the double integral of a continuous function by iterated integrals over general polar regions, we 
consider two types of regions, analogous to Type I and Type II as discussed for rectangular coordinates in 
Double Integrals over General Regions. It is more common to write polar equations as r = f (0) than 

6 = f (r), so we describe a general polar region as R = {(r, 0)|a < 0 < B,hi (0) <r < hg (0)} (see the 
following figure). 


r=h,(0) 


A general polar region between 
a<6< Bandh; (6) <r < ho (8). 


Note: 
Double Integrals over General Polar Regions 
If f (r, 0) is continuous on a general polar region D as described above, then 


Equation: 
0=6 r=h2(0) 
| f (r, 0)r dr dd = ii / f (r, 6)r dr dO 
D 6=a r=h,(@) 
Example: 


Exercise: 


Problem: 
Evaluating a Double Integral over a General Polar Region 


Evaluate the integral i i r’sin Or dr dO where D is the region bounded by the polar axis and the 


D 
upper half of the cardioid r = 1 + cos 6. 


Solution: 


We can describe the region D as {(r,0)|0 < 6 < 2,0 < r < 1+ cos 9} as shown in the following 
figure. 


r=1+cos¢ 


05 10 15 2,0 


The region D is the top half of a cardioid. 


Hence, we have 
Equation: 


6=n r=1+cosé 


// r’sin 6rdrd@ = / i (r?sin 0)r dr dO 
D 


6=0 r=0 
9 r=1-+cos 6 


= + i Leal sin 6 dé 
=0 
g= 


(1+cos 6)? a 8 
0 


5 ah 


AIH 
| ea | 


Note: 
Exercise: 


Problem: Evaluate the integral 


Equation: 
[[Psx26r dr d0 where D = { (r,)|0 <0< 7,0 <r < 2Vcos 20}. 
D 
Solution: 
1/8 
Hint 


Graph the region and follow the steps in the previous example. 


Polar Areas and Volumes 


As in rectangular coordinates, if a solid S' is bounded by the surface z = f (r, 6), as well as by the surfaces 
r=a,r = b,0 =a, and 0 = £, we can find the volume V of S by double integration, as 


Equation: 
9=8 
V= [[ #00.9)rara9 = i: 
R 6=a T 


If the base of the solid can be described as D = {(r, 0)|a < 6 < B,h, (0) <r < hg (6)}, then the double 
integral for the volume becomes 


Equation: 
6=B 
V= [[ #00.9)rara9 = / 
D 0=a 


We illustrate this idea with some examples. 


r=b 
f (r, 9)r dr dé. 


=a 


r=h,(0) 
f (r, 9)r dr dé. 
r=h, (8) 


Example: 
Exercise: 


Problem: 
Finding a Volume Using a Double Integral 


Find the volume of the solid that lies under the paraboloid z = 1 — 2? — y? and above the unit circle 
on the zy-plane (see the following figure). 


1.0 


The paraboloidz =) — «> — 47, 


Solution: 


By the method of double integration, we can see that the volume is the iterated integral of the form 
if ( = y’)dA where R = {(r,0)|0 <r < 1,0 < 6 < 2z}. 


This integration was shown before in [link], so the volume is + cubic units. 


Example: 
Exercise: 
Problem: 
Finding a Volume Using Double Integration 


Find the volume of the solid that lies under the paraboloid z = 4 — x? — y? and above the disk 


(x - ly + y” = 1 on the zy-plane. See the paraboloid in [link] intersecting the cylinder 
(x — 1)? + y* = 1 above the ay-plane. 


(x-1)%+y*=1 


x 
Finding the volume of a solid with a 
paraboloid cap and a circular base. 
Solution: 


First change the disk (2 — Dy: + y” = 1 to polar coordinates. Expanding the square term, we have 
x? — 22 +1+y? = 1. Then simplify to get x? + y? = 22, which in polar coordinates becomes 

r? = 2r cos @ and then either r = 0 or r = 2 cos @. Similarly, the equation of the paraboloid changes 
to z = 4 — r*. Therefore we can describe the disk (x — 1)? + y? = 1 on the xy-plane as the region 
Equation: 


D= {(r,9)|0 <0< 7,0 <r < 2cos6}. 


Hence the volume of the solid bounded above by the paraboloid z = 4 — x? — y? and below by 
r= 2cos@is 
Equation: 


6=n r=2cosé 


|] #0.9)rarao = / / (4 —r?)r dr dO 
D 6=0 r=0 
=n 
D) 4 |2cos0 
i a = — |a 


0 
= / [8 cos?0 — 4 cos”6] dé = [30 ++ sin 6 cos 0 — sin 6 cos*6] ¢ = on. 
0 


=< 
I 


Notice in the next example that integration is not always easy with polar coordinates. Complexity of 
integration depends on the function and also on the region over which we need to perform the integration. 
If the region has a more natural expression in polar coordinates or if f has a simpler antiderivative in polar 
coordinates, then the change in polar coordinates is appropriate; otherwise, use rectangular coordinates. 


Example: 
Exercise: 


Problem: 
Finding a Volume Using a Double Integral 


Find the volume of the region that lies under the paraboloid z = x” + y? and above the triangle 
enclosed by the lines y = x, x = 0, and x + y = 2 in the zy-plane ([link)). 


Solution: 


First examine the region over which we need to set up the double integral and the accompanying 
paraboloid. 


Finding the volume of a solid under a paraboloid and above a given triangle. 


The region D is {(z, y)|(0 < # <1,4 < y < 2 —<z}. Converting the lines y = x, x = 0, and 

x + y = 2 in the xy-plane to functions of r and 0, we have 6 = 7/4, 6 = 7/2, and 

r = 2/ (cos @ + sin 8), respectively. Graphing the region on the zy-plane, we see that it looks like 
D= {(r,6)|7/4 < 6 < 2/2,0 <r < 2/ (cos6 + sin 6)}. Now converting the equation of the 
surface gives z = x? + y? = r?. Therefore, the volume of the solid is given by the double integral 
Equation: 


6=n/2 r=2/(cos 6+sin 8) na, BUCS b Cstne) 
V = ff f (r,0)r dr dO = / | rr dr dO = i B do 
0= a 1/4 0 
n/2 ; n/2 , = n/2 : = 
=if (sats) one (sors) a=4 | (<crrus) oe 
n/4 n/4 /4 


As you can see, this integral is very complicated. So, we can instead evaluate this double integral in 
rectangular coordinates as 
Equation: 


1 
v=) (a? + y*)dy dz. 
0 a 


Evaluating gives 


Equation: 
1 2-2 1 2-2 
ye 

vf fi a? + y? Jayae = f |arv+ 4] dx 
0 0 - 
a 3 

= [zt t 40 —— ar 
0 

_ | 8x Ag? a4 oe 

= |% -20 wees =||,=- 


To answer the question of how the formulas for the volumes of different standard solids such as a sphere, a 
cone, or a cylinder are found, we want to demonstrate an example and find the volume of an arbitrary cone. 


Example: 
Exercise: 


Problem: 
Finding a Volume Using a Double Integral 


Use polar coordinates to find the volume inside the cone z = 2 — ,/x? + y? and above the 
xy-plane. 


Solution: 


The region D for the integration is the base of the cone, which appears to be a circle on the xy-plane 
(see the following figure). 


Finding the volume of a solid inside the cone and above the ry- 
plane. 


We find the equation of the circle by setting z = 0: 
Equation: 

0 = 2-V22+y¥ 
oe 
aa = A. 


bo 
II 


This means the radius of the circle is 2, so for the integration we have 0 < 0 < 27 and0 <r < 2. 


Substituting z = r cos @ and y = rsin 0 in the equation z = 2 — \/z? + y? we have z =2 —r. 
Therefore, the volume of the cone is 


O=2n r=2 
(2 —r)rdr d= ans = — cubic units. 
6=0 r=0 
Analysis 


Note that if we were to find the volume of an arbitrary cone with radius a units and height A units, then the 
equation of the cone would be z = h — a Jr? + y?. 


h 
We can still use [link] and set up the integral as / i (1 = +) r dr dé. 
a 


. ‘4 1 2 
Evaluating the integral, we get > 7a*h. 


Note: 
Exercise: 


Problem: 


Use polar coordinates to find an iterated integral for finding the volume of the solid enclosed by the 
paraboloids z = x? + y* and z = 16 — x? — 7. 


Solution: 
Qn 2/2 
v= i / (16 — rl dr d@ = 647 cubic units 
0 0 
Hint 


Sketching the graphs can help. 


As with rectangular coordinates, we can also use polar coordinates to find areas of certain regions using a 
double integral. As before, we need to understand the region whose area we want to compute. Sketching a 
graph and identifying the region can be helpful to realize the limits of integration. Generally, the area 
formula in double integration will look like 

Equation: 


B hal) 
Area A = / Ir dr dé. 


hi(9) 


Example: 
Exercise: 


Problem: 
Finding an Area Using a Double Integral in Polar Coordinates 


Evaluate the area bounded by the curve r = cos 40. 


Solution: 


Sketching the graph of the function r = cos 40 reveals that it is a polar rose with eight petals (see the 
following figure). 


Finding the area of a polar rose with 
eight petals. 


Using symmetry, we can see that we need to find the area of one petal and then multiply it by 8. 
Notice that the values of @ for which the graph passes through the origin are the zeros of the function 
cos 40, and these are odd multiples of 7/8. Thus, one of the petals corresponds to the values of @ in 
the interval [—7/8, 7/8]. Therefore, the area bounded by the curve r = cos 46 is 

Equation: 


6=7/8 r=cos 40 
A =& / i Ir dr dO 


0=—7/8 r=0 
7/8 7/8 
1 9) cos 40 1 2 1 il 6 1/8 T T 
= / Ei lo Ja=s / 7 49 d9 = 8 [40 + 4, sin 48 cos 40|""* = 5/4] =F: 
—1/8 —1/8 


Example: 
Exercise: 


Problem: 


Finding Area Between Two Polar Curves 


Find the area enclosed by the circle r = 3 cos @ and the cardioid r = 1 + cos 0. 


Solution: 


First and foremost, sketch the graphs of the region ([link]). 


o~ » \ 


r=1 + cosd 


Finding the area enclosed by both a circle and a cardioid. 


We can from see the symmetry of the graph that we need to find the points of intersection. Setting the 
two equations equal to each other gives 
Equation: 


3cos?é=1+cos 0. 


One of the points of intersection is 9 = 7/3. The area above the polar axis consists of two parts, with 
one part defined by the cardioid from 9 = 0 to 6 = 7/3 and the other part defined by the circle from 
6 = 7/3 to 0 = 7/2. By symmetry, the total area is twice the area above the polar axis. Thus, we 
have 

Equation: 


0=17/3 r=1+cos 0 0=n/2 r=3cos6 


A =? / / Ir dr dé + i / Ir dr d@}|. 


6=0 r=0 d=1/3 r=0 


Evaluating each piece separately, we find that the area is 


Equation: 
1 9 ~- 3 9 - 5 5 
FAG —e2) 7+ V3+4+—n V3) =2 a ) = —7 square units. 
4 16 8 16 8 4 
Note: 
Exercise: 
Problem: 


Find the area enclosed inside the cardioid r = 3 — 3 sin 8 and outside the cardioid r = 1 + sin 0. 


Solution: 
7/6 3-3sin0 
A= 2 i / r dr dO = 84 + 9V3 
—n/2 1+sin0 
Hint 


Sketch the graph, and solve for the points of intersection. 


Example: 
Exercise: 


Problem: 
Evaluating an Improper Double Integral in Polar Coordinates 


Evaluate the integral | / e 10(2"+4") dp dy. 
R2 


Solution: 


This is an improper integral because we are integrating over an unbounded region R?. In polar 
coordinates, the entire plane R? can be seen as 0 < 0 < 27,0 <r<o. 


Using the changes of variables from rectangular coordinates to polar coordinates, we have 
Equation: 


[fee Mae dy = i il e 0» dr db = / 
=0 


R? 6 


=27r] lim [as 


Note: 
Exercise: 


Problem: Evaluate the integral | / e 4(@"+¥") dx dy. 
R2 


Solution: 


®|y 


Hint 


Convert to the polar coordinate system. 


Key Concepts 


e To apply a double integral to a situation with circular symmetry, it is often convenient to use a double 
integral in polar coordinates. We can apply these double integrals over a polar rectangular region or a 
general polar region, using an iterated integral similar to those used with rectangular double integrals. 


e The area dA in polar coordinates becomes r dr dé. 


e Use x = rcos6, y= rsin 9, and dA = r dr dé to convert an integral in rectangular coordinates to an 


integral in polar coordinates. 


° User? = 27+ y? and@ = tan ! (+) to convert an integral in polar coordinates to an integral in 


rectangular coordinates, if needed. 


¢ To find the volume in polar coordinates bounded above by a surface z = f (r, 8) over a region on the 


zy-plane, use a double integral in polar coordinates. 


Key Equations 


¢ Double integral over a polar rectangular region R 
| (i f(r,@)dA = lim SOS OF (ri*, 07°) AA = lim SONOS (ris*, 0ij* )rig* Ar M0 
eS ey gal MONEE RN Egy 
R 


¢ Double integral over a general polar region 


6=8 r=h2(6) 
[fs 6)r dr dO = i il f (r,4)r dr d6 
D 6=a_ r=h,(6) 


In the following exercises, express the region D in polar coordinates. 
Exercise: 


Problem: D is the region of the disk of radius 2 centered at the origin that lies in the first quadrant. 
Exercise: 


Problem: 


D is the region between the circles of radius 4 and radius 5 centered at the origin that lies in the 
second quadrant. 


Solution: 


D= {(r,)|4<r<5,5 <0<nx} 


Exercise: 


Problem: D is the region bounded by the y-axis and x = \/1 — y?. 


Exercise: 


Problem: D is the region bounded by the z-axis and y = 2 — 2?. 


Solution: 


D= {(r,6)|0 <r < v2,0<0<-n} 


Exercise: 


Problem: D = {(z, y)|x? + y? < 4x} 
Exercise: 

Problem: D = {(x, y)|x? + y* < 4y} 

Solution: 


D= {(r,0)|0 <r < 4sin0,0 < 0< 7} 


In the following exercises, the graph of the polar rectangular region D is given. Express D in polar 
coordinates. 
Exercise: 


Problem: 


—4 —2 0 2 4 


Exercise: 


Problem: 


Solution: 


D={(,0)[3 <r < 5,7 <0< 


d|A 


} 


Exercise: 


Problem: 


Exercise: 


Problem: 


Solution: 

D= {(r,0)|3 <r <5, 3% <6< & 
Exercise: 

Problem: 


In the following graph, the region D is situated below y = = and is bounded by # = 1, = 5, and 
y=—0. 


Exercise: 


Problem: In the following graph, the region D is bounded by y = x and y = 2?. 


y 
5 


Solution: 


D = {(r,9)|0 <r < tanOsec0,0<0< F} 


In the following exercises, evaluate the double integral i / f (x, y)dA over the polar rectangular region D. 
R 


Exercise: 


Problem 


Exercise: 


Problem 


: f(,y) =2? + y?, D = {(r,0)|3 <7 <5,0 <0 < 2m} 


: f(x,y) =a+y, D={(r,0)|3 <r <5,0<0< 2m} 


Solution: 


0 


Exercise: 


Problem 


Exercise: 


Problem 


> f(z,y) =2?+2ay,D= {(r,A)|1 <r <2,7 <0 < 27r} 


fag) =2t+44,D= {0 |bxr<2, + <0 < 2a} 


Solution: 


637 
16 


Exercise: 


Problem 


Exercise: 


Problem 


: f (z,y) = Vx? + y*, where D = {(r,0)|0 <r <1,5 <O<7n}. 


: f(x,y) = 24 + 2x?y? + y*, where D = {(r,0)|3 <r <4,% <O< 2}. 


Solution: 


33677 
18 


Exercise: 


Problem 


Exercise: 


Problem 


: f(x,y) = sin (arctan ), where D = {(r,0)|1 <r<2,% <0< F} 


: f(x,y) = arctan (%), where D = {(r, 4)|2 SPa3) 4 SOS t} 


Solution: 


351? 
576 


Exercise: 


Problem: 


en [1 + 2arctan (4) ]aa,D = {(r, 1 <r<2,% <0< F} 
D 


Exercise: 


Problem: a Cre + ¢4 + 2n7y? + y') arctan (4)dA, D= {(r,)|1<r<2,3 <a< F} 
D 


Solution: 


er (21 —e+ e*) 


In the following exercises, the integrals have been converted to polar coordinates. Verify that the identities 
are true and choose the easiest way to evaluate the integrals, in rectangular or polar coordinates. 
Exercise: 


2 ©&£ 
Problem: [ [ (c? + y?)dyde =a i r°dr do 
1 0 0 secé 


Exercise: 


w/4 tan @sec0 
dydz = [ if r cos 0 dr d6 
0 


3 2 
Problem: i: / Z 
Solution: 
$In (3 + 2v2) 
Exercise: 


7/4 tan @sec0 


LP wt 
Problem: / i =e dz = ; / dr d0 
0 a2 Y wy? 0 0 


Exercise: 
1 2 7/4 tan @sec 0 
Problem: i; / a dz = ‘, / r sin 6 dr d@ 
2 2 
O° as MY 0 0 
Solution: 


3(2- v2) 


In the following exercises, convert the integrals to polar coordinates and evaluate them. 
Exercise: 


3 V9-¥ 
Problem: / (x? + y”)dx dy 
0 0 


Exercise: 


fo) 
eens 
— 
, 
d| 
 gtiaie. 
= 
a 
m4 
YS 
i) 
Q 
8 
Q 
< 


Problem: 
fag 
Solution: 
x 2 
: / r°dr dO = Bee 
0 0 
Exercise: 


1 V1-2x? 
Problem: i / (x + y)dy dx 
0 0 


Exercise: 
4 V16—a2 
Problem: / / sin (2” + y”)dy dz 
0 _V16—x2 
Solution: 
m/2 4 
il Jr sin (r?) dr dO = msin28 
—n/2 0 
Exercise: 
Problem: 


Evaluate the integral | / r dA where D is the region bounded by the polar axis and the upper half of 


D 
the cardioid r = 1+ cos 0. 


Exercise: 
Problem: 


Find the area of the region D bounded by the polar axis and the upper half of the cardioid 
r=1+cos@. 


Solution: 


3m 
4 


Exercise: 


Problem: 


Evaluate the integral i / r dA, where D is the region bounded by the part of the four-leaved rose 


D 
r = sin 26 situated in the first quadrant (see the following figure). 


Exercise: 


Problem: 


Find the total area of the region enclosed by the four-leaved rose r = sin 26 (see the figure in the 
previous exercise). 


Solution: 


dR 


Exercise: 
Problem: 
Find the area of the region D, which is the region bounded by y = V4— 27,2 = V3, x = 2, and 
y=0. 
Exercise: 
Problem: 


Find the area of the region D, which is the region inside the disk x? + y? < 4 and to the right of the 
line x = 1. 


Solution: 


+ (40 — 3v3) 
Exercise: 


Problem: 


Determine the average value of the function f (a, y) = x” +? over the region D bounded by the 
polar curve r = cos 26, where — 7 < @ < (see the following graph). 


Exercise: 


Problem: 


Determine the average value of the function f (2, y) = \/x? + y? over the region D bounded by the 
polar curve r = 3 sin 26, where 0 < 6 < 4 (see the following graph). 


Solution: 


16 
3a 


Exercise: 


Problem: 


Find the volume of the solid situated in the first octant and bounded by the paraboloid 
z= 1-42? — 4y? and the planes xz = 0, y = 0, and z = 0. 


Exercise: 


Problem: 


Find the volume of the solid bounded by the paraboloid z = 2 — 9x? — 9y? and the plane z = 1. 


Solution: 


ea 


Exercise: 


Problem: 


a. Find the volume of the solid S; bounded by the cylinder x? + y? = 1 and the planes z = 0 and 
z= 1. 

b. Find the volume of the solid Sz outside the double cone z? = x? + y’, inside the cylinder 
x? + y? = 1, and above the plane z = 0. 

c. Find the volume of the solid inside the cone z? = x” + y? and below the plane z = 1 by 
subtracting the volumes of the solids $ and S9. 


Exercise: 


Problem: 


a. Find the volume of the solid S; inside the unit sphere x? + y? + z? = 1 and above the plane 


z=0. 
b. Find the volume of the solid $2 inside the double cone (z — 1) = x? + y? and above the plane 
z=0. 


c. Find the volume of the solid outside the double cone (z — 1)" = x? + y? and inside the sphere 
eo a? 2" = 1, 


Solution: 


For the following two exercises, consider a spherical ring, which is a sphere with a cylindrical hole cut so 
that the axis of the cylinder passes through the center of the sphere (see the following figure). 


Exercise: 


Problem: 


If the sphere has radius 4 and the cylinder has radius 2, find the volume of the spherical ring. 
Exercise: 

Problem: 

A cylindrical hole of diameter 6 cm is bored through a sphere of radius 5 cm such that the axis of the 

cylinder passes through the center of the sphere. Find the volume of the resulting spherical ring. 

Solution: 


208. cm? 


Exercise: 
Problem: 
Find the volume of the solid that lies under the double cone z? = 4a? + 4y?, inside the cylinder 
az? + y* = a, and above the plane z = 0. 


Exercise: 


Problem: 


Find the volume of the solid that lies under the paraboloid z = x? + y?, inside the cylinder 
a? + y” = a, and above the plane z = 0. 
Solution: 


3 
32 


Exercise: 
Problem: 
Find the volume of the solid that lies under the plane x + y+ z = 10 and above the disk 
2 D2 
a+y’ =4e. 
Exercise: 
Problem: 


Find the volume of the solid that lies under the plane 2” + y + 2z = 8 and above the unit disk 
ge ty? =1. 


Solution: 

Ar 
Exercise: 

Problem: 


A radial function f is a function whose value at each point depends only on the distance between that 
point and the origin of the system of coordinates; that is, f(a, y) = g(r), where r = \/x? + y?. 


Show that if f is a continuous radial function, then // f (z,y)dA = (02 — 01) [(G(R2) — G(R1)], 


D 
where Gi(r) = rg(r) and (z,y) € D= {(r,6)|Ri <r < Ro,0 < 6 < 27}, withO < Ri < Ry 
and 0 < @, < 0) < 2r. 


Exercise: 


Problem: 


Use the information from the preceding exercise to calculate the integral i / (x? + y’) °dA, where D 
D 

is the unit disk. 

Solution: 


qT 


4 
Exercise: 


Problem: 


Let f(z, y) = ae be a continuous radial function defined on the annular region 


D={(r,0)|Ri <7 < Ro,0 < 0 < 2m}, where r = \/a? + y?, 0 < Ri < Ry, and Fisa 
differentiable function. Show that | / f(z, y)dA = 27 |F(R,) — F(R,)). 


D 
Exercise: 
Problem: 
evaty 
Apply the preceding exercise to calculate the integral —— dx dy, where D is the annular 
JJ a+ y? 
region between the circles of radii 1 and 2 situated in the third quadrant. 
Solution: 
+me(e — 1) 
5 Tele 
Exercise: 
Problem: 


Let f be a continuous function that can be expressed in polar coordinates as a function of @ only; that 
is, f(z, y) = h(6), where (x, y) € D= {(r,0)|Ri <r < Ro, 1 < 6 < 45}, withO < Ri < Rp 


and 0 < 6, < 62 < 27. Show that [[#@vaa = + (RB — R?) [H (62) — H(01)], where H is an 


D 
antiderivative of h. 


Exercise: 


Problem: 


2 
Apply the preceding exercise to calculate the integral i / a dA, where 
x 
D 


D= {(r,O)|1<r<2,% <0< F}. 


Solution: 


us 
J3—* 
Exercise: 
Problem: 
Let f be a continuous function that can be expressed in polar coordinates as a function of @ only; that 


is, f(z, y) = g(r)h(9), where (x,y) € D= {(r,0)|Ri <r < Ro, 01 < 6 < 62} with 
0< R, < Ro andO < 0, < O_ < 27. Show that 


// f (a, y)dA = [G(R2) — G(R1)] [H (82) — H(61)|, where G and H are antiderivatives of g and 


D 
h, respectively. 


Exercise: 


Problem: Evaluate // arctan (=) V2? + y2dA, where D = {(r,0)|2 22,7 205 ae 
x 
D 
Solution: 


1332 
864 


Exercise: 


Problem: A spherical cap is the region of a sphere that lies above or below a given plane. 


a. Show that the volume of the spherical cap in the figure below is quh (3a? + h?) : 


b. A spherical segment is the solid defined by intersecting a sphere with two parallel planes. If the 
distance between the planes is h, show that the volume of the spherical segment in the figure 
below is emh (3a? + 3b? + i), 


Exercise: 


Problem: 


In statistics, the joint density for two independent, normally distributed events with a mean 4 = O and 
2 y2 


at 


a standard distribution o is defined by p(z,y) = 5-e ? . Consider (X, Y), the Cartesian 
coordinates of a ball in the resting position after it was released from a position on the z-axis toward 
the xy-plane. Assume that the coordinates of the ball are independently normally distributed with a 
mean jt = 0 and a standard deviation of o (in feet). The probability that the ball will stop no more 


than a feet from the origin is given by P |X? + Y? < a?] = | [vc y)dy dz, where D is the disk 


D 
of radius a centered at the origin. Show that P [x Fae a’ a ee ae 


Exercise: 
Problem: 
The double improper integral i / e (2 +9"/2) dy dx may be defined as the limit value of the 


ney: ; 2 ea : 
double integrals | / e(-*+¥'/2) dA over disks D, of radii a centered at the origin, as a increases 


oOo wo 


Da 
without bound; that is, / i e (2? +9"/2) dy dx = lim fe dA. 
—00 —0O D: 


a. Use polar coordinates to show that i / e+?) dy da = Qn. 


=00 =0coO 


oe) 


b. Show that i edgy = /2n, by using the relation 


- [eoray ae | ree fees 


Glossary 


polar rectangle 
the region enclosed between the circles r = a and r = b and the angles 6 = a and 0 = {; it is 
described as R = {(r,0)|a<r<ba<0< B} 


Triple Integrals 


e Recognize when a function of three variables is integrable over a rectangular box. 

e Evaluate a triple integral by expressing it as an iterated integral. 

e Recognize when a function of three variables is integrable over a closed and bounded region. 
e Simplify a calculation by changing the order of integration of a triple integral. 

¢ Calculate the average value of a function of three variables. 


In Double Integrals over Rectangular Regions, we discussed the double integral of a function f(z, y) of two 
variables over a rectangular region in the plane. In this section we define the triple integral of a function f(z, y, z) 
of three variables over a rectangular solid box in space, IR°. Later in this section we extend the definition to more 
general regions in R°. 


Integrable Functions of Three Variables 


We can define a rectangular box Bin R? as B= {(z,y,z)|a<a<bc<y<de<z< f}. We followa 
similar procedure to what we did in Double Integrals over Rectangular Regions. We divide the interval |a, b] into J 
subintervals [x;_1, 2;] of equal length Az = “—"+, divide the interval [c, d] into m subintervals [y;_1, y;] of 
equal length Ay = 4% 
Az = +—*" , Then the rectangular box B is subdivided into Imn subboxes 


Bix = |@i-1, 24] * [yi—1, ys] X [2;-1, 2%], a8 shown in [link]. 


, and divide the interval |e, f] into n subintervals [z;_1, z;] of equal length 


A rectangular box in R* divided into subboxes by planes parallel to 
the coordinate planes. 


* * * 
For each i, j, and k, consider a sample point (x ik Yijko Zn) in each sub-box B;;,. We see that its volume is 
AV = AzAyAz. Form the triple Riemann sum 
Equation: 


3 


I m * * * 
S> » F(rizn) Vijko 2 jp) AtAydz. 


We define the triple integral in terms of the limit of a triple Riemann sum, as we did for the double integral in 
terms of a double Riemann sum. 


Note: 

Definition 

The triple integral of a function f(z, y, z) over a rectangular box B is defined as 
Equation: 


Ff (Liskes Yijns Zig ATAYAZ = |i f (x,y, z)dV 
z B 


if this limit exists. 


When the triple integral exists on B, the function f(z, y, z) is said to be integrable on B. Also, the triple integral 
exists if f(a, y, z) is continuous on B. Therefore, we will use continuous functions for our examples. However, 
continuity is sufficient but not necessary; in other words, f is bounded on B and continuous except possibly on the 
boundary of B. The sample point (z;,z, Vises Zign) can be any point in the rectangular sub-box B;,, and all the 


properties of a double integral apply to a triple integral. Just as the double integral has many practical applications, 
the triple integral also has many applications, which we discuss in later sections. 


Now that we have developed the concept of the triple integral, we need to know how to compute it. Just as in the 
case of the double integral, we can have an iterated triple integral, and consequently, a version of Fubini’s thereom 
for triple integrals exists. 


Note: 
Fubini’s Theorem for Triple Integrals 
If f(x, y, z) is continuous on a rectangular box B = [a,b] x [c,d] x |e, f], then 


Equation: 
g ad ob 
I fenaav= fff sen2tedyde 


This integral is also equal to any of the other five possible orderings for the iterated triple integral. 


For a, b,c, d,e, and f real numbers, the iterated triple integral can be expressed in six different orderings: 
Equation: 


f db 


/ / / Le aie) dnay a i ( / ( i f(x,y, z)da)dy)dz = / ( / ( / f(x,y, z)dx)dz)dy 


£ 


/ (ft / f(,y,2)dy)dz)dx = / i / f(x,y, 2)dy)dx)dz 
i ( / ( / f(z, y, z)dz)dx)dy = / ( / ( / f(a, y, z)dz)dy)dz. 


For a rectangular box, the order of integration does not make any significant difference in the level of difficulty in 
computation. We compute triple integrals using Fubini’s Theorem rather than using the Riemann sum definition. 
We follow the order of integration in the same way as we did for double integrals (that is, from inside to outside). 


Example: 


Exercise: 


Problem: 
Evaluating a Triple Integral 


al y=4 a=5 
Evaluate the triple integral it Hi / (a + yz”)da dy dz. 
z=) vy =—1 


Solution: 


The order of integration is specified in the problem, so integrate with respect to first, then y, and then z. 
Equation: 


= y=4 a=5 
a / i (x + yz”)dx dy dz 
z=0 y=2 L=— 
z= y=4 2 z=5 
= i) / — + xyz? dy dz Integrate with respect to z. 
z=0 y=2 2 ¢=— 
= y=4 
= i i [12 + 6yz”| dy dz Evaluate. 
z=0 y=2 
— y? y=4 
= / 12y + cas ze dz Integrate with respect to y. 
=) y=2 
= 
= / [24 sF 362° | dz Evaluate. 
z=0 
31 2=1 
= [242 of 365 | Nae 36. Integrate with respect to z. 
ae 


Example: 
Exercise: 


Problem: 
Evaluating a Triple Integral 


Evaluate the triple integral [fev dV where B = {(z,y,z)| -2<2<1,0<y<3,1<2z<5has 
B 


shown in the following figure. 


Evaluating a triple integral over a given 
rectangular box. 


Solution: 


The order is not specified, but we can use the iterated integral in any order without changing the level of 
difficulty. Choose, say, to integrate y first, then x, and then z. 


i 8 1 
[fea =f | [iorvidaess= [Jf | Jax dz 
B =2 


0 =2 


5 
i ic zZdndz— ae SA), = 
1 =2 1 


Now try to integrate in a different order just to see that we get the same answer. Choose to integrate with 
respect to z first, then z, and then y. 


mz dz— "#|® — 162. 


i 


Equation: 
S 8 i 3 § ‘i 
we 
[ff ewav =|) [tevleaa= f [Sv |acay 
B 0 1 -2 0 1 a 
3 5 3 3 
2 |? 
=f sededy= [ ws |= [ say — 90% = 8(9 — 0) = 162. 
0 1 0 z 0 
Note: 


Exercise: 


Problem: 
Evaluate the triple integral | I ( zsin x cos ydV where 


B 
B= {(z,y,z)|0 <a <7, # SUS Ug SENS 


Solution: 


| ]f -n2 cosyav = 8 


B 
Hint 


Follow the steps in the previous example. 


Triple Integrals over a General Bounded Region 


We now expand the definition of the triple integral to compute a triple integral over a more general bounded region 
E in R®. The general bounded regions we will consider are of three types. First, let D be the bounded region that 
is a projection of E onto the xy-plane. Suppose the region F in R? has the form 

Equation: 


E= {(2,y, 2)|(t,y) € D,ui(z,y) < 2 < u2(z, y)}- 


For two functions z = u;(a, y) and z = ua(z, y), such that ui (x, y) < ue(z, y) for all (a, y) in.D as shown in 
the following figure. 


Z = Ux, y) 


Z= Uy (x, y) 


We can describe region & as the space between 
ui(a, y) and u(x, y) above the projection D of 
F onto the xy-plane. 


Note: 

Triple Integral over a General Region 

The triple integral of a continuous function f(a, y, z) over a general three-dimensional region 
Equation: 


B={(#, y,z)|\(e,y) © D,ui(e, y) < z < vale, y)} 


in R*, where D is the projection of E onto the xy-plane, is 
Equation: 


Uo(zx,y) 
y) 


If Heri 7)4V — // / f(a, y, 2)dz| dA. 


(x 


Similarly, we can consider a general bounded region D in the xy-plane and two functions y = ui(z, z) and 
y = U2(z, 2) such that u(x, z) < w2(z, z) for all (x, z) in D. Then we can describe the solid region F in R? as 
Equation: 


E= {(z,y, z)|(x, z) € D,ui(z, z) sys ua(z, z)} 


where D is the projection of & onto the xy-plane and the triple integral is 


Equation: 
u2(x,z) 
|]fsevaa= ff] [ teweau)aa. 
E D uy (x,z) 


Finally, if D is a general bounded region in the yz-plane and we have two functions x = ui(y, z) and 
x = u2(y, z) such that ui(y, z) < we(y, z) for all (y, z) in D, then the solid region E in R? can be described as 
Equation: 


E= {(2,y,2z)|(y,z) € Duly, 2) < @ < ualy, z)} 


where D is the projection of & onto the yz-plane and the triple integral is 


Equation: 
u2(y,2) 
|]f t@vaav = |f : f(x,y, z)dx| dA. 
E D ur(y,z) 


Note that the region D in any of the planes may be of Type I or Type II as described in Double Integrals over 
General Regions. If D in the xy-plane is of Type I ({link]), then 
Equation: 


E= {(z,y, z)|a See b, gi (x) sys g2(x), ui(z, y) Ses u2(z, y)}. 


Z = Ux, Y) 


ween 


z= U(x, y) 


y = gi(x) 
: y = g(x) 


A box E where the projection D in the xy-plane is of Type I. 


Then the triple integral becomes 
Equation: 


g2(z) u2(x,y) 


[ff tensav= | / / f(a, y, z)dz dy dz. 


1(z) un (x,y) 


If D in the zy-plane is of Type II ([link]), then 
Equation: 


B= {(2,y,z)le <a <d,hi(a) <y < ha(e),wi(o,y) <2 < wale,y)}. 


Z = U2lx, y) 


A box E where the projection D in the xy-plane 
is of Type II. 


Then the triple integral becomes 


Equation: 
y=d x=ho(y) Z=Uy(z,y) 
/I/ f(x,y, z)dV =i / f(z, y, z)dz dz dy. 
_ y=c Ja=hi(y) J z=ui(z,y) 
Example: 
Exercise: 
Problem: 


Evaluating a Triple Integral over a General Bounded Region 


Evaluate the triple integral of the function f(z, y, z) = 5a — 3y over the solid tetrahedron bounded by the 
Dies 2 = 0,97 = 0,2 = 0, ancl a -- gb 2 = I. 


Solution: 


[link] shows the solid tetrahedron F and its projection D on the xy-plane. 


The solid E& has a projection D on the xy-plane of Type I. 


We can describe the solid region tetrahedron as 
Equation: 


E={(a,y,2|0<2<1,0<y<1-—2,0<z<1-2-y}. 


Hence, the triple integral is 


Equation: 
a=1 y=1—a z=1—2—-y 
/I/ Jif, y,2)aV — / / / (5x — 3y)dz dy dz. 
S z=0 y=0 z=0 


Zl ot) 

To simplify the calculation, first evaluate the integral / (5x2 — 3y)dz. We have 
2=0 

Equation: 


iP Ge — 3y)dz = (5a — 3y) (1—a2—y). 


(al 
Now evaluate the integral / (5a — 3y) (1 — x — y)dy, obtaining 
y=0 
Equation: 


[e390 — 2 - shdy = He — 1)*(62 - 1). 


Finally, evaluate 
Equation: 


tL 
/ 1(a— 1)?(6a —1)dz =5. 


Putting it all together, we have 


Equation: 
a=) y=1—2 z=1—-2-y 
/I/ f(x,y, z)dV = / / i, (5a — 3y)dzdydz =. 
3 z=0 Jy=0 Z=0 


Just as we used the double integral | / 1dA to find the area of a general bounded region D, we can use | i I 1dV 


D E 
to find the volume of a general solid bounded region E. The next example illustrates the method. 


Example: 
Exercise: 


Problem: 
Finding a Volume by Evaluating a Triple Integral 


Find the volume of a right pyramid that has the square base in the zy-plane [—1,1] x [—1, 1] and vertex at 
the point (0, 0, 1) as shown in the following figure. 
z 


y 


Finding the volume of a pyramid with a square base. 


Solution: 


In this pyramid the value of z changes from 0 to 1, and at each height z, the cross section of the pyramid for 
any value of z is the square [—-1+ z,1—z] x [-1+z,1-— 2]. Hence, the volume of the pyramid is 


[fer 


E 
Equation: 


B= ({(2,y,z)|0<z2<1,-1l+z<y<1-2,-l+z<a¢<1-z2}. 


Thus, we have 
Equation: 


z=1 y=1—z a2=1—z z=1 y=1-z z=1 ‘ 4 
|ffpav-] / h Ide dy de = [ / (2- 22)dyde = [ (2 — 2z)"dz = —. 
S z=0 y=14+z z=1+z z=0 y=14+z z=0 3 


Hence, the volume of the pyramid is + cubic units. 


Note: 
Exercise: 


Problem: 
Consider the solid sphere E = {(a, y, z)|z? + y? + 2? = 9}. Write the triple integral /I/ f(z, y, z)dV 
E 


for an arbitrary function f as an iterated integral. Then evaluate this triple integral with f(z, y, z) = 1. 
Notice that this gives the volume of a sphere using a triple integral. 


Solution: 


z=3 py=V9—2? — pz=/9—2?-y? 
|[fav-sf i, i ldzdy dz = 36r. 
3 2=—3 J y=—- 9-2? Jz J/9 ry? 


Hint 


Follow the steps in the previous example. Use symmetry. 


Changing the Order of Integration 


As we have already seen in double integrals over general bounded regions, changing the order of the integration is 
done quite often to simplify the computation. With a triple integral over a rectangular box, the order of integration 
does not change the level of difficulty of the calculation. However, with a triple integral over a general bounded 
region, choosing an appropriate order of integration can simplify the computation quite a bit. Sometimes making 
the change to polar coordinates can also be very helpful. We demonstrate two examples here. 


Example: 
Exercise: 


Problem: 


Changing the Order of Integration 


Consider the iterated integral 
Equation: 


g=1 y=z? 


= i 
/ / / f (a, y, z)dz dy dz. 
z=0 y=0 2=0 


The order of integration here is first with respect to z, then y, and then x. Express this integral by changing 
the order of integration to be first with respect to x, then z, and then y. Verify that the value of the integral is 
the same if we let f (x, y, z) = xyz. 


Solution: 


The best way to do this is to sketch the region E and its projections onto each of the three coordinate planes. 
Thus, let 


Equation: 
jo {(x,y, 2) |0 aca l0ayor Ua749)) 
and 
Equation: 
g=1 y=2 Z=y 
/ a f f(a, wadedvas= ff 1 L,Y, 2 
z=0 y=0 2=0 


We need to express this triple integral as 
Equation: 


y=d 2z=v2x(y) z=u2(y,z) 


i f (a, y, z)dx dz dy. 


y=C z=vui(y) z=ui(y,z) 


Knowing the region & we can draw the following projections ([link]): 


on the zy-plane is Di = {(z,y)|0< a <1,0<y<27}={(2,y)|0<y<1,Vy<2x <1}, 


on the yz-plane is Dy = {(y, A <r 1 yt}, and 


on the xz-plane is D3 = {(x, z)|0 <a <1,0<z<2’}. 


ys (1, 1) ‘ (1, 1) 74 (1, 1) 


(1,0) (0,0) 


The three cross sections of on the three coordinate planes. 


Now we can describe the same region E as tiga. z)|0 Syn 1 a2 50,79 2 05 1}, and 
consequently, the triple integral becomes 


Equation: 
y=d  2=v2(y) z=u2(y,z) y=l z=2? = x=1 
i / / f (a, y, z)dx dz dy = i / f (a, y, z)dx dz dy. 
y=C z=vi1(y) z=ur(y,z) y=0 2-0 r=V/y 


Now assume that f (x, y, z) = xyz in each of the integrals. Then we have 
Equation: 


e=1 y=x? oe z=1 y=a? z=1 
= os | J dy dz = ps dydz = ol 
2 2 12 
=) a=W os x=0 


/ yz dx dzdy 
y=0 2z=0 cA) 0 
y=1 zy" pal 
@ 
= i / [es |e dy 
25|\e= 
y=0 z=0 vy 
yee — ye y=1 2 y=1 
2 2 2,2 |2=Y 5 6 
-f f G-Pleae f PP EE |e fF -2)aa 
y=0 2=0 y=0 z=0 p20 


The answers match. 


Note: 
Exercise: 


Problem: Write five different iterated integrals equal to the given integral 


Equation: 
i I fs (x, y, z)dax dy dz. 
z=0 y=0 2=0 
Solution: 
z=4 g=V4—-z y=4-z y=4 z=4-y t=,/y 
(i) i / ik f (a, y, z)dy dx dz, (ii) I | fr f (a, y, z)dx dz dy, (iii) 
z=0 z=0 =r? 2=0 c= 
y=4 t= fy 2=4- g=2  y=4 2z=4-y 


v¢ f (x, y, z)dz dx dy, (iv) i fs f (2, y, z)dz dy dz, (v) 


Hint 


Follow the steps in the previous example, using the region F as 
{ (a, Uns 2z— 405 ys 4—-7205e5 Vy}; and describe and sketch the projections onto each of the 


three planes, five different times. 


Example: 
Exercise: 


Problem: 
Changing Integration Order and Coordinate Systems 


Evaluate the triple integral | \) / J x? + z?dV, where FE is the region bounded by the paraboloid 


E 
y = x? + 2? ((link]) and the plane y = 4. 


y 


Integrating a triple integral over a paraboloid. 


Solution: 


The projection of the solid region F onto the zy-plane is the region bounded above by y = 4 and below by 
the parabola y = a? as shown. 


Cross section in the zy-plane of the paraboloid 
in [link]. 


Thus, we have 
Equation: 


E={(2,y,2)|-25 252,07 <y<4-Vy-a@ <z< Vy— 2}. 


The triple integral becomes 
Equation: 


eo 


[ore 


z=2  y=4 J yx? 
/ / / Va? + 2%dz dy dz. 
z 


@=-2 y=n? ,_ J y—22 


This expression is difficult to compute, so consider the projection of & onto the xz-plane. This is a circular 
disc x? + z? < 4. So we obtain 


Equation: 
a=2 y=4 z=y/y-2? @=2 z=vV4-2? y=4 
|i V2? + 2dV = i / V2? + 2dzdydz = ‘| iI i Vx? + z2dy dz 
B 2 yaa? 2 ./y=a? 2 yo V4-a? yan? +? 


Here the order of integration changes from being first with respect to z, then y, and then z to being first with 
respect to y, then to z, and then to z. It will soon be clear how this change can be beneficial for computation. 
We have 


Equation: 
w=2 0 z=//4—2? y=A4 r= z=V4—22 
i il i Va? + 22dydzdz = i ‘| (4— = aN V 22 + 22dz de. 
G=—2 pei Waae ers C22 ga a 


Now use the polar substitution z = r cos 0, z = r sin 6, and dz dx = r dr dé in the xz-plane. This is 
essentially the same thing as when we used polar coordinates in the zy-plane, except we are replacing y by z. 
Consequently the limits of integration change and we have, by using r? = x? + z?, 


Equation: 
a=2 z=V 4-27 6=2n7 r=2 
it (4-2? — z*) Vaer2+z22dzedx = / / (4—r?)rr dr do 
C2 ey eer? 6=0  r=0 


: f 64 
|a-/ a6 = ee. 


0 


Average Value of a Function of Three Variables 


Recall that we found the average value of a function of two variables by evaluating the double integral over a 
region on the plane and then dividing by the area of the region. Similarly, we can find the average value of a 
function in three variables by evaluating the triple integral over a solid region and then dividing by the volume of 
the solid. 


Note: 


Average Value of a Function of Three Variables 
If f (x, y, z) is integrable over a solid bounded region EF with positive volume V (£), then the average value of 


the function is 


Equation: 


fave = an ff sem z)dV. 


Note that the volume is V (EZ) = /I/ 1dV. 
E 


Example: 
Exercise: 


Problem: 
Finding an Average Temperature 


The temperature at a point (2, y, z) of a solid E bounded by the coordinate planes and the plane 
ety+z=1isT (a, y, z) = (cy + 8z+ 20) °C. Find the average temperature over the solid. 


Solution: 


Use the theorem given above and the triple integral to find the numerator and the denominator. Then do the 
division. Notice that the plane x + y + z = 1 has intercepts (1,0, 0), (0, 1,0), and (0,0, 1). The region E 
looks like 
Equation: 


E=({(z,y,2z)|(0<2<1,0<y<1-2,0<z<1-—2-y}. 


Hence the triple integral of the temperature is 
Equation: 


=i Weilwe 74 


[[fteuaw- f 
E £=0) -ya0 z 
Ps 1 
The volume evaluation is V (EZ) = /I/ 1dV = / i ‘A ldzdy dx = = 
=o 


E z=0 y=—0 z 


ty 


1 

147 
/ (ay + 8z+ 20)dzdydzx = ag 
=0 


147/40 6(147) 441 , 
1/6 70 59 degrees Celsius. 


Hence the average value is Taye = 


Note: 
Exercise: 


Problem: 


Find the average value of the function f (x, y, z) = xyz over the cube with sides of length 4 units in the first 
octant with one vertex at the origin and edges parallel to the coordinate axes. 


Solution: 
ik ave — 8 
Hint 


Follow the steps in the previous example. 


Key Concepts 


¢ To compute a triple integral we use Fubini’s theorem, which states that if f (x, y, z) is continuous on a 
rectangular box B = [a,b] x [c,d] x [e, f], then 


Equation: 
f db 
|[fseusa= fof [ teu2ardya: 
B e ca 


and is also equal to any of the other five possible orderings for the iterated triple integral. 
e To compute the volume of a general solid bounded region F we use the triple integral 


Equation: 
= fffoav. 
E 


Interchanging the order of the iterated integrals does not change the answer. As a matter of fact, interchanging 
the order of integration can help simplify the computation. 
e To compute the average value of a function over a general three-dimensional region, we use 


Equation: 
Fave = rw I lf penne 
E 


Key Equations 


¢ Triple mastiaies 


I n 
fs Dade aslo? gk) weetuds fff (2, y, 2 


In the following exercises, evaluate the triple integrals over the rectangular solid box B. 
Exercise: 


Problem: /I/ (2a + 3y + 42°) dV, where B= {(z,y,z)|O0<ax<10<y<2,0<z< 3} 
B 


Solution: 


192 


Exercise: 


Problem: /I/ (xy + yz + xz)dV, where B= {(z,y,z)|1 << 2 <2,0<y<2,1<2z< 3} 
B 


Exercise: 


Problem: /I/ (x cos y+ z)dV, where B= {(z,y,z)|0< ¢<1,0<y<a,-l<z<1} 


Solution: 


0 


Exercise: 


Problem: /I/ (zsina + y”)dV, where B= {(z,y,2)|0< a <2,0<y<1,-1<2z< 2} 


In the following exercises, change the order of integration by integrating first with respect to z, then x, then y. 
Exercise: 


1 2 3 
Problem: / i. J a? +Iny+ z) dz dy dz 
Or d2 


Solution: 


203% 4 
[of fet rimus ddedeay= 8 +2102 
LE, 2 


0 


Exercise: 
t 73 
Problem: , i fi (ze” + 2y)dx dy dz 
0 -1 0 
Exercise: 
2 3 4 
Problem: [ff z+ «de dy dz 
-1 1 0 
Solution: 
3 04 2 
1 
J/ [ (222+ =) aede dy = 044 12103 
1 0 -1 e 
Exercise: 
2 -1 1 
Problem: / / _ rag dy dz 
1 0 


Exercise: 


Problem: 


Let F’, G, and H be continuous functions on [a, 8], [c, d], and [e, f], respectively, where a, b,c, d,e, and f 
are real numbers such that a < b,c < d, ande < f. Show that 


Equation: 
b d f b d f 
/ / [Fos (y)H (z)dzdydz = [Fae [ewes [#@a 
Exercise: 
Problem: 


Let F', G, and H be differential functions on [a, 6], [c, d], and [e, f], respectively, where a, b,c, d,e, and f 
are real numbers such that a < b,c < d, ande < f. Show that 
Equation: 


b d f 
) / F' ()G! (y)H! (2)dz dy dx = [F (6) — F(a)] [(G(@) — G(o)] [H() — H(0))- 


In the following exercises, evaluate the triple integrals over the bounded region 
E={(a,y,z)|a<a<b,hi (2) <y<ho(2),e<2z< fh. 
Exercise: 


Problem: /I/ (2a + 5y + 7z)dV, where F = {(2,y,z)|0< a2 <10<y<-r#+11<z<2} 
E 


Solution: 


11 
12 


Exercise: 


Problem: /I/ (ylna + z)dV, where FE = {(z,y,z)|l <2 <e,0<y<Inz,0<z<1} 
E 
Exercise: 


Problem: /I/ (sin x + sin y)dV, where E = { (x,y, z)\O<a< 4, cosz Sy <cosz,-l<z< 1} 
E 


Solution: 


2 


Exercise: 


Problem: /I/ (cy + yz + xz)dV, where E = {(z,y,2z)|0< a <1,-2? <y<2?,0<z<1} 
E 


In the following exercises, evaluate the triple integrals over the indicated bounded region EF. 
Exercise: 


Problem: [fe 2ve)av, where £ = {(2,y,2)|0< @ < 1O<y<2,0<27<5-2-y} 
E 


Solution: 


439 
120 


Exercise: 


Problem: /I/ (a? + y° + 2°)dV, where E = {(z,y,2)|0 <a <2,0<y< 22,0<2<4-2-y} 
E 


Exercise: 


Problem: 


[ff vav, wnece 2 = {(e,y,2)|-1 <0 <1,-vI= a8 <y< VI=a80 <z<1-a?—y'} 
B 


Solution: 


0 
Exercise: 


Problem: 


[ff eav where B= {(2,y,2)| -2 <2 <2,-4vI— a < ys Yaw <2 <4-2t-y} 
E 


In the following exercises, evaluate the triple integrals over the bounded region F of the form 


E= {(2,y,2)|g (y) <2 <go(y),ce<y<de<z< fy. 
Exercise: 


Problem: |[feav, where EF = {(x,y,2)|1 -y<a2<y-1,-l<y<11<z< 2} 
E 
Solution: 


_ 64 
105 


Exercise: 


Problem: [ff ine av, where B = {(z,y,2)| -y4<a<yt0<y<2,0<2< 4} 
E 


Exercise: 


Problem: |]fe-woav, where & = { (2, y, z)| —y$<a<VWy¥0<y<le,-1l<z<1} 
B 


Solution: 


a 
26 


Exercise: 


Problem: [ff av net B {(x,y,2)|2-2y<e2<2+ /¥,0<y<1e,2<z< 3} 
E 


In the following exercises, evaluate the triple integrals over the bounded region 
Equation: 


E= {(a, y, 2)|g1 (y) <S2U2<cgo (y),¢ <y<dju (x,y) SzSug (x, y)}- 


Exercise: 


Problem: [ff av net B= { (x,y, 2)| —-y<a<y,0<y<1,0<2<1-24 —y*} 
B 
Solution: 


113 
450 


Exercise: 


Problem: |]ftee+ 1)dV, where E = { (x,y, z)|0 £@<4/9,0 < 752,028 1— 2? —y?} 
E 


Exercise: 


Problem: 


[ff e—2av, where 2 = {(0,y,2)|- Vimw <2 <y0<y< 42,0<2<1-2?-y} 
E 


Solution: 


ak (6v3 a 41) 


Exercise: 


Problem: [ff e+ wav, were £ = {(z,y,2))0<2< VJ1—y2,0<y<12,0<z< 1-2} 
B 


In the following exercises, evaluate the triple integrals over the bounded region 


E= {(a, y, z)| (x,y) € D, uy (x,y) < z < Ue (a, y)}, where D is the projection of E onto the xy-plane. 
Exercise: 


2 
Problem: // fe + z)dz | dA, where D = {(z,y)|a? + y? < 1} 
D \I 


Solution: 


3H 
2 


Exercise: 


3 
Problem: // [eer 1)dz } dA, where D = {(z,9)|2” -y >l1a< v5} 


D 1 
Exercise: 
10—x—y 
Problem: // / (a + 2z)dz | dA, where D = {(z, y)|y > 0,4 > 0,4 +y< 10} 
D 
Solution: 
1250 
Exercise: 
4x? +4y? 
Problem: // / ydz | dA, where D = {(z,y)|z? + y? <4,y> 1,2 >0} 
D 0 
Exercise: 
Problem: 


The solid E bounded by y? + z? = 9, z = 0, and x = 5 is shown in the following figure. Evaluate the 
integral | if / zdV by integrating first with respect to z, then y, and then z. 
E 


Solution: 


5 9-y? 


3 v/ 
/ / zdzdydzx = 90 
3 0 


0 
Exercise: 


Problem: 


The solid E bounded by y = »/z, « = 4, y = 0, and z = 1 is given in the following figure. Evaluate the 
integral | i} / xyz dV by integrating first with respect to x, then y, and then z. 


Exercise: 
Problem: 
0 0 a?+y? 
[T] The volume of a solid F is given by the integral / / i dz dy dx. Use a computer algebra system 
—2 £ 0 


(CAS) to graph & and find its volume. Round your answer to two decimal places. 


Solution: 


V = 5.33 


ty 
Mttrh 


u 
= hf 


Exercise: 


Problem: 


0 0 l4+y/a?+y? 
[T] The volume of a solid F is given by the integral / / i dz dy dz. Use a CAS to graph & and 
-1 _,? 0 


find its volume V. Round your answer to two decimal places. 


In the following exercises, use two circular permutations of the variables x, y, and z to write new integrals whose 
values equal the value of the original integral. A circular permutation of x, y, and z is the arrangement of the 
numbers in one of the following orders: y, z, and x or z,x, and y. 


Exercise: 
1 3 4 
Problem: ff [ (x*s* + 1)de dy dz 
0 1 2 
Solution 
13 4 1 3 4 
[ [ f@errvacacay fff (tv +iayacde 
0 i 2 0 i 2 
Exercise: 
3 1 -a+1 
Problem: | / / (2x + 5y + 7z)dy dz dz 
1 0 0 


Exercise: 


1 y 1-at—y! 


Problem: / / - In x dz dz dy 


0 -y 0 

Exercise: 
1 1 Wy 

Problem: / / / (x + yz)dx dy dz 
= oe eee 

Exercise: 

Problem: 


Set up the integral that gives the volume of the solid E bounded by y? = x? + z? and y = a”, where a > 0. 


Solution: 
a Vee @ 
V= / / ‘I dy dz dz 


Exercise: 


Problem: 


Set up the integral that gives the volume of the solid E bounded by x = y? + z? and x = a”, wherea > 0. 
Exercise: 
Problem: 


Find the average value of the function f (x, y, z) = x + y + z over the parallelepiped determined by 
z=0,c¢=1,y=0,y=3,z=0,andz=5. 


Solution: 


9 


2 
Exercise: 


Problem: 
Find the average value of the function f (a, y, z) = xyz over the solid FE = [0,1] x [0,1] x [0,1] situated 
in the first octant. 

Exercise: 


Problem: 


Find the volume of the solid & that lies under the plane x + y + z = 9 and whose projection onto the zy- 
plane is bounded by z = Jy —l,x=0,andz+y=7. 


Solution: 


156 
5 


Exercise: 


Problem: 


Find the volume of the solid E that lies under the plane 2x + y+ z = 8 and whose projection onto the ry- 
plane is bounded by « = sin~'y, y = 0, and x = 2. 
Exercise: 


Problem: 
Consider the pyramid with the base in the xy-plane of [—2, 2] x [—2, 2] and the vertex at the point (0, 0, 8). 
a. Show that the equations of the planes of the lateral faces of the pyramid are 4y + z = 8, 4y — z = —8, 
4c +2=8,and—4r%+2z=8. 
b. Find the volume of the pyramid. 


Solution: 


a. Answers may vary; b. — 
Exercise: 
Problem: 


Consider the pyramid with the base in the xy-plane of [—3,3] x [—3, 3] and the vertex at the point (0, 0, 9). 


a. Show that the equations of the planes of the side faces of the pyramid are 3y + z = 9, 3y+z=9, 
y =Oandz =0. 
b. Find the volume of the pyramid. 


Exercise: 


Problem: 


The solid E bounded by the sphere of equation x? + y? + z? = r? with r > 0 and located in the first octant 
is represented in the following figure. 


2.0 


1.5 


0.5 


0.0 


0.0 0.0 


0.5 


1.0 2 15 P 
y “2.920 


a. Write the triple integral that gives the volume of & by integrating first with respect to z, then with y, and 
then with x. 
b. Rewrite the integral in part a. as an equivalent integral in five other orders. 


Solution: 


4 Vri—a? ray? 2 VPP VPP Vie Ve 
) / / dz dy dz; b. / / / dz dz ay, | / / dy dz dz, 
0 0 0 

r — Vr—a—2 r — Jr-y-2 Tas r VrPa—y JrPa-y2 

/ | / dy dzdz, / | / dx dy dz, / / i. dx dz dy 
0 0 0 0 
Exercise: 

Problem: 


The solid E bounded by the equation 927 + 4y? + z? = 1 and located in the first octant is represented in the 
following figure. 


0.2 93 04 


0.5 


y 


a. Write the triple integral that gives the volume of F by integrating first with respect to z, then with y, and 
then with z. 


b. Rewrite the integral in part a. as an equivalent integral in five other orders. 


Exercise: 
Problem: 
Find the volume of the prism with vertices (0, 0,0), (2,0, 0), (2, 3,0), (0, 3,0), (0,0, 1), and (2,0, 1). 
Solution: 


3 
Exercise: 
Problem: 


Find the volume of the prism with vertices (0, 0,0), (4, 0, 0), (4,6, 0), (0, 6,0), (0,0, 1), and (4,0, 1). 
Exercise: 


Problem: 


The solid & bounded by z = 10 — 2a” — y and situated in the first octant is given in the following figure. 
Find the volume of the solid. 


10.0 


Solution: 


250 
3 


Exercise: 
Problem: 


The solid E bounded by z = 1 — z? and situated in the first octant is given in the following figure. Find the 
volume of the solid. 


Exercise: 


Problem: 
The midpoint rule for the triple integral | Ih / f (x, y, z)dV over the rectangular solid box Bisa 
B 


generalization of the midpoint rule for double integrals. The region B is divided into subboxes of equal sizes 
loom 


n 
and the integral is approximated by the triple Riemann sum x se S> f (&i, yj, ZR) AV, where (3, yj, 2%) 
i=1 j=l k=l 
is the center of the box B,;, and AV is the volume of each subbox. Apply the midpoint rule to approximate 


/I/ a” dV over the solid B = {(a, y, z)|0 <# <1,0<y<1,0 < z< 1} by using a partition of eight 


cubes of equal size. Round your answer to three decimal places. 
Solution: 
aw 
76 © 0.313 
Exercise: 


Problem: [T] 


a. Apply the midpoint rule to approximate | i} / e* dV over the solid 


B= {(a,y,z)|0< a2 <1,0<y<1,0 < z< 1} by using a partition of eight cubes of equal size. 
Round your answer to three decimal places. 

b. Use a CAS to improve the above integral approximation in the case of a partition of n° cubes of equal 
size, where n = 3,4,..., 10. 


Exercise: 
Problem: 


Suppose that the temperature in degrees Celsius at a point (2, y, z) of a solid & bounded by the coordinate 
planes andx +y+z=5isT (a,y, z) = 2z+5z-+ 10. Find the average temperature over the solid. 


Solution: 


35 
2 


Exercise: 
Problem: 
Suppose that the temperature in degrees Fahrenheit at a point (2, y, z) of a solid E bounded by the coordinate 
planes andx+y+z=5isT (z,y,z) = 2+ y+ cy. Find the average temperature over the solid. 
Exercise: 


Problem: 


Show that the volume of a right square pyramid of height h and side length a is v = hat by using triple 
integrals. 


Exercise: 


Problem: 


Show that the volume of a regular right hexagonal prism of edge length a is Sa°v3 by using triple integrals. 
Exercise: 


Problem: 


a? 


v3 by using triple integrals. 


Show that the volume of a regular right hexagonal pyramid of edge length a is —} 


Exercise: 


Problem: 
If the charge density at an arbitrary point (2, y, z) of a solid FE is given by the function p(z, y, z), then the 


total charge inside the solid is defined as the triple integral | i} / p(x, y, z)dV. Assume that the charge density 


of the solid F enclosed by the paraboloids x = 5 — y* — z? and x = y? + z” — 5 is equal to the distance 
from an arbitrary point of & to the origin. Set up the integral that gives the total charge inside the solid E. 


Glossary 
triple integral 


the triple integral of a continuous function f (x, y, z) over a rectangular solid box B is the limit of a Riemann 
sum for a function of three variables, if this limit exists 


Triple Integrals in Cylindrical and Spherical Coordinates 


e Evaluate a triple integral by changing to cylindrical coordinates. 
e Evaluate a triple integral by changing to spherical coordinates. 


Earlier in this chapter we showed how to convert a double integral in rectangular coordinates into a double 
integral in polar coordinates in order to deal more conveniently with problems involving circular symmetry. 
A similar situation occurs with triple integrals, but here we need to distinguish between cylindrical 
symmetry and spherical symmetry. In this section we convert triple integrals in rectangular coordinates into 
a triple integral in either cylindrical or spherical coordinates. 


Also recall the chapter opener, which showed the opera house |’ Hemispheric in Valencia, Spain. It has four 
sections with one of the sections being a theater in a five-story-high sphere (ball) under an oval roof as long 
as a football field. Inside is an IMAX screen that changes the sphere into a planetarium with a sky full of 
9000 twinkling stars. Using triple integrals in spherical coordinates, we can find the volumes of different 
geometric shapes like these. 


Review of Cylindrical Coordinates 


As we have seen earlier, in two-dimensional space R?, a point with rectangular coordinates (x, y) can be 
identified with (r, 0) in polar coordinates and vice versa, where z = rcos 0, y=rsin0,r? = 27 +y? 
and tan 0 = (=) are the relationships between the variables. 

In three-dimensional space R°, a point with rectangular coordinates (x, y, z) can be identified with 
cylindrical coordinates (r, 6, z) and vice versa. We can use these same conversion relationships, adding z as 
the vertical distance to the point from the zy-plane as shown in the following figure. 


P(x, y, Z) 
P(r, 0, 2) 


ee i | 


¥ 


Cylindrical coordinates are similar to polar 
coordinates with a vertical z coordinate 
added. 


To convert from rectangular to cylindrical coordinates, we use the conversion z = r cos @ and y = rsin@. 
To convert from cylindrical to rectangular coordinates, we use r? = 2? + y* and@ = tan! (+). The z- 
coordinate remains the same in both cases. 


In the two-dimensional plane with a rectangular coordinate system, when we say z = k (constant) we mean 
an unbounded vertical line parallel to the y-axis and when y = | (constant) we mean an unbounded 
horizontal line parallel to the z-axis. With the polar coordinate system, when we say r = c (constant), we 
mean a circle of radius c units and when # = a (constant) we mean an infinite ray making an angle a with 
the positive z-axis. 


Similarly, in three-dimensional space with rectangular coordinates (x, y, z), the equations z = k, y = l, 
and z = m, where k, 1, and m are constants, represent unbounded planes parallel to the yz-plane, xz-plane 
and xy-plane, respectively. With cylindrical coordinates (r, 6, z), by r = c, 9 = a, and z = m, where c, a, 
and m are constants, we mean an unbounded vertical cylinder with the z-axis as its radial axis; a plane 
making a constant angle a with the zy-plane; and an unbounded horizontal plane parallel to the zy-plane, 
respectively. This means that the circular cylinder x? + y? = c? in rectangular coordinates can be 
represented simply as r = c in cylindrical coordinates. (Refer to Cylindrical and Spherical Coordinates for 
more review.) 


Integration in Cylindrical Coordinates 


Triple integrals can often be more readily evaluated by using cylindrical coordinates instead of rectangular 
coordinates. Some common equations of surfaces in rectangular coordinates along with corresponding 
equations in cylindrical coordinates are listed in [link]. These equations will become handy as we proceed 
with solving problems using triple integrals. 


Circular 
cylinder Circular cone Sphere Paraboloid 


Rectangular et+y=c? fae (o* + y’) gat ae? 6" zZ=C io + y’) 


Cylindrical r=c Z=cr ries ¢ z=cr* 


Equations of Some Common Shapes 


As before, we start with the simplest bounded region B in R?, to describe in cylindrical coordinates, in the 
form of a cylindrical box, B = {(r,0,z)|a<r<b,a<6< B,c < z< d} ([link]). Suppose we divide 


each interval into 1, m and n subdivisions such that Ar = baa Ad = fa and Az = dae . Then we can 


state the following definition for a triple integral in cylindrical coordinates. 


Bie = (ia FX [0j-1, 9) X [2-15 24] 


A cylindrical box B described by 
cylindrical coordinates. 


Note: 

Definition 

Consider the cylindrical box (expressed in cylindrical coordinates) 
Equation: 


BO, aor a oO one 2 dy, 


If the function f (7,6, z) is continuous on B and if (risks Ges Zu) is any sample point in the cylindrical 
subbox By, = [ri-1,7i] < [0;-1,9;] < [Z-1, Ze] ([link]), then we can define the triple integral in 
cylindrical coordinates as the limit of a triple Riemann sum, provided the following limit exists: 
Equation: 


n 


lm * * * * 
lim » DOO Sie et pO eae 


j=1 j=l k=1 


Note that if g(a, y, z) is the function in rectangular coordinates and the box B is expressed in rectangular 


coordinates, then the triple integral i I / g(x, y, z)dV is equal to the triple integral 


B 
Il g(rcos 6,r sin 0, z)r dr d0 dz and we have 
B 


Equation: 
[[fsceuarav = fff a (7 cos0,rsind, 2)rar ap de = fff #(r,0,2)r dr ad az 
B B B 


As mentioned in the preceding section, all the properties of a double integral work well in triple integrals, 
whether in rectangular coordinates or cylindrical coordinates. They also hold for iterated integrals. To 
reiterate, in cylindrical coordinates, Fubini’s theorem takes the following form: 


Note: 

Fubini’s Theorem in Cylindrical Coordinates 

Suppose that g (x, y, z) is continuous on a rectangular box B, which when described in cylindrical 
coordinates looks like B = {(r,0,z)|a<r<ba<0<B,c<z<d}. 

Then g (z, y, z) = g(rcos6,rsin 6, z) = f (r, 90, z) and 


Equation: 
d B b 
|[fseue )dV = || | peeancae. 
B 


The iterated integral may be replaced equivalently by any one of the other five iterated integrals obtained 
by integrating with respect to the three variables in other orders. 


Cylindrical coordinate systems work well for solids that are symmetric around an axis, such as cylinders 
and cones. Let us look at some examples before we define the triple integral in cylindrical coordinates on 
general cylindrical regions. 


Example: 
Exercise: 


Problem: 
Evaluating a Triple Integral over a Cylindrical Box 


Evaluate the triple integral | i) / (zr sin 6)r dr d6 dz where the cylindrical box B is 
B 


B=A (re. 20S 20 8 7/2, 0s 4A 


Solution: 


As stated in Fubini’s theorem, we can write the triple integral as the iterated integral 


Equation: 
O=7/2 pr=2 pz=4 
iff (zr sin 0) \rdr dade = f i i. (zr sin 0)r dz dr dé. 
= r=0 Jz=0 


The evaluation of the iterated integral is straightforward. Each variable in the integral is independent 
of the others, so we can integrate each variable separately and multiply the results together. This 
makes the computation much easier: 

Equation: 


=n /2 z=4 
iL, ie | (zr sin 6)r dz dr dO 
d= r= z=0 
m/2 2 4 
= ul sade) (/ rar) ( zdz) = (- cos 60”) (5 
0 0 0 


bo 


4\ 64 
oJ 3° 


0) (3 


Note: 
Exercise: 


6= 
Problem: Evaluate the triple integral i) 
= 


i ii rzsin Or dz dr dé. 
=) 2=0) 


Solution: 


8 
Hint 


Follow the same steps as in the previous example. 


If the cylindrical region over which we have to integrate is a general solid, we look at the projections onto 
the coordinate planes. Hence the triple integral of a continuous function f (r, 6, z) over a general solid 
region E = {(r, 6, z)| (r,0) € D,ui(r, 0) < z < us (r,6)} in R3, where D is the projection of E onto 
the r6-plane, is 

Equation: 


u2(r,8) 


{ff r,6 Ones), [i Ga nae aeal 


ui(r,0) 


In particular, if D = {(r,6)|gi (0) <r < g2(0),a <6 < G}, then we have 
Equation: 


0=8 r=92(0) z=u2(r,0) 


[[[tosarara= [ / : f (r, 0, z)r dz dr dO. 


0=a r=gi\(0) z=u1 (7,0) 


Similar formulas exist for projections onto the other coordinate planes. We can use polar coordinates in 
those planes if necessary. 


Example: 
Exercise: 


Problem: 
Setting up a Triple Integral in Cylindrical Coordinates over a General Region 


Consider the region F& inside the right circular cylinder with equation r = 2 sin 0, bounded below by 
the r@-plane and bounded above by the sphere with radius 4 centered at the origin ([link]). Set up a 


triple integral over this region with a function f (7, 6, z) in cylindrical coordinates. 
z 


0,0, 4 
( ) r= 2sing 


r+z%=16 


Setting up a triple integral in cylindrical coordinates 
over a cylindrical region. 


Solution: 


First, identify that the equation for the sphere is r? + z? = 16. We can see that the limits for z are 


from 0 to z = V 16 — r?. Then the limits for r are from 0 to r = 2 sin @. Finally, the limits for 6 are 
from 0 to 7. Hence the region is 
Equation: 


B= {(r,8,z)|0 <6<7,0<7 < 2smn6,0<2< V'16 = 7}. 


Therefore, the triple integral is 
Equation: 


O=n r=2sin0 z=/16—r?2 
i / / f(r, 6, z)r dz dr dé. 
=0 0 


70 Z| 


II f (r, 0, z)r dz dr dO = 


0 


Note: 
Exercise: 


Problem: 
Consider the region F inside the right circular cylinder with equation r = 2 sin 0, bounded below by 


the r6-plane and bounded above by z = 4 — y. Set up a triple integral with a function f (7, , z) in 
cylindrical coordinates. 


Solution: 
O=n r=2sin@ z=4—rsin@ 
[[f 8.2 dedras = i / f (r, 6, z)r dz dr dé. 
E 6=0 7—0) = 
Hint 


Analyze the region, and draw a sketch. 


Example: 
Exercise: 


Problem: 
Setting up a Triple Integral in Two Ways 


Let & be the region bounded below by the cone z = vy x? + y? and above by the paraboloid 
z= 2-2? — y’. ([link]). Set up a triple integral in cylindrical coordinates to find the volume of the 
region, using the following orders of integration: 


a. dz dr dé 
b. dr dz dé. 


Setting up a triple integral in cylindrical coordinates over a 
conical region. 


Solution: 


a. The cone is of radius 1 where it meets the paraboloid. Since z = 2 — 2? — y* = 2— r” and 
z= \/x2 + y? =r (assuming r is nonnegative), we have 2 — r? = r. Solving, we have 
r? +7 —2= (r+ 2)(r—1) =0. Since r > 0, we have r = 1. Therefore z = 1. So the 
intersection of these two surfaces is a circle of radius 1 in the plane z = 1. The cone is the lower 
bound for z and the paraboloid is the upper bound. The projection of the region onto the ry- 
plane is the circle of radius 1 centered at the origin. 
Thus, we can describe the region as 
Equation: 


E={(r,6,z)|0<0<27,0<r<1r<z<2-r'}. 


Hence the integral for the volume is 
Equation: 


(=e geil ea 


= / i) i r dz dr dé. 


6=0 n=0) Z2=P 


b. We can also write the cone surface as r = z and the paraboloid as r? = 2 — z. The lower bound 
for r is zero, but the upper bound is sometimes the cone and the other times it is the paraboloid. 
The plane z = 1 divides the region into two regions. Then the region can be described as 
Equation: 


=A G2) |0 0 = 27,0 lO eh 
Uf (r,0,2)10 <8 < 2m,1<2<2,0<r< v2—2}. 


Now the integral for the volume becomes 


Equation: 
0—2 ae — 0=2n 2z=2 r=V2—-z 
v= / il [ rar dedi / / / r dr dz dé. 
6=0 =O 7=0 6=0 c= r=0 
Note: 
Exercise: 


Problem: Redo the previous example with the order of integration d6 dz dr. 


Solution: 


Hint 


Note that @ is independent of r and z. 


Example: 
Exercise: 


Problem: 
Finding a Volume with Triple Integrals in Two Ways 


Let E be the region bounded below by the r6-plane, above by the sphere x? + y* + z? = 4, and on 
the sides by the cylinder x? + y? = 1 ({link]). Set up a triple integral in cylindrical coordinates to 
find the volume of the region using the following orders of integration, and in each case find the 
volume and check that the answers are the same: 


a. dz dr d@ 
b. dr dz dé. 


x+y+227=4 


Finding a cylindrical volume with a triple integral in 
cylindrical coordinates. 


Solution: 


a. Note that the equation for the sphere is 
Equation: 


oe eee =4orr’?+2=4 
and the equation for the cylinder is 
Equation: 
x? +y? =lorr? =1. 


Thus, we have for the region # 
Equation: 


B= {(r,6,2)0 <2< V4—r?2,0<r< 1,0<6<2n} 


Hence the integral for the volume is 
Equation: 


V(E) = i r dz dr d6 
G—Omer—O) 2=0 
G=—2 7 o—l (=e  p=il 
— | bea | dr d§ = / / (rv 4— r?) dr d0 
6=0 r=0 6=0 r=0 
21 


b. Since the sphere is x? + y? + z” = 4, which is r? + z” = 4, and the cylinder is x? + y? = 1, 
which is r? = 1, we have 1 + z? = 4, that is, 2? = 3. Thus we have two regions, since the 


sphere and the cylinder intersect at (1, J 3) in the rz-plane 


Equation: 

Bi { (r,8,2)|0 <r<vV¥4—-73,V3<2<2,0<0< an} 
and 
Equation: 


i — {(r,8,2)0 <r < 1,0<2< V3,0<0< 2h. 


Hence the integral for the volume is 


Equation: 
0=2r 22 HV 4Sr" 6=27 z=V3 r=1 
VCE) i if ii r dr dzd0+ i i il r dr dz dO 
6-0 2-/3 7=0 (mec —() f=0) 
= V3n+ (2 _ 3V3)m = Oi (3 — v3) cubic units. 
Note: 
Exercise: 


Problem: Redo the previous example with the order of integration d6 dz dr. 


Solution: 


Ey = { (r,8,2)|0 SO eT Ona? lara = va—7?} and 
r=1 2z=V4—r? 0=20 
Vi= / / / rd dzdr. 
r=0-z=r ss 9=0 
Hint 


A figure can be helpful. Note that # is independent of r and z. 


Review of Spherical Coordinates 


In three-dimensional space R® in the spherical coordinate system, we specify a point P by its distance p 
from the origin, the polar angle 6 from the positive x-axis (same as in the cylindrical coordinate system), 
and the angle y from the positive z-axis and the line OP ([link]). Note that p > 0 and0 < » < a. (Refer 
to Cylindrical and Spherical Coordinates for a review.) Spherical coordinates are useful for triple integrals 
over regions that are symmetric with respect to the origin. 


P(x, y, 2) 
P(p, 9, ¢) 


The spherical coordinate system locates 
points with two angles and a distance from 
the origin. 


Recall the relationships that connect rectangular coordinates with spherical coordinates. 


From spherical coordinates to rectangular coordinates: 
Equation: 


x = psiny cos 0,y = psinysin 6, and z = pcos y. 


From rectangular coordinates to spherical coordinates: 
Equation: 


p=a?t+y?+27,tand = yy = arccos ey 
x VeP+yt+2 


Other relationships that are important to know for conversions are 


Equation: 
e r=psing 
§=8 These equations are used to convert from 
e = 
spherical coordinates to cylindrical coordinates 
°e z=pcosy 
and 
Equation: 
e p=Vvre4+22 
These equations are used to convert from 
e 60=0 cylindrical coordinates to spherical 


coordinates. 


e y=arccos (=) 
tf z 


The following figure shows a few solid regions that are convenient to express in spherical coordinates. 


o<0<> 5<c<e@ 


Sphere p = c (constant) Half plane @ = c (constant) Half cone ¢ = c (constant) 


Spherical coordinates are especially convenient for working with solids bounded by these types of 
surfaces. (The letter c indicates a constant.) 


Integration in Spherical Coordinates 


We now establish a triple integral in the spherical coordinate system, as we did before in the cylindrical 
coordinate system. Let the function f (~, 0, y) be continuous in a bounded spherical box, 
B= {(p,9,p)|a < p<ba<O0< B,y< y < qd}. We then divide each interval into 1,mandn 


subdivisions such that Ap = ba, Ab = fe Ay = oy 


Now we can illustrate the following theorem for triple integrals in spherical coordinates with 


(Dis 9 sus Pik) being any sample point in the spherical subbox B;;,. For the volume element of the 
subbox AV in spherical coordinates, we have AV = (Ap) (pAy) (p sin pAQ), , as shown in the 
following figure. 


p sin @ Ad 


AV = p* sin d Ap Ad AO 


The volume element of a box in spherical coordinates. 


Note: 


Definition 
The triple integral in spherical coordinates is the limit of a triple Riemann sum, 
Equation: 


I m 
4 * * * * » 
lim SS De i (Pijks Dies Lise) (Pisn) sin pApAdAy 


provided the limit exists. 


As with the other multiple integrals we have examined, all the properties work similarly for a triple integral 
in the spherical coordinate system, and so do the iterated integrals. Fubini’s theorem takes the following 
form. 


Note: 

Fubini’s Theorem for Spherical Coordinates 

If f (p, 9, p) is continuous on a spherical solid box B = [a,b] x [a, 6] x [y, a], then 
Equation: 


p=b 0=6 p=b 


/I/ f (0,9, p)p’sin v dp dy dd = ii / / f (0,9, p)p’sin v dp dy dé. 
B p= O=a p=a 


This iterated integral may be replaced by other iterated integrals by integrating with respect to the three 
variables in other orders. 


As stated before, spherical coordinate systems work well for solids that are symmetric around a point, such 
as spheres and cones. Let us look at some examples before we consider triple integrals in spherical 
coordinates on general spherical regions. 


Example: 
Exercise: 


Problem: 
Evaluating a Triple Integral in Spherical Coordinates 


0=2n y=r/2 p=l1 


Evaluate the iterated triple integral i / p’sin y dp dy dé. 


Solution: 


As before, in this case the variables in the iterated integral are actually independent of each other and 
hence we can integrate each piece and multiply: 


Equation: 
Qn m/2 1 Qn m/2 1 
QQ: ¢ 2 1 21 
psinydpdydé= | dé | sinydy | p*dp = (2n) (1) eae 
0 0 0 0 0 0 


The concept of triple integration in spherical coordinates can be extended to integration over a general 
solid, using the projections onto the coordinate planes. Note that dV and dA mean the increments in 
volume and area, respectively. The variables V and A are used as the variables for integration to express 
the integrals. 


The triple integral of a continuous function f (p, 8, y) over a general solid region 
Equation: 


E = {(,, ”)| (0,9) € D,ui (p, 9) < p < U2 (p, )} 


in R?, where D is the projection of E onto the p6-plane, is 
Equation: 


u2(p,0 


[ff soe [| P roeeaelaa 
# D | u(o,) 


In particular, if D = {(p,)|gi (8) < p < go (9),a < 6 < B}, then we have 
Equation: 


B  92(0) — u2(p,0) 


iI] le / / / f (p, 9, p)p*sin vp dep dp dé. 


a gi(O)  u1(,8) 


Similar formulas occur for projections onto the other coordinate planes. 


Example: 
Exercise: 


Problem: 
Setting up a Triple Integral in Spherical Coordinates 


Set up an integral for the volume of the region bounded by the cone z = ,/3 (x? + y2) and the 
hemisphere z = 1/4 — x? — y? (see the figure below). 
z 


z=\4-x-y 


A region bounded below by a cone and above by a 
hemisphere. 


Solution: 


Using the conversion formulas from rectangular coordinates to spherical coordinates, we have: 


For the cone: z = \/3 (x? + y2) or pcosy = V3psing or tan y = a orp=%- 


For the sphere: z = 4/4 — 22 — y2orz?+27+y? = 4orp? = 4orp=2. 
0=2n y=n7/6 p=2 
Thus, the triple integral for the volume is V (EZ) = il p’sin y dp dy dé. 


6=0 ¢=0 p=0 


Note: 
Exercise: 


Problem: 


Set up a triple integral for the volume of the solid region bounded above by the sphere p = 2 and 
bounded below by the cone y = 7/3. 


Solution: 


viey= f / | esin pd ap a8 


Hint 


Follow the steps of the previous example. 


Example: 
Exercise: 


Problem: 
Interchanging Order of Integration in Spherical Coordinates 


Let & be the region bounded below by the cone z = vy x? + y? and above by the sphere 
z= 27+ y? + 2? ([link]). Set up a triple integral in spherical coordinates and find the volume of the 
region using the following orders of integration: 


a. dp d¢ dé, 
b. dy dp dé. 


Sphere z = x? + y* + 2? 


Cone z = \x? + y? 


A region bounded below by a cone and above by a 
sphere. 


Solution: 


a. Use the conversion formulas to write the equations of the sphere and cone in spherical 


coordinates. 
For the sphere: 
Equation: 
te ee 
p> = pcos 
p = cosy. 
For the cone: 
Equation: 


z= vV2rre+y 
pcosp = V p2sin’y cos? + p2sin’y sin” 


V p’sin”y (cos? + sin*¢) 


pcosp = 

pcospy = psing 
cosy = sing 
a) ae aye 


Hence the integral for the volume of the solid region # becomes 
Equation: 


0=2n y=m/4 p=cosy 
Vie | i p’sin y dp dy dé. 
6=0 y=0 p=0 
b. Consider the yp-plane. Note that the ranges for y and p (from part a.) are 
Equation: 
0<y<77/4 
0<p<cosy. 


The curve p = cos y meets the line y = 7/4 at the point (7/4, v2/2) . Thus, to change the 


order of integration, we need to use two pieces: 
Equation: 


0<p<v2/2 V2/2 


p<l 
0<y<n/4 0 : 


at 
=P cose p: 


Hence the integral for the volume of the solid region & becomes 


Equation: 
6=20 p=V2/2 y=n/4 0=20 p= lero —CoOsmup 
a- f | p’sin y dy dp dé + i, | i p’sin y dy dp dé. 
6= p=0 y=0 C=) paw =v 


In each case, the integration results in V (EB) = 4. 


Before we end this section, we present a couple of examples that can illustrate the conversion from 
rectangular coordinates to cylindrical coordinates and from rectangular coordinates to spherical 
coordinates. 


Example: 
Exercise: 


Problem: 
Converting from Rectangular Coordinates to Cylindrical Coordinates 


Convert the following integral into cylindrical coordinates: 
Equation: 


=1 r=/1-7? z=r/ 2? +y? 


/ / i xyz dz dx dy. 
z=0 


< 


Solution: 


The ranges of the variables are 
Equation: 


AR 
IA 


oy sol 
< 2<Vl-y¥ 


gety? < z2</ar24+y?. 


fo) 
/\ 


The first two inequalities describe the right half of a circle of radius 1. Therefore, the ranges for 0 and 
r are 


Equation: 
1 T 
—— <@< —and0<r<l. 
2 2 
The limits of z are r2 < z< r, hence 
Equation: 
y=1 r= VJ1-y? 2=/a? +y? O=n/2 r=1  z=r 
/ i i xyz dz dx dy = i / r (rcos 8) (r sin 0)z dz dr dé. 
y=-1 z=0 z=a?+y? 6=—1n/2 r=0 z=r? 
Example: 
Exercise: 
Problem: 


Converting from Rectangular Coordinates to Spherical Coordinates 


Convert the following integral into spherical coordinates: 


Equation: 
y=3 w=/9-y? z= 4/18—2?-y? 
i | / (2 ahey ste 2) dz dx dy. 
v0 20 geal yy 
Solution: 


The ranges of the variables are 
Equation: 


IA 


0 y<3 
0 a</9-y 


Jert+y < z<V/18—2?—y?. 


IA 


The first two ranges of variables describe a quarter disk in the first quadrant of the zy-plane. Hence 
the range for 0 isO <0 < F. 


The lower bound z = af x? + y? is the upper half of a cone and the upper bound 


z = \/18 — x? — y? is the upper half of a sphere. Therefore, we have 0 < p < 18, which is 
0<p<3v2. 


For the ranges of y, we need to find where the cone and the sphere intersect, so solve the equation 
Equation: 


oe: lS 
2 
(v2? +9") +27 = 18 
ls 
2 aes 
oo 
p= & 
This gives 
Equation: 
3/2cosy = 3 
=, Ws 
cosp = a 
eae 
Putting this together, we obtain 
Equation: 
y=3 w=/9-y? z= 18-2? -y? y=n/4 0=n/2 p=3V2 
/ / i ie + y+ zy dz dz dy = i / / p'sin y dp dé dy. 
y=0 z=0 z=/erry? y=0 6=0 p=0 
Note: 
Exercise: 
Problem: 


Use rectangular, cylindrical, and spherical coordinates to set up triple integrals for finding the volume 
of the region inside the sphere x? + y” + z? = 4 but outside the cylinder x? + y? = 1. 


Solution: 


a2 y=V 4-2? t= y/4-o w=1 0 y=V 1a? z= 4-0? =? 
Rectangular: / / dz dy dx— i / dz dy dz 
w=-2 yo Vag? z=-y/4-2)-? e=—-1l y-— Vlg? z= /4-2?-y? 
0=2n  r=2 z=V/4—=7? 
Cylindrical: / r dz dr do. 
6=0 r=l 22 /ga7? 
y=5r/6 6=27 = 
Spherical: ih p’sin y dp dé dy. 
y=n/6 6=0 p=cscy 


Now that we are familiar with the spherical coordinate system, let’s find the volume of some known 
geometric figures, such as spheres and ellipsoids. 


Example: 
Exercise: 


Problem: 
Chapter Opener: Finding the Volume of Il’ Hemispheric 


Find the volume of the spherical planetarium in 1’ Hemispheric in Valencia, Spain, which is five 
stories tall and has a radius of approximately 50 ft, using the equation x” + y? + 22 =r? 


(credit: modification of work by 
Javier Yaya Tur, Wikimedia 
Commons) 


Solution: 


We calculate the volume of the ball in the first octant, where x > 0, y > 0, and z > 0, using 
spherical coordinates, and then multiply the result by 8 for symmetry. Since we consider the region D 
as the first octant in the integral, the ranges of the variables are 

Equation: 


0<p<5,0<p<r0<0< 


nw| 3 


Therefore, 
Equation: 


O=n/2  p=n p=n/2 


V = |[fazdyde=8 i / i p’sin 0 dy dp do 
D 6=0 p=0 y=0 


p=n/2 p=r d=1/2 
=8 | dp p'dp sin 6 dé 
y=0 p=0 0=0 
Tv r3 
=8 (5) (S) (1) 
3 


This exactly matches with what we knew. So for a sphere with a radius of approximately 50 ft, the 
volume is 47(50)° ~~ 523,600 ft®. 


For the next example we find the volume of an ellipsoid. 


Example: 
Exercise: 


Problem: 
Finding the Volume of an Ellipsoid 


2 


Find the volume of the ellipsoid _ + - + a = il, 


Solution: 


We again use symmetry and evaluate the volume of the ellipsoid using spherical coordinates. As 
before, we use the first octant z > 0, y > 0, and z > 0 and then multiply the result by 8. 


In this case the ranges of the variables are 
Equation: 


0< pS 7,0<p<>,0<p<1, and0<0< 


w/a 


Also, we need to change the rectangular to spherical coordinates in this way: 
Equation: 


x =apcosypsin#,y = bosin ysin 0, and z = cp cos 6. 


Then the volume of the ellipsoid becomes 


Equation: 
Y= | i} / dx dy dz 
D 
O0=n/2  p=1 y=n/2 
=8 ik / / abcp’sin 6 dy dp dé 
6=0 p=0 y=0 
y=n/2 pai 6=n/2 
= sabe f ap ff odp | sin 0 d0 
y=0 p=0 6=0 
= 8abe (F) (3) (1) 
= + nabe 
Example: 
Exercise: 
Problem: 


Finding the Volume of the Space Inside an Ellipsoid and Outside a Sphere 


22 
2 


90 


Find the volume of the space inside the ellipsoid aS ar & ae = 1 and outside the sphere 


pose see 5 0e 
Solution: 


This problem is directly related to the |’ Hemispheric structure. The volume of space inside the 
ellipsoid and outside the sphere might be useful to find the expense of heating or cooling that space. 
We can use the preceding two examples for the volume of the sphere and ellipsoid and then substract. 


First we find the volume of the ellipsoid using a = 75 ft, b = 80 ft, and c = 90 ft in the result from 


[link]. Hence the volume of the ellipsoid is 
Equation: 


4 
Vaid = 3 7(75)(80)(90) = 2,262,000 ft?. 


From [link], the volume of the sphere is 
Equation: 


Vephere © 523,600 ft?. 


x? y" 2 
757 1 80? + 90? 


Therefore, the volume of the space inside the ellipsoid = 1 and outside the sphere 


a? + y* + z? = 50? is approximately 
Equation: 


Vitemisferic = Vellipsoid = Vephere ar 1,738,400 it. 


Note: 

Hot air balloons 

Hot air ballooning is a relaxing, peaceful pastime that many people enjoy. Many balloonist gatherings take 
place around the world, such as the Albuquerque International Balloon Fiesta. The Albuquerque event is 
the largest hot air balloon festival in the world, with over 500 balloons participating each year. 


Balloons lift off at the 2001 Albuquerque International Balloon Fiesta. (credit: David Herrera, Flickr) 


As the name implies, hot air balloons use hot air to generate lift. (Hot air is less dense than cooler air, so 
the balloon floats as long as the hot air stays hot.) The heat is generated by a propane burner suspended 
below the opening of the basket. Once the balloon takes off, the pilot controls the altitude of the balloon, 
either by using the burner to heat the air and ascend or by using a vent near the top of the balloon to release 
heated air and descend. The pilot has very little control over where the balloon goes, however—balloons 
are at the mercy of the winds. The uncertainty over where we will end up is one of the reasons balloonists 
are attracted to the sport. 

In this project we use triple integrals to learn more about hot air balloons. We model the balloon in two 
pieces. The top of the balloon is modeled by a half sphere of radius 28 feet. The bottom of the balloon is 
modeled by a frustum of a cone (think of an ice cream cone with the pointy end cut off). The radius of the 
large end of the frustum is 28 feet and the radius of the small end of the frustum is 6 feet. A graph of our 
balloon model and a cross-sectional diagram showing the dimensions are shown in the following figure. 


(a) (b) 


(a) Use a half sphere to model the top part of the balloon and a frustum of a cone to model the 
bottom part of the balloon. (b) A cross section of the balloon showing its dimensions. 


We first want to find the volume of the balloon. If we look at the top part and the bottom part of the 
balloon separately, we see that they are geometric solids with known volume formulas. However, it is still 
worthwhile to set up and evaluate the integrals we would need to find the volume. If we calculate the 
volume using integration, we can use the known volume formulas to check our answers. This will help 
ensure that we have the integrals set up correctly for the later, more complicated stages of the project. 


1. Find the volume of the balloon in two ways. 


a. Use triple integrals to calculate the volume. Consider each part of the balloon separately. 
(Consider using spherical coordinates for the top part and cylindrical coordinates for the bottom 


part.) 
b. Verify the answer using the formulas for the volume of a sphere, V = +nrs, and for the volume 


of acone, V = $77r7h. 


In reality, calculating the temperature at a point inside the balloon is a tremendously complicated 
endeavor. In fact, an entire branch of physics (thermodynamics) is devoted to studying heat and 
temperature. For the purposes of this project, however, we are going to make some simplifying 
assumptions about how temperature varies from point to point within the balloon. Assume that just 
prior to liftoff, the temperature (in degrees Fahrenheit) of the air inside the balloon varies according to 
the function 

Equation: 


To(r, 6,2) = = a0e 


N 


. What is the average temperature of the air in the balloon just prior to liftoff? (Again, look at each part 
of the balloon separately, and do not forget to convert the function into spherical coordinates when 
looking at the top part of the balloon.) 


Now the pilot activates the burner for 10 seconds. This action affects the temperature in a 12-foot- 
wide column 20 feet high, directly above the burner. A cross section of the balloon depicting this 
column in shown in the following figure. 


Activating the bummer heats the air in a 20- 
foot-high, 12-foot-wide column directly 
above the burner. 


Assume that after the pilot activates the burner for 10 seconds, the temperature of the air in the 
column described above increases according to the formula 
Equation: 


H(r, 6, z) = —2z — 48. 


Then the temperature of the air in the column is given by 
Equation: 


Ti(r,0,z) = = “210-1 (= 22 — 45), 


while the temperature in the remainder of the balloon is still given by 
Equation: 


To(r, 6,2) = ae ce0s 


ice) 


. Find the average temperature of the air in the balloon after the pilot has activated the burner for 10 
seconds. 


Key Concepts 


e To evaluate a triple integral in cylindrical coordinates, use the iterated integral 
Equation: 


0=B8 r=g2(0) z=ua(r,0) 
f (r, 9, z)r dz dr dé. 


0=a r=gi(0)  z=u1(r,0) 


e To evaluate a triple integral in spherical coordinates, use the iterated integral 
Equation: 


0=B p=g92(9) p=ua(r,8) 
f (0,9, 9) p’ sin p dip dp dé. 


O=a p=g9i(9) y=ui(r,0) 


Key Equations 


¢ Triple integral in cylindrical coordinates 


Gane = | [ai aera naia | Gerais 


¢ Triple integral in spherical coordinates 
g=$ 0=8 p=b 


Il f (p,9, )p?sin y dp dy dO = / i / f (p,8, y)p*sin v dp dy dé 
B ~=7 b=a p=a 


In the following exercises, evaluate the triple integrals ‘ (i (i f(x, y, z)dV over the solid EB. 
E 


Exercise: 


Problem: f(z, y, 2) = z, B= {(2,y,2)|@’? +y? <9,2>0,y>0,0<2z<1} 


Exercise: 


Problem: f(z, y, 2) = x2’, B= {(a,y,z)|2’?>+y’ < 16,2 >0,y<0,-1<z< 1} 


Exercise: 


Problem: F(z, y, 2) =ay,B= {(@,9)2) 2? ey? a Le¢20,¢2y-lezcs 1} 


Solution: 


ool 


Exercise: 


Problem: f(z, y,z) = «7+ y’, B= {(x,y,z)|a? +? eae 0.290 a2 3} 


Exercise: 
Problem: f(x,y, z) =eV" ‘”, B= {(,y,2)|1 <a ty <4,y<0,¢ <yV3,2<2< 3} 


Solution: 


me” 


6 


Exercise: 


Problem: f(x,y, z) = \/z? + y?, B= { (oy, 2) [1 <ety<9,y<0,0<z<1} 
Exercise: 


Problem: 


a. Let B be a cylindrical shell with inner radius a, outer radius b, and height c, where 0 <a < 6 
and c > 0. Assume that a function F’ defined on B can be expressed in cylindrical coordinates as 


b 
F (x,y,z) = f(r) + h(2), where f and hare differentiable functions. If ‘| f(r)dr = 0 and 


h(0) = 0, where f and hare antiderivatives of f and h, respectively, show that 
Equation: 


[ff Peay = 2 (70) ~afla)) +26 able 


b. Use the previous result to show that II] (2 + siny/a2 + v) dx dy dz = 6n? (x — 2), where 
B 
Bisa cylindrical shell with inner radius 7, outer radius 27, and height 2. 


Exercise: 


Problem: 


a. Let B be a cylindrical shell with inner radius a, outer radius 6, and height c, where 0 < a < b 
and c > 0. Assume that a function Ff defined on B can be expressed in cylindrical coordinates as 


b 
F(a, y,z) = f(r)g(@)h(z), where f, 9, and h are differentiable functions. If [Fever = 0; 


where f is an antiderivative of f, show that 
Equation: 


[[f F @u2av = [670) — aFla)] (927) — 910) [ile) - 0), 


where g and A are antiderivatives of g and h, respectively. 


b. Use the previous result to show that | i} / zsiny/22 + y2dz dy dz = —12n?, where Bisa 


cylindrical shell with inner radius 7, outer radius 27, and height 2. 


In the following exercises, the boundaries of the solid & are given in cylindrical coordinates. 


a. Express the region & in cylindrical coordinates. 


b. Convert the integral f (x,y, z)dV to cylindrical coordinates. 


Exercise: 


Problem: 


E is bounded by the right circular cylinder r = 4 sin 6, the r@-plane, and the sphere r? + z? = 16. 


Solution: 


a B= {(r,0,2)|0<0< 0,0 <r<4sin8,0<2< VI6—r hb, 
a 4sin0 V16—r? 
f (r, 9, z)r dz dr d0 
0 60 0 
Exercise: 


Problem: 


E is bounded by the right circular cylinder r = cos 6, the r@-plane, and the sphere r? + z? = 9. 
Exercise: 


Problem: 


E is located in the first octant and is bounded by the circular paraboloid z = 9 — 3r?, the cylinder 
r = V3, and the plane r (cos 6 + sin 0) = 20 — z. 


Solution: 


= { (7,4, 2)]0 << 7,027S V3,9—-r<z< 10 ~ r (cos 6 + sin 6) }; b. 
1/2 4/3 10—r(cos 6+sin @) 


| / i Ae Bn dedr ab 


Exercise: 


Problem: 


Eis located in the first octant outside the circular paraboloid z = 10 — 2r? and inside the cylinder 
r = V5 and is bounded also by the planes z = 20 and 6 = 7. 


In the following exercises, the function f and region F are given. 


a. Express the region F and the function f in cylindrical coordinates. 


b. Convert the integral f (x,y, z)dV into cylindrical coordinates and evaluate it. 


Exercise: 


Problem: f (x, y, 2) = a, B= {(z,y,2)|0 <a? +y<9,2>0,y>0,0<z2<2+3} 
Solution: 


a. B= {(r,0,2)|0 <7 <3,0<0< $,0<2z<rcosd +3}, f(r,0,2) = 
3 1/2 rcos6+3 


r On 
ap : roi ae 
0 0 0 


1 : 
rcos0+3? b. 


Exercise: 


Problem: f (x,y,z) =a? +y’, B= {(2,y,2)|0< ae +y? <4,y>0,0<z< 3-2} 


Exercise: 


Problem: f (x,y,z) =2,E= {(z,y,2|1<yt+2<90<2<1-y-2} 
Solution: 


rcos#,z=rsin@,x = 2, 
(r,0,2)|L<r<3,0<0<27,0<2< ae? pF. 2) = ee be 


1-r? 
/ zrdzdédr = a 


0 0 


ll 
ys Il 


ee ese 


Exercise: 
Problem: f (x,y,z) = y, E= {(x, y, 2)|1 <a? 4+27<90<y<1-2 — 27} 


In the following exercises, find the volume of the solid F whose boundaries are given in rectangular 
coordinates. 


Exercise: 


Problem: £ is above the xy-plane, inside the cylinder x? + y? = 1, and below the plane z = 1. 


Solution: 


Tv 


Exercise: 


Problem: £ is below the plane z = 1 and inside the paraboloid z = x? + y?. 


Exercise: 


Problem: E is bounded by the circular cone z = \/x? + y? and z = 1. 


Solution: 


Coe) 


Exercise: 
Problem: 
F is located above the xy-plane, below z = 1, outside the one-sheeted hyperboloid 
x? + y? — z? = 1, and inside the cylinder x? + y? = 2. 

Exercise: 


Problem: 


2 


E is located inside the cylinder x” + y? = 1 and between the circular paraboloids z = 1 — x? — y? 


and z= a? + y?. 
Solution: 


TV 
Exercise: 


Problem: 


E is located inside the sphere x? + y? + z? = 1, above the xy-plane, and inside the circular cone 
z=vVarr2+y?. 
Exercise: 


Problem: 
E is located outside the circular cone x? + y? = (z— 1)? and between the planes z = 0 and z = 2. 


Solution: 


An 


3 


Exercise: 


Problem: 


F is located outside the circular cone z = 1 — J/ x? + y*, above the zy-plane, below the circular 
paraboloid, and between the planes z = 0 and z = 2. 


Exercise: 


Problem: 
[T] Use a computer algebra system (CAS) to graph the solid whose volume is given by the iterated 
m/2 1 


integral in cylindrical coordinates / / / r dz dr d@. Find the volume V of the solid. Round your 


—n/2 0 r 
answer to four decimal places. 


Solution: 


Exercise: 


Problem: 


[T] Use a CAS to graph the solid whose volume is given by the iterated integral in cylindrical 


m/2 1 or 
coordinates : / r dz dr dé. Find the volume V of the solid Round your answer to four decimal 
0 0 rf 
places. 
Exercise: 
Problem: 


1 Vi-y Var+y? 
Convert the integral / / i xz dz dz dy into an integral in cylindrical coordinates. 


0 —/1-¥ x+y? 


Solution: 


17 fF 
/ 7 i zr’cos 6 dz dé dr 
0 O Fr 


Exercise: 


2 #« it 
Problem: Convert the integral / / / (ay + z)dz dz dy into an integral in cylindrical coordinates. 
0 0 0 


In the following exercises, evaluate the triple integral i} (i (i f(a, y, z)dV over the solid B. 
B 
Exercise: 


Problem: f(z, y, z) = 1, B= {(z,y, z) |x? +4? 4+ 27 < 90,2> o} 


Solution: 


1807/10 


Exercise: 


Problem: f(z, y, z) =1—/2?+y?+22,B= {(x,y,2)|a? +4? +22 <9,y>0,2>0} 


Exercise: 


Problem: 


f(x,y, z) = \/x? + y?, B is bounded above by the half-sphere x? + y? + z? = 9 with z > 0 and 
below by the cone 22” = x? + y’. 
Solution: 


817(1—2) 
16 


Exercise: 
Problem: 
f(x, y, z) = 2, Bis bounded above by the half-sphere x” + y? + 27 = 16 with z > 0 and below by 
the cone 2z” = a2? + y’. 
Exercise: 
Problem: 


Show that if F'(p,0, ¢) = f (p)g (@)h (¢) is a continuous function on the spherical box 
B={(p,6,y)la<p<ba<O0<B,y<9 < }, then 


Equation: 
b B y 
[[fev=| [eroa)| | foe} | [nosing ae 
B a Qa Y 
Exercise: 
Problem: 


a. A function F is said to have spherical symmetry if it depends on the distance to the origin only, 
that is, it can be expressed in spherical coordinates as F'(x, y,z) = f(p), where 


p=v/ «2? + y? + z?. Show that 


Equation: 
b 
Il F(z,y, z)dV oo an f 0° Fo)de, 
B a 


where B is the region between the upper concentric hemispheres of radii a and 6 centered at the 
origin, with 0 < a < band F'a spherical function defined on B. 


b. Use the previous result to show that Il (e + y" + z*) Vx? + y? + 22 dV = 217, where 
B 
Equation: 


B= {(2,y,2)|1 <a? +y +2 <2,2>0}. 


Exercise: 


Problem: 


a. Let B be the region between the upper concentric hemispheres of radii a and b centered at the 
origin and situated in the first octant, where 0 < a < b. Consider F a function defined on B 
whose form in spherical coordinates (p, 0, y) is F (x,y, z) = f(p)cos y. Show that if 

b 


ageySoene fe pd 2 O:tieh 


I[[ Fe v,y2)4V =“ [ah(a) — bh(0) 


where g is an antiderivative of f and h is an antiderivative of g. 


Equation: 


Jety+2 2 
b. Use the previous result to show that / ii / ee Veter? wae dV = =, where B is the 
far ty? + 2? 


region between the upper concentric nena iene of radii 7 and 27 centered at the origin and 
situated in the first octant. 


In the following exercises, the function f and region E are given. 


a. Express the region F and function f in cylindrical coordinates. 


b. Convert the integral | i i f (x,y, z)dV into cylindrical coordinates and evaluate it. 


Exercise: 


Problem: f (x,y,z) = z; E= {(zx, y, 2) |0 <etyt+2<1,z>0} 


Exercise: 
Problem: f (x,y,z) =zr+y;E= { (x,y, 2) [1 Soo ep te a 9.2 > 0 > 0} 


Solution: 


a. f (p,9, 9) = psiny (cos@ + sin0), E = {(p,0,y)|1< p< 2,0<0<7,0<9K< Fh5b. 
an m/2 2 
3 : — 157 
j [e cos y sin y dp dip dO = =" 
0 0 1 


Exercise: 


Problem: f (x, y, z) = 2ry; E = {(2,y,2)| Va? +9? <z<vVl1l-2?-y,¢>0,y> 0} 


Exercise: 


Problem: f (x, y, z) = z; H = { (x,y, 2)? +y? +27 -22< 0,22 +y? < z} 


Solution: 


a. f (9,9, 9) = pcos y; E = {(p,9, )|0 < p< 2cosy,0<0< F,0<y< F}sb. 
m/2 1/4 2cosp 


psin ycos y dp dy dé = — 
0.2 <0. <0 


In the following exercises, find the volume of the solid F whose boundaries are given in rectangular 
coordinates. 
Exercise: 


Problem: £ = {(2.y,2)| Va +H <z2< J/16-2?—y?,2>0,y> 0} 


Exercise: 


Problem: E = {(2.y,2)|2” +y° + y —2z<0, V ae? + y? Ss z} 


Solution: 


ENE) 


Exercise: 


Problem: 


Use spherical coordinates to find the volume of the solid situated outside the sphere p = 1 and inside 
the sphere p = cos y, with y € [0, =. 


Exercise: 


Problem: 


Use spherical coordinates to find the volume of the ball » < 3 that is situated between the cones 


p=Fandy= fz. 


Solution: 
On (v2 _ 1) 


Exercise: 


Problem: 


4 /16—y? VJ/16—22-¥? 


Convert the integral (az ty? + z”) dz dz dy into an integral in spherical 


-4 /16 y? 4/16 x2—y? 


coordinates. 
Exercise: 


Problem: 


4 Vi6—2? V/16—a?-y? 
Convert the integral i, i / (x? Sea z’) dz dy dz into an integral in spherical 
eo —/16—2?—y? 


coordinates. 


Solution: 


w/2 m/2 4 
/ / [efi y dp dy dé 
0 0 0 


Exercise: 


Problem: 


2 Via? 16-2? 
Convert the integral dz dy dz into an integral in spherical coordinates and 
2 yee Vee 
evaluate it. 
Exercise: 


Problem: 


[T] Use a CAS to graph the solid whose volume is given by the iterated integral in spherical 


xn 7/6 2 
coordinates : i i. p’sin y dp dy dé. Find the volume V of the solid. Round your answer to 
m/2 57/6 0 


three decimal places. 


Solution: 


V = 408 ~ 7.255 


Exercise: 
Problem: 
[T] Use a CAS to graph the solid whose volume is given by the iterated integral in spherical 
2r o/4 
coordinates as / / / p’sin y dp dy dé. Find the volume V of the solid. Round your answer to 


0 37/4 0 
three decimal places. 


Exercise: 


Problem: 


[T] Use a CAS to evaluate the integral i (i / (a? +y’ )dV where F& lies above the paraboloid 
E 
z= a«* +? and below the plane z = 3y. 


Solution: 


3437 
32 


Exercise: 
Problem: [T] 
a. Evaluate the integral i (} evar dV, where EF is bounded by the spheres 
E 


4a? + 4y? + 42? = anda? +y?4 2? =1. 


b. Use a CAS to find an approximation of the previous integral. Round your answer to two decimal 
places. 


Exercise: 
Problem: 


Express the volume of the solid inside the sphere x? + y? + z? = 16 and outside the cylinder 
az? + y” = 4as triple integrals in cylindrical coordinates and spherical coordinates, respectively. 


Solution: 

2 4 Vi6—r? 5n/6 Qn 

ey i rdedras; / | e%sin pdp d9 dg 
0 2." Sige ee n/6 0 2cscp 
Exercise: 

Problem: 


Express the volume of the solid inside the sphere x? + y? + z? = 16 and outside the cylinder 
a? + y” = 4 that is located in the first octant as triple integrals in cylindrical coordinates and spherical 
coordinates, respectively. 


Exercise: 


Problem: 


The power emitted by an antenna has a power density per unit volume given in spherical coordinates 
by 


P(p, 9, 9) = rcos"6 sin*y, 


where Pp is a constant with units in watts. The total power within a sphere B of radius r meters is 


defined as P = |[f- (p, 9, y)dV. Find the total power P. 
B 


Solution: 


P= 22Ror. watts 


Exercise: 


Problem: 


Use the preceding exercise to find the total power within a sphere B of radius 5 meters when the 
power density per unit volume is given by p(p, 0, y) = 37 cos” sin*y. 


Exercise: 


Problem: 


A charge cloud contained in a sphere B of radius r centimeters centered at the origin has its charge 


density given by q (a, y, z) = ky/x? + y? + 2? ue , where k > 0. The total charge contained in B is 


given by Q = I q(x, y, z)dV. Find the total charge Q. 
B 


Solution: 

Q = kr*mpC 
Exercise: 

Problem: 


Use the preceding exercise to find the total charge cloud contained in the unit sphere if the charge 
density is q (2, y, z) = 20./x? + y? + pee 


Glossary 


triple integral in cylindrical coordinates 
the limit of a triple Riemann sum, provided the following limit exists: 
Equation: 


l min 

* * * * 
1 jim | S Frisk: Os its ZijK)T ijghATAOAz 
as i=1 j=1 k=1 


triple integral in spherical coordinates 
the limit of a triple Riemann sum, provided the following limit exists: 
Equation: 


l m 
‘ * * * * 2 . 
lim » Ss f(Pijns Dies Pizr) (ijn) sin pApAbAy 


Calculating Centers of Mass and Moments of Inertia 


e Use double integrals to locate the center of mass of a two-dimensional object. 
e Use double integrals to find the moment of inertia of a two-dimensional object. 
e Use triple integrals to locate the center of mass of a three-dimensional object. 


We have already discussed a few applications of multiple integrals, such as finding areas, volumes, and the 
average value of a function over a bounded region. In this section we develop computational techniques for 
finding the center of mass and moments of inertia of several types of physical objects, using double integrals 
for a lamina (flat plate) and triple integrals for a three-dimensional object with variable density. The density 
is usually considered to be a constant number when the lamina or the object is homogeneous; that is, the 
object has uniform density. 


Center of Mass in Two Dimensions 


The center of mass is also known as the center of gravity if the object is in a uniform gravitational field. If the 
object has uniform density, the center of mass is the geometric center of the object, which is called the 
centroid. [link] shows a point P as the center of mass of a lamina. The lamina is perfectly balanced about its 
center of mass. 


A lamina is perfectly balanced on a spindle if the 
lamina’s center of mass sits on the spindle. 


To find the coordinates of the center of mass P(z, y) of a lamina, we need to find the moment M, of the 
lamina about the z-axis and the moment M, about the y-axis. We also need to find the mass m of the 
lamina. Then 

Equation: 


_ M, _ 
x= —andy = 
m 


Refer to Moments and Centers of Mass for the definitions and the methods of single integration to find the 
center of mass of a one-dimensional object (for example, a thin rod). We are going to use a similar idea here 
except that the object is a two-dimensional lamina and we use a double integral. 


M, 
m 


If we allow a constant density function, then z = Ms and y = give the centroid of the lamina. 


Suppose that the lamina occupies a region R in the zy-plane, and let p (z, y) be its density (in units of mass 
per unit area) at any point (a, y). Hence, p(x, y) = Jim, oe where Am and AA are the mass and area of 
> 


a small rectangle containing the point (z, y) and the limit is taken as the dimensions of the rectangle go to 0 
(see the following figure). 


(x, Y) 


Lamina 


The density of a lamina at a point is the limit of its 
mass per area in a small rectangle about the point 
as the area goes to zero. 


* * 
Just as before, we divide the region R into tiny rectangles R;; with area A.A and choose (x5, Vij) as sample 


points. Then the mass m,; of each R;; is equal to p (2;,, y;;) AA ([link]). Let & and I be the number of 


subintervals in x and y, respectively. Also, note that the shape might not always be rectangular but the limit 
works anyway, as seen in previous sections. 


(xj. Yi) Ri 


Subdividing the lamina into tiny rectangles R;;, 
* * 


each containing a sample point (z;;, y;;). 


Hence, the mass of the lamina is 
Equation: 


kool k ol 
m= fim Dm = in Dents )OA= ff oa 
q=1 7=1 i=l 9=1 
Let’s see an example now of finding the total mass of a triangular lamina. 


Example: 
Exercise: 


Problem: 
Finding the Total Mass of a Lamina 


Consider a triangular lamina R with vertices (0, 0), (0,3), (3,0) and with density 
Ce, y)—2y kg/m’. Find the total mass. 


Solution: 


A sketch of the region R is always helpful, as shown in the following figure. 


x+y=3 


dm = p(x, y)dA 


A lamina in the zy-plane with density 
p(z,y) = zy. 


Using the expression developed for mass, we see that 


Equation: 
2=3 y=3-zr p=8) 
y? 
[fam |f elena = / / cy dydz = / Ee 
R z=0 y=0 0 
3 


3 
| 


R 


1 =) 


2 _ | 9x? 3 De 
5 2(3 — 2) dx = |e -2 += 


I 


The computation is straightforward, giving the answer m = aL kg. 


Note: 
Exercise: 


Problem: 


Consider the same region R as in the previous example, and use the density function p(x, y) = ,/ry. 
Find the total mass. 


Solution: 


is 
“3 kg 


Now that we have established the expression for mass, we have the tools we need for calculating moments 
and centers of mass. The moment M,, about the x-axis for R is the limit of the sums of moments of the 
regions R,; about the z-axis. Hence 

Equation: 


Ma= mn 99> (H,)m= im D292 (7 )o (ones) a= [[ ween 


Similarly, the moment M, about the y-axis for R is the limit of the sums of moments of the regions R,; 
about the y-axis. Hence 


Equation: 
k ol ‘ k ol 
My = fim, 2297 (25) = tie, 29 (se (Poss) O4= ff soon 
Example: 
Exercise: 
Problem: 


Finding Moments 


Consider the same triangular lamina R with vertices (0,0), (0,3), (3,0) and with density 
p(x, y) = xy. Find the moments M, and My. 


Solution: 


Use double integrals for each moment and compute their values: 


Equation: 
2=3 y=3—-2 
9 81 
IM. = yp (z,y)dA = zy dy dz = 30” 
R e=0 G0 
Equation: 
z=3 y=3—-2 
9 81 
i zp(z,y)dA = x ydy dz = 5. 
R z=0 70 


The computation is quite straightforward. 


Note: 
Exercise: 


Problem: 


Consider the same lamina R as above, and use the density function p(x, y) = ,/xy. Find the moments 
M, and My. 


Solution: 


M, = 5 and M, = 3 


Finally we are ready to restate the expressions for the center of mass in terms of integrals. We denote the x- 


coordinate of the center of mass by x and the y-coordinate by y. Specifically, 


Equation: 
[[ro(waa [[wvenaa 
_ M - M, 
g=— = and y = 
m m 
[[eocouaa [[ecuaa 
R R 
Example: 
Exercise: 
Problem: 


Finding the Center of Mass 


Again consider the same triangular region R with vertices (0,0), (0,3), (3,0) and with density 
function p (x, y) = xy. Find the center of mass. 


Solution: 


Using the formulas we developed, we have 


Equation: 
fro z,y)dA 
- M, _ 81/20 6 
m NG 8 ie 
i i p(x, y)dA / 
R 
Equation: 


Therefore, the center of mass is the point (g, £), 
Analysis 


If we choose the density p (2, y) instead to be uniform throughout the region (i.e., constant), such as the 
value 1 (any constant will do), then we can compute the centroid, 
Equation: 


5? 9) is not exactly the same as the centroid (1, 1) of the triangular region. 
This is due to the variable density of R. If the density is constant, then we just use p (x, y) = c (constant). 
This value cancels out from the formulas, so for a constant density, the center of mass coincides with the 
centroid of the lamina. 


Notice that the center of mass (t 5 


Note: 
Exercise: 


Problem: 


Again use the same region R as above and the density function p (x, y) = ,/zy. Find the center of 
mass. 


Solution: 
—- My, _ 8in/64_ — _M, _ 81n/64 _ 9 
Cee Said 9n/8 = 3 andy = m 9n/8 8° 


Once again, based on the comments at the end of [link], we have expressions for the centroid of a region on 


the plane: 
Equation: 

ut fl xcdA wl / ydA 
and yo = 
ue 7 fl dA ™ fu faa dA 


We should use these formulas and verify the centroid of the triangular region R referred to in the last three 
examples. 


Le = 


Example: 
Exercise: 


Problem: 
Finding Mass, Moments, and Center of Mass 


Find the mass, moments, and the center of mass of the lamina of density p (x, y) = x + y occupying 


the region R under the curve y = 2” in the interval 0 < x < 2 (see the following figure). 
y 


Locating the center of mass of a 
lamina R with density 


p(z,y)=a2t+y. 


Solution: 


First we compute the mass m. We need to describe the region between the graph of y = x? and the 
vertical lines = 0 and x = 2: 


Equation: 
z=2 y=2? 2=2 2 y=2? 
m ce ae ff @+viyar= f Jev+¥ |i 
z=0 y=0 2=0 TM 
- 4 5] |e=2 
z=0 


Now compute the moments M, and M,: 
Equation: 


M, 


Equation: 


Finally, evaluate the center of mass, 
Equation: 


80 
y (a + y)dy dx = ae 


176 
dy dx = —. 
x(x + y)dydz ie 


[[eoemaa 


Pi M, _ 176/15 44 
Siri = ays — ore 
[[ocemaa / 
R 
[[ veya 
eo ie _, Sf _ i100 
US a = 36/5 — 63° 


Hence the center of mass is (x, y) = ( a 


Note: 
Exercise: 


Problem: 


Calculate the mass, moments, and the center of mass of the region between the curves y = x and 


44 100 
27? 63 


I play)dA 


i 


y = x’ with the density function p (a, y) = a in the interval 0 < x <1. 


Solution: 


Example: 
Exercise: 


Problem: 
Finding a Centroid 


Find the centroid of the region under the curve y = e” over the interval 1 < x < 3 (see the following 


figure). 


Finding a centroid of a region below the 
curve y = e”. 


Solution: 


To compute the centroid, we assume that the density function is constant and hence it cancels out: 
Equation: 


3|5 
r 
8 
= 
= 
< 
= 


c— = and y, Ms _ 
m 
// dA // dA 
z=3 y=e” 2=3 
[faa / / a dy dx / ze*dx 
Re = My, _ oR _ #=1_y=0 _ a= Je? 2e? 
Cc =, t= es = 2=3 a re ~~ PST? 
m [fea 3 y e'—e e2—1 
/ i dy dz e*dx 
oly) g=1 
z=3 y=e” 2=3 
2a 
[fvea / i) y dy dx ele 
12(.4 
_ Mrz _ RB __ 2=1_ y=0 meee _ (e4-1) _i 9 
c= a = $=0 yee" “  2=3 e(e?—1) ae (e + 1) 
dA 
/ dy dz e*dx 
cl 0: g=1 


Thus the centroid of the region is 
Equation: 


Yee «AN 
(Ze, Ye) = (= = ee (e ar i)). 


Note: 
Exercise: 


Problem: 


Calculate the centroid of the region between the curves y = az and y = ./z with uniform density in the 
interval0 < a2 <1. 


Solution: 
My, _ 1/18 _ 2 Mz _ 1/2 _ 4 
Vee a 1/6 5 and ye = =, 1/6 2 


Moments of Inertia 


For a clear understanding of how to calculate moments of inertia using double integrals, we need to go back 


to the general definition in Section 6.6. The moment of inertia of a particle of mass m about an axis is mr?, 


where r is the distance of the particle from the axis. We can see from [link] that the moment of inertia of the 
subrectangle R;; about the x-axis is (y;;)° p(X;;, yi;) MA. Similarly, the moment of inertia of the 


rae *\2 ‘ io . 
subrectangle R;; about the y-axis is (z,;) p(x;;, y;;) 4A. The moment of inertia is related to the rotation of 


the mass; specifically, it measures the tendency of the mass to resist a change in rotational motion about an 
axis. 


The moment of inertia 7, about the z-axis for the region R is the limit of the sum of moments of inertia of 
the regions R;; about the x-axis. Hence 
Equation: 


2 
I, = tim 3 (v:;) Miz = tim > > (v.,) p (2is.443) AA= i / yp (x, y)dA. 
R 


Similarly, the moment of inertia I, about the y-axis for FR is the limit of the sum of moments of inertia of the 
regions R,,; about the y-axis. Hence 
Equation: 


Sometimes, we need to find the moment of inertia of an object about the origin, which is known as the polar 
moment of inertia. We denote this by Jp and obtain it by adding the moments of inertia J, and I. Hence 
Equation: 


igp=i ed = ff + v)ole.waa. 
R 


All these expressions can be written in polar coordinates by substituting « = r cos 0, y = rsin @, and 
dA = r dr d@. For example, Ip = // r’p(rcos 6,7 sin 0)dA. 
R 


Example: 
Exercise: 


Problem: 
Finding Moments of Inertia for a Triangular Lamina 


Use the triangular region R with vertices (0, 0), (2, 2), and (2,0) and with density p(x, y) = ryas in 
previous examples. Find the moments of inertia. 


Solution: 


Using the expressions established above for the moments of inertia, we have 


Equation: 
a |[vote.naa= / | evtdyde = §, 
R z=0 y= 
C2 = 
a |[eeeaa = i | Pudyde = 3, 
R z=0 y=0 
7 iQ 
i — // (0? + y*)ole,y)da= ff (2? +v)aydyds 
R 0 0 
= I,+1,=8. 
Note: 
Exercise: 
Problem: 


Again use the same region R as above and the density function p(x, y) = ,/zy. Find the moments of 
inertia. 


Solution: 
a=2 y=" a=2 y=" 
z=0 y=0 z=0 y=0 
r=2 y=a 
iy (a? + y*),/aydy dx = 8. 


8 
ll 
= 
4 
ll 
i= 


As mentioned earlier, the moment of inertia of a particle of mass m about an axis is mr? where r is the 
distance of the particle from the axis, also known as the radius of gyration. 


Hence the radii of gyration with respect to the z-axis, the y-axis, and the origin are 


Equation: 
Ty [I Ty 
eg ee —¥ and Ry = 4/—, 
m m m 


respectively. In each case, the radius of gyration tells us how far (perpendicular distance) from the axis of 
rotation the entire mass of an object might be concentrated. The moments of an object are useful for finding 
information on the balance and torque of the object about an axis, but radii of gyration are used to describe 
the distribution of mass around its centroidal axis. There are many applications in engineering and physics. 
Sometimes it is necessary to find the radius of gyration, as in the next example. 


Example: 
Exercise: 


Problem: 
Finding the Radius of Gyration for a Triangular Lamina 


Consider the same triangular lamina R with vertices (0,0), (2, 2), and (2,0) and with density 
p(x, y) = xyas in previous examples. Find the radii of gyration with respect to the x-axis, the y-axis, 
and the origin. 


Solution: 
If we compute the mass of this region we find that m = 2. We found the moments of inertia of this 


lamina in [link]. From these data, the radii of gyration with respect to the x-axis, y-axis, and the origin 
are, respectively, 


Equation: 
R [Ie a= = (2-48 
16/3 2/6 
n= fe-\E=f8-¥ 
Note: 
Exercise: 
Problem: 


Use the same region R from [link] and the density function p (x, y) = ,/zy. Find the radii of gyration 
with respect to the z-axis, the y-axis, and the origin. 


Solution: 


_ 6/35 _ 6/15 _ 4/2 
Bett ar aN nS 


Hint 


Follow the steps shown in the previous example. 


Center of Mass and Moments of Inertia in Three Dimensions 


All the expressions of double integrals discussed so far can be modified to become triple integrals. 


Note: 
Definition 
If we have a solid object Q with a density function p (2, y, z) at any point (x, y, z) in space, then its mass is 


Equation: 
m= fffoeusav 
Q 


Its moments about the zy-plane, the xz-plane, and the yz-plane are 


Equation: 
Moy = [fou z)dV, M,, = [ff (x, y, z)dV, 
Q Q 
Ve = /I/ rp (x,y, z)dV. 
Q 


If the center of mass of the object is the point (2, Y; 2) , then 


Equation: 


Also, if the solid object is homogeneous (with constant density), then the center of mass becomes the 
centroid of the solid. Finally, the moments of inertia about the yz-plane, the xz-plane, and the ry-plane 


joe /If (y? + 27)o(2,y, 2)dV, 
Q 


Equation: 
i= /If (a? + 27) p(z,y, z)dV, 
Q 


— /If (x? + y*)p (zx, y, z)dV. 
Q 


Example: 
Exercise: 


Problem: 
Finding the Mass of a Solid 


Suppose that Q is a solid region bounded by x + 2y + 3z = 6 and the coordinate planes and has 
density p (x,y,z) = xyz. Find the total mass. 


Solution: 


The region Q is a tetrahedron ([link]) meeting the axes at the points (6, 0,0), (0, 3,0), and (0, 0,2). To 


find the limits of integration, let z = 0 in the slanted plane z = + (6 — g — 2y). Then for x and y find 


the projection of Q onto the zy-plane, which is bounded by the axes and the line x + 2y = 6. Hence 
the mass is 


Equation: 


z=6 y=1/2(6—-x) z=1/3(6—x—2y) 
108 


m= fff o(e.u.2av = i i / xyz dzdydx = a 
Q x 


=0 y=0 z=0 


= 3.086. 


Zz 


xX+2y+3z=6 


Finding the mass of a three- 
dimensional solid Q. 


Note: 
Exercise: 


Problem: 
Consider the same region Q ([link]), and use the density function p (2, y, z) = xy*z. Find the mass. 
Solution: 


54 __ 
54 — 1.543 


Hint 


Follow the steps in the previous example. 


Example: 
Exercise: 


Problem: 
Finding the Center of Mass of a Solid 


Suppose Q is a solid region bounded by the plane x + 2y + 3z = 6 and the coordinate planes with 
density p(x, y, z) = xyz (see [link]). Find the center of mass using decimal approximation. 


Solution: 


We have used this tetrahedron before and know the limits of integration, so we can proceed to the 


computations right away. First, we need to find the moments about the zy-plane, the xz-plane, and the 


yz-plane: 
Equation: 


z=6 y=1/2(6—z) 


May = [ff 20 (2.v,24V = 
Q 


z=0 


y=0 


z=1/3(6—z—2y) 
4 
a yz*dz dy dx = = = 1.543, 


z=0 


z=6 y=1/2(6—z) z=1/3(6—x—2y) 


Maz = [ff ve(e.u2av = / 
Q 


z=0 


y=0 


81 
xyz dz dy dz = ne w 2.314, 


2=0 


2z=6 y=1/2(6—z) z=1/3(6—x—2y) 


My. = fff 20 (2,y,24Vv = / 
Q x 


243 
a> yz dz dy dz = ae = 6.943. 


=0 y=0 z=0 
Hence the center of mass is 
Equation: 
ae My. Ts Mz, ee Myy 
| m a m z= m? 
= M, 243/35 243 
— ye — — 
i= f0s/35 «08 2.25, 
5 ey, — Se _ gt — 
Y= m= 108/35 — 108 = 0.75, 
= mM 54/35 54 
= Ly — — 
a= a = 10s/as 7 108 0.5. 


The center of mass for the tetrahedron Q is the point (2.25, 0.75, 0.5). 


Note: 


Exercise: 


Problem: 


Consider the same region Q ([link]) and use the density function p (x, y, z) = ry*z. Find the center of 
mass. 


Solution: 
3 ® il 
(ee) 
Hint 


Check that M,, = a. Me — — and M,, = = Then use m from a previous checkpoint question. 


We conclude this section with an example of finding moments of inertia J,, I, and I,. 


Example: 
Exercise: 


Problem: 
Finding the Moments of Inertia of a Solid 


Suppose that Q is a solid region and is bounded by x + 2y + 3z = 6 and the coordinate planes with 
density p(x, y, 2) = xyz (see [link]). Find the moments of inertia of the tetrahedron Q about the 
yz-plane, the xz-plane, and the ry-plane. 


Solution: 
Once again, we can almost immediately write the limits of integration and hence we can quickly 


proceed to evaluating the moments of inertia. Using the formula stated before, the moments of inertia 
of the tetrahedron @ about the xy-plane, the xz-plane, and the yz-plane are 


Equation: 
I, = II (y? Se 2’) p (x,y, z)dV, 
Q 
i= /I/ (x? =p 2 )p(a,y, z)dV, 
Q 
and 
Equation: 


I,= Ill (x? + y”)p (x, y, z)dV with p (a, y, z) = x7 yz. 
Q 


Proceeding with the computations, we have 
Equation: 


z=6 y=7(6—-2) 2=3(6—-2—2y) 


I, = I] (y? + z’)xyzdV = i / i (y? cc 2) a yz dz dy dx a alti ~ 3.343, 
Q xz=0 y=0 z=0 


z=6 y=4(6—-2) 2=%(6—2-2y) 


j= Il G. = 2) x? yz aV = / il / («? a 2) xyz dz dy dx = 684 ~ 19.543, 
Q xz=0 y=0 z=0 


a=6 y=z(6—a) 2=3(6—2—2y) 


I, = /I/ (ee = y’) a? yz dV = / / / (x? & y?) 2? yz dz dy dx = ie) ~ 20.829. 
Q z=0 y=0 z=0 


Thus, the moments of inertia of the tetrahedron Q about the yz-plane, the xz-plane, and the ry-plane 
are 117/35, 684/35, and 729/35, respectively. 


Note: 
Exercise: 


Problem: 


Consider the same region Q ({link]), and use the density function p (x, y, z) = xy”z. Find the moments 
of inertia about the three coordinate planes. 


Solution: 


The moments of inertia of the tetrahedron Q about the yz-plane, the xz-plane, and the xy-plane are 
99/35, 36/7, and 243/35, respectively. 


Key Concepts 
Finding the mass, center of mass, moments, and moments of inertia in double integrals: 
¢ Fora lamina R with a density function p (x, y) at any point (2, y) in the plane, the mass is 


m= [fo(anaa 


e The moments about the z-axis and y-axis are 
Equation: 


Me= [f vo(evtasnast = ff so(ex)aa 


i se — My ~ M, 
¢ The center of mass is given by x = —-,y= =. 


The center of mass becomes the centroid of the plane when the density is constant. 
The moments of inertia about the z — axis, y — axis, and the origin are 


Equation: 
LS |[Pote.naa, i= |[Pocnaa, and Ip = I, + I. -|/ (2" + y")p (a, y)dA. 
R R 


Finding the mass, center of mass, moments, and moments of inertia in triple integrals: 


¢ For a solid object Q with a density function p (2, y, z) at any point (z, y, z) in space, the mass is 


m= fff o(.waav. 


e¢ The moments about the zy-plane, the xz-plane, and the yz-plane are 


Equation: 
May = [ff 20(.4,24V, Mes = fff ve (e.4.24V, My = fff 20 (2,y,2aV 
Q Q Q 
e The center of mass is given by x = sale yr Me 2 ies 


The center of mass becomes the centroid of the solid when't the density is constant. 
The moments of inertia about the yz-plane, the xz-plane, and the xy-plane are 


Equation: 
te= fff w+ 2)o@uaavity= fff (+ 2)o@uaav, 
Q Q 


L= II (x? + y”)p(z,y, z)dV. 
Q 


Key Equations 


e Mass of a rae 
k ol 
m= im, 9) omy = tim. 790 play 43)4 = [fev 
4=1 79> tal st 
¢ Moment eel wes -axis 
k ol 

Mz = ay (vis) ma = Jim 37 (us) olin) = / yo(e, ud 

e Moment peu sa y-axis 


k I 
My, = gm OY (= i) my = im DS 


(xis) (wip Yi; )AA =, // rp(z,y)dA 
=1 
R 
e Center of mass of a lamina 


ee 
™ // p(z,y)dA m il oe 


In the following exercises, the region R occupied by a lamina is shown in a graph. Find the mass of R with 
the density function p. 
Exercise: 


Problem: RF is the triangular region with vertices (0,0), (0,3), and (6,0); p(z, y) = zy. 


Solution: 


2 
Exercise: 


Problem: R is the triangular region with vertices (0,0), (1,1), (0,5); p(z,y) =x+y. 


Oo 0204 0608 1%* 
Exercise: 


Problem: R is the rectangular region with vertices (0,0), (0,3), (6,3), and (6,0); p(z, y) = ./zy. 


Solution: 


24,/2 


Exercise: 


Problem: R is the rectangular region with vertices (0, 1), (0, 3), (3, 3), and (3, 1); p(x, y) = xy. 


y 
4 


Exercise: 


Problem: 


R is the trapezoidal region determined by the lines y = — ra + 34y =0,y = 2, and z = 0; 
p(z,y) = 3zy. 


Solution: 


76 
Exercise: 
Problem: 


RF is the trapezoidal region determined by the lines y = 0, y = 1, y = a, and 
y=—2+3;p(z,y) =2r+y. 


Exercise: 


Problem: R is the disk of radius 2 centered at (1,2); p(z,y) =2? + y? — 2x — 4y+ 5. 


A 
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Solution: 


8r 


Exercise: 


Problem: R is the unit disk; p (z, y) = 3x* + 6x7y? + 3y*. 


Exercise: 


Problem: R is the region enclosed by the ellipse x? + 4y? = 1; p(z,y) =1. 


Solution: 


wT 


2 


Exercise: 


Problem: R = {(x,y)|9x? + y? < 1,2 >0,y> 0}; p(x, y) = 9x? + 9. 


Exercise: 


Problem: R is the region bounded by y = z,y = —z,y=2+2,y=—2+2; p(z,y) =1. 


os 


Solution: 


2 


Exercise: 


Problem: RF is the region bounded by y ayy iy 1, and y = 2;p(z,y) =4(x+y). 


In the following exercises, consider a lamina occupying the region R and having the density function p given 
in the preceding group of exercises. Use a computer algebra system (CAS) to answer the following questions. 


a. Find the moments M, and M, about the z-axis and y-axis, respectively. 
b. Calculate and plot the center of mass of the lamina. 
c. [T] Use a CAS to locate the center of mass on the graph of R. 


Exercise: 
Problem: [T] F is the triangular region with vertices (0,0), (0,3), and (6,0); p(x, y) = zy. 


Solution: 


Exercise: 


Problem: [T] F is the triangular region with vertices (0,0), (1,1), and (0,5); p(z,y)=x+y. 
Exercise: 


Problem: 
[T] F is the rectangular region with vertices (0,0), (0,3), (6,3), and (6,0); p(x, y) = ./zy. 
Solution: 


_ 21672 oo ABI Pee. 2 AB See 
aM, = == My = Si bce y= 
C. 


ao 


? 


ws 


Exercise: 


Problem: [T] R is the rectangular region with vertices (0, 1), (0, 3), (3,3), and (3,1); p(z, y) = xy. 
Exercise: 


Problem: 


[T] FR is the trapezoidal region determined by the lines y = —f2 + 3, y=0,y = 2,andz = 0; 
p(z,y) = 3zy. 


Solution: 


_ 368 — 1552, Bice a Oa te coe DOO t 
a. M, = —=,M, = —";b.2= Fy = 3 
Gy 


Exercise: 


Problem: 


[T] FR is the trapezoidal region determined by the lines y = 0, y = 1, y = a, and 
y=—a+3;p(z,y) =2r+y. 


Exercise: 


Problem: [T] R is the disk of radius 2 centered at (1,2); p(z,y) = 2? +y? — 2x —4y+5. 
Solution: 


aM, = 1670; = 87s b: 2 1, 92; 
oF 


> — 
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Exercise: 


Problem: [T] R is the unit disk; p (x, y) = 3x4 + 6x7y? + 3y+. 


Exercise: 


Problem: [T] R is the region enclosed by the ellipse 2? + 4y” = 1; p(x, y) = 1. 


Solution: 


0.2 04 06 0.8 


Exercise: 


Problem: [T] R = {(x,y)|9x? + y? < 1,2 >0,y > 0}; p(z,y) = /9x? + y?. 


Exercise: 


Problem: [T] F is the region bounded by y = z,y = —2,y=2+2, andy = —x+4 2; p(z,y) = 1. 
Solution: 


a Mo 2, M0; bie = 0,1 
Cy 


Noss 
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Exercise: 


Problem: [T] # is the region bounded by y i, y ay 1,and y = 2; p(z,y) =4(@+ y). 


In the following exercises, consider a lamina occupying the region R and having the density function p given 
in the first two groups of Exercises. 


a. Find the moments of inertia I,, I, and Ip about the x-axis, y-axis, and origin, respectively. 
b. Find the radii of gyration with respect to the x-axis, y-axis, and origin, respectively. 


Exercise: 
Problem: RF is the triangular region with vertices (0,0), (0,3), and (6,0); p(z, y) = zy. 
Solution: 


al, = 22.7 


243. J, — 486. and Ip = 8; b. R, = 2¥5, R, = 83, and Ry =3 


Exercise: 


Problem: £ is the triangular region with vertices (0,0), (1,1), and (0,5); p(z,y) =a+y. 


Exercise: 


Problem: R is the rectangular region with vertices (0,0), (0,3), (6,3), and (6, 0); p (a, y) = ./zy. 


Solution: 


ani UAE eee 648/2 vand Ip = = 2?» R, = EE i 3V21 ,and Ry 3Vv105 
Exercise: 
Problem: R is the rectangular region with vertices (0, 1), (0,3), (3, 3), and (3, 1); p(x, y) = xy. 


Exercise: 


Problem: 


RF is the trapezoidal region determined by the lines y = — 42 + 3 ,y=0,y=2, and 
t= 0;p(x,y) = 3zy. 


Solution: 


a. I, = 88, Iy = 1560, and Ip = 1648; b. Rr = “#8, Ry = YO and Ry = 2¥4982 
Exercise: 
Problem: 
RF is the trapezoidal region determined by the lines y = 0, y = 1, y = a, and 
y=—a+3;p(z,y) =2r+y. 


Exercise: 


Problem: R is the disk of radius 2 centered at (1,2); p(z, y) = a? + y? — 2a —4y+5. 


Solution: 


_ 128 _ 56 _ 1847, _ 4v3 _ v2i _ vV69 
al, = 3 dy = 737, and Ip = —3*; b. Ry = =~, Ry = 4, and Ro = > 


Exercise: 


Problem: R is the unit disk; p (x, y) = 3x4 + 6x?y? + 3y%. 


Exercise: 


Problem: R is the region enclosed by the ellipse x” + 4y? = 1; p(a, y) = 1. 


Solution: 
Seep pa =) bre = i Ale — v5 
aI, = 31, = F,andIp = 5b. R, = 7, Ry = 7, and Ry = FY 
Exercise: 


Problem: R = {(zx,y)|9x? + y? < 1,2 >0,y > 0}; p(z,y) = /9x? + y?. 


Exercise: 


Problem: R is the region bounded by y = z,y = —2,y= 2+ 2,andy = —x+ 2; p(az,y) =1. 


Solution: 


Exercise: 


Problem: RF is the region bounded by y ayy yy 1,and y = 2;p(z,y) =4(x+ y). 


Exercise: 


Problem: 


Let Q be the solid unit cube. Find the mass of the solid if its density is equal to the square of the 
distance of an arbitrary point of Q to the xy-plane. 


Solution: 


a3, hl 
ih 3 
Exercise: 


Problem: 


Let Q be the solid unit hemisphere. Find the mass of the solid if its density p is proportional to the 
distance of an arbitrary point of Q to the origin. 


Exercise: 


Problem: 


The solid Q of constant density 1 is situated inside the sphere x? + y” + z” = 16 and outside the 
sphere 2? + y? + z? = 1. Show that the center of mass of the solid is not located within the solid. 


Exercise: 


Problem: 


Find the mass of the solid Q = {(x, y, z)|1 < a7 + 2? <5, <1 — 2? — z’\ whose density is 
p(x, y, z) = k, where k > 0. 


Exercise: 


Problem: 


[T] The solid Q = {(z, y, z)|z? + y? < 9,0 < z<1,2>0,y > 0} has density equal to the distance 
to the zy-plane. Use a CAS to answer the following questions. 


a. Find the mass of Q. 

b. Find the moments M,,, M,z,and M,, about the zy-plane, xz-plane, and yz-plane, respectively. 
c. Find the center of mass of Q. 

d. Graph Q and locate its center of mass. 


Solution: 


am=b My = >, Ma=  My,= Ste= 2, 9= 2.2 = 7d the solid Q and iis 


center of mass are shown in the following figure. 


Exercise: 


Problem: 


Consider the solid Q = {(z, y, z)|0 << @<1,0<y< 2,0 < z < 3} with the density function 
p(2,Y; 2) =ar+yt+1. 


a. Find the mass of Q. 
b. Find the moments M,,, M,z,and M,, about the xy-plane, xz-plane, and yz-plane, respectively. 
c. Find the center of mass of Q. 


Exercise: 


Problem: 
1. ob 
[T] The solid Q has the mass given by the triple integral 7 J / r’dr d0 dz. Use a CAS to answer 
-1 0 0 
the following questions. 


a. Show that the center of mass of Q is located in the zy-plane. 
b. Graph Q and locate its center of mass. 


Solution: 

- yg - __ 3(2-v2) ee : 7 : j 
at = 3>,y = zy, z = 0; b. the solid Q and its center of mass are shown in the following 
figure 


Exercise: 
Problem: 
The solid Q is bounded by the planes x + 4y + z= 8,2 = 0,y = 0, and z = 0. Its density at any 


point is equal to the distance to the xz-plane. Find the moments of inertia I, of the solid about the 
xz-plane. 


Exercise: 
Problem: 
The solid Q is bounded by the planes x + y+ z= 3, 4 = 0,y = 0, and z = 0. Its density is 


p(x, y, Z) = x + ay, where a > 0. Show that the center of mass of the solid is located in the plane 
z= 2 for any value of a. 


Exercise: 


Problem: 


oe Q pe the pone situated outside the sphere x? + y? + z? = z and inside the upper hemisphere 


x? +y? +z? = R?, where R > 1. If the density of the solid is p (2, y, z) = TS find R such 


that the mass of the solid is = 


Exercise: 
Problem: 
2 V4-a2 V/16—-2?-¥ 
The mass of a solid Q is given by i i, : (2? + y%+ Ze) "dz dy dx, where n is an 


0 0 ze +y2 
integer. Determine n such the mass of the solid is (2 — v2) Tr. 


Solution: 


n=-2 
Exercise: 
Problem: 
Let Q be the solid bounded above the cone x? + y? = z” and below the sphere 


x? + y* + 27 — 4z = 0. Its density is a constant k > 0. Find k such that the center of mass of the solid 
is situated 7 units from the origin. 


Exercise: 


Problem: 


The solid Q = {(a,y, z)|(0 < 27 +y < 16,2 >0,y > 0,0 < z < z} has the density 
p(x, y, z) = k. Show that the moment M,, about the zy-plane is half of the moment M,, about the 
yz-plane. 


Exercise: 
Problem: 


The solid Q is bounded by the cylinder x? + y? = a?, the paraboloid b? — z = x? + y’”, and the 
i es where 0 < a < b. Find the mass of the solid if its density is given by 


(x,y,z) = fa? +y?. 
Exercise: 


Problem: 


Let Q be a solid of constant density k, where k > 0, that is located in the first octant, inside the circular 
cone x? + y? = 9(z — 1)”, and above the plane z = 0. Show that the moment Mozy about the zy-plane 
is the same as the moment M,, about the xz-plane. 


Exercise: 


1 1/2 
Problem: The solid Q has the mass given by the triple integral i J rt r) \dz dé dr. 
0 0 0 


a. Find the density of the solid in rectangular coordinates. 
b. Find the moment M,,, about the zy-plane. 


Exercise: 


Problem: 


The solid Q has the moment of inertia J, about the yz-plane given by the triple integral 


2 view Vere 
/ / / (y° + Zz) (2 + y”)dz dx dy. 
0 Vay Het?) 


a. Find the density of Q. 
b. Find the moment of inertia J, about the zy-plane. 


Solution: 


a. p (x,y,z) = 2? + y";b. Bt 
Exercise: 


Problem: 


m/4 2sec@ 1 
The solid @ has the mass given by the triple integral / i (r° cos # sin 86 + 2r) dz dr dé. 
0 oOo oO 


a. Find the density of the solid in rectangular coordinates. 
b. Find the moment M,, about the xz-plane. 


Exercise: 
Problem: 


Let Q be the solid bounded by the xy-plane, the cylinder x? + y? = a”, and the plane z = 1, where 


a > 1is areal number. Find the moment Mz, of the solid about the zy-plane if its density given in 


cylindrical coordinates is p (r, 0, z) = = (r), where f is a differentiable function with the first and 


second derivatives continuous and differentiable on (0, a). 


Solution: 


Mz, = 7 (f (0) — f (a) + af! (a) 
Exercise: 


Problem: 


b 
A solid Q has a volume given by | / / dA dz, where D is the projection of the solid onto the 
D a 


xy-plane and a < bare real numbers, and its density does not depend on the variable z. Show that its 
center of mass lies in the plane z = arb 


Exercise: 
Problem: 
Consider the solid enclosed by the cylinder x? + z? = a? and the planes y = b and y = c, where a > 0 
and b < care real numbers. The density of Q is given by p (a, y, z) = f(y), where f is a differential 


function whose derivative is continuous on (6, c). Show that if f (b) = f (c), then the moment of inertia 
about the xz-plane of Q is null. 


Exercise: 


Problem: 


[T] The average density of a solid Q is defined as Pave = vig Il p(z,y,z)dV = Vo)? where 
Q 


V (Q) and mare the volume and the mass of Q, respectively. If the density of the unit ball centered at 


the origin is p(x, y, z) = ee case use a CAS to find its average density. Round your answer to three 


decimal places. 
Exercise: 


Problem: 


Show that the moments of inertia J, Iy, and I, about the yz-plane, xz-plane, and xy-plane, 


respectively, of the unit ball centered at the origin whose density is p(x, y, z) = e* —¥'—*" are the 


same. Round your answer to two decimal places. 
Solution: 


I, =I, = 1, ~ 0.84 


Glossary 


radius of gyration 
the distance from an object’s center of mass to its axis of rotation 


Change of Variables in Multiple Integrals 


e Determine the image of a region under a given transformation of variables. 
¢ Compute the Jacobian of a given transformation. 

e Evaluate a double integral using a change of variables. 

e Evaluate a triple integral using a change of variables. 


Recall from Substitution Rule the method of integration by substitution. When evaluating an integral such as 


3 
/ a(x? — 4) de, we substitute u = g(x) = x? — 4. Then du = 22 dz oraz dz = +du and the limits change to 
2 


5] 
u = g(2) = 2? —4= 0 andu = g(3) = 9 —4=5. Thus the integral becomes / audu and this integral is 
0 
much simpler to evaluate. In other words, when solving integration problems, we make appropriate substitutions to 


obtain an integral that becomes much simpler than the original integral. 


We also used this idea when we transformed double integrals in rectangular coordinates to polar coordinates and 
transformed triple integrals in rectangular coordinates to cylindrical or spherical coordinates to make the 
computations simpler. More generally, 

Equation: 


/ a / f (9(u))g"(u)du, 


Where z = g(u), dx = gi(u)du, and u = c and u = d satisfy c = g(a) andd = g(b). 


A similar result occurs in double integrals when we substitute z = f (r,0) = rcos 0, y = g(r,@) = rsin 9, and 
dA = dz dy = r dr dO. Then we get 


Equation: 
|] tena = ff t 0086, sin8)r dra 
R s 


where the domain R is replaced by the domain S in polar coordinates. Generally, the function that we use to 
change the variables to make the integration simpler is called a transformation or mapping. 


Planar Transformations 


A planar transformation T is a function that transforms a region G' in one plane into a region R in another plane 
by a change of variables. Both G and R are subsets of R?. For example, [link] shows a region G in the wv-plane 
transformed into a region R in the zy-plane by the change of variables x = g (u,v) and y = h (u,v), or 
sometimes we write z = x (u,v) and y = y(u, v). We shall typically assume that each of these functions has 
continuous first partial derivatives, which means g,,, gy, h,, and hy, exist and are also continuous. The need for this 
requirement will become clear soon. 


xX = g(u, Vv) 
y = hy, v) 


Cartesian uv-plane Cartesian xy-plane 


The transformation of a region G in the wv-plane into a region F in the zy-plane. 


Note: 

Definition 

A transformation T: G — R, defined as T (u, v) = (a, y), is said to be a one-to-one transformation if no two 
points map to the same image point. 


To show that T is a one-to-one transformation, we assume T' (ui, v1) = T (ue, v2) and show that as a 
consequence we obtain (wi, v1) = (we, v2). If the transformation T is one-to-one in the domain G, then the 
inverse T~' exists with the domain R such that T~! o T and T o T~? are identity functions. 


[link] shows the mapping T (u,v) = (a, y) where x and y are related to u and v by the equations x = g (u,v) and 
y = h(u, v). The region G is the domain of T' and the region R is the range of T’, also known as the image of G 
under the transformation T’. 


Example: 
Exercise: 


Problem: 
Determining How the Transformation Works 


Suppose a transformation T is defined as T (r, 0) = (x, y) where = rcos 0, y = r sin 9. Find the image 
of the polar rectangle G = {(r,0)|0 <r < 1,0 < 0 < 7/2} in the r@-plane to a region R in the zy-plane. 
Show that T is a one-to-one transformation in G and find T’“! (a, y). 


Solution: 


Since r varies from 0 to 1 in the r6-plane, we have a circular disc of radius 0 to 1 in the xy-plane. Because 
6 varies from 0 to 7/2 in the r6-plane, we end up getting a quarter circle of radius 1 in the first quadrant of 
the zy-plane ([link]). Hence R is a quarter circle bounded by xz” + y? = 1 in the first quadrant. 


xX =rcosé 
y=rsine 
7 
r 
Cartesian ré-plane Cartesian xy-plane 


A rectangle in the r@-plane is mapped into a quarter circle in the zy-plane. 


In order to show that T is a one-to-one transformation, assume T' (71,01) = T (r2, 02) and show as a 
consequence that (71,01) = (r2, 92). In this case, we have 
Equation: 


T (ri, 91) => T (r2, 92), 
(x1, 41) = (er yi); 


(ricos 61, 71sin 01) = (racos 6, rosin 02), 
rycos0; = r2cos 69,7 sin 6, = resin 9». 
Dividing, we obtain 
Equation: 
ricos@,; __ = rgcos 0 
risin@,; — rgsin@ 
cos6; _ cos 
sin 0, sin 0 
tan@; = tan6, 
0, = 4 


since the tangent function is one-one function in the interval 0 < 0 < 7/2. Also, since 0 < r < 1, we have 
11 = 12,01 = Oy. Therefore, (71,61) = (172,62) and T is a one-to-one transformation from G into R. 


To find T+ (z, y) solve for r, @ in terms of a, y. We already know that r? = x? + y? and tan 0 = x, Thus 
T~* (,y) = (r, 9) is defined as r = »/x? + y? and 6 = tan (4) 


a)* 


Example: 
Exercise: 


Problem: 
Finding the Image under T 


Let the transformation T' be defined by T (u,v) = (x,y) where z = u? — v” and y = uv. Find the image of 
the triangle in the wv-plane with vertices (0,0), (0, 1), and (1, 1). 


Solution: 


The triangle and its image are shown in [link]. To understand how the sides of the triangle transform, call the 
side that joins (0,0) and (0, 1) side A, the side that joins (0,0) and (1, 1) side B, and the side that joins 
(1,1) and (0,1) side C. 


A triangular region in the wv-plane is transformed into an image in the zy-plane. 


For the side A: u = 0,0 < v < 1 transforms to z = —v”, y = 0 so this is the side A/ that joins (—1, 0) and 
(0,0). 


For the side B: u = v, 0 < u < 1 transforms to z = 0, y = u” so this is the side B’ that joins (0, 0) and 
(0,1). 


For the side C:0 < u < 1,v = 1 transforms to x = u? — 1, y = u (hence x = y” — 1) so this is the side 
C’ that makes the upper half of the parabolic arc joining (—1, 0) and (0, 1). 


All the points in the entire region of the triangle in the uv-plane are mapped inside the parabolic region in 
the xy-plane. 


Note: 
Exercise: 


Problem: 
Let a transformation T be defined as T (u, v) = (x, y) where « = u+ v, y = 3v. Find the image of the 


rectangle G = {(u,v):0 <u < 1,0 < v < 2} from the wv-plane after the transformation into a region R 
in the zy-plane. Show that T is a one-to-one transformation and find T~! (BB). 


Solution: 


ees 
and v = 3 


—1 -, a 325 
T(x, y) = (u,v) where u = = 
Hint 


Follow the steps of [link]. 


Jacobians 


Recall that we mentioned near the beginning of this section that each of the component functions must have 
continuous first partial derivatives, which means that gy, gv, hy, and h, exist and are also continuous. A 
transformation that has this property is called a C’! transformation (here C’ denotes continuous). Let 

T (u,v) = (g (u,v), h (u,v)), where z = g(u,v) and y = h (u,v), be a one-to-one C’ transformation. We want 
to see how it transforms a small rectangular region S, Aw units by Av units, in the uv-plane (see the following 
figure). 


A small rectangle S in the wv-plane is transformed into a region R in the zy-plane. 


Since x = g(u,v) and y = h(u, v), we have the position vector r (u,v) = g(u, v)i + h (u, v)j of the image of 
the point (wu, v). Suppose that (uo, vo) is the coordinate of the point at the lower left corner that mapped to 

(x0, yo) = T (uo, vo). The line v = vp maps to the image curve with vector function r (w, vo), and the tangent 
vector at (2, yo) to the image curve is 

Equation: 


: . Ox, Oy, 
Tu = Gu (uo, vo)i + hy (uo, vo)J = Aa + Ses: 


Similarly, the line w = uo maps to the image curve with vector function r (wo, v), and the tangent vector at 
(0, yo) to the image curve is 

Equation: 

Ox Oy 


pees 


Py = Qy (uo, vo)i + hy (uo, Vo)j = Do Bu 


Now, note that 


Equation: 
r (up + Au, v9) — r (uo, vo) 
ry = dim, Ae sor (uo + Au, vo) — r (uo, vo) © Aury. 
Similarly, 
Equation: 
Ax) = 
r, = li nua Aye (959) sor (uo, Yo + Av) — r (uo, v0) © Avry. 


~~ A330 Av 


This allows us to estimate the area AA of the image R by finding the area of the parallelogram formed by the 
sides Avr, and Aur,,. By using the cross product of these two vectors by adding the kth component as 0, the area 
AA of the image R (refer to The Cross Product) is approximately |Aur,, x Avr,| = |ru x ry|AuAv. In 
determinant form, the cross product is 


Equation: 

i j k 
Ox oy 

: gee dx dy Ox 

xa le o -[e Se (Fm soe) 

iL 

ac yg) bw ne 

Ov Ov 


Since |k| = 1, we have AA & |r, x r,|AuAv = (= OY 2 Be St) AuAv, 


Note: 

Definition 

The Jacobian of the C’ transformation T (u,v) = (g (u,v), h (u,v)) is denoted by J (u,v) and is defined by the 
2 x 2 determinant 


Equation: 
Jur) = |2@o|_ |e 3 -(% Oy — Ow a) 
"1 A(ujv)} ja &) du dv dv du)’ 
Ov Ov 
Using the definition, we have 
Equation: 
O(a, 
AA = J(u, v)AuAv = (2,9) AuAv. 
O (u,v) 
Note that the Jacobian is frequently denoted simply by 
Equation: 
d(z,y) 
J (u,v) = . 
a O (u,v) 
Note also that 
Equation: 
ha 0 Ox Ox 
ie oe 7 Ox Oy Ox Oy _ oe av 
de Oy Ou dv Ov Ou dy oy 
Ov ov dv 


(x,y) 
O(u,v) 
first row and partials of y in the second row. 


Hence the notation J (u,v) = suggests that we can write the Jacobian determinant with partials of x in the 


Example: 


Exercise: 


Problem: 
Finding the Jacobian 


Find the Jacobian of the transformation given in [link]. 
Solution: 


The transformation in the example is T (r, 0) = (r cos 6, r sin 8) where « = r cos @ and y = r sin 6. Thus 
the Jacobian is 


Equation: 
Teo A(a,y) _ a aa _ |jcos6 = =—rsin# 
. ~ O(r,8) ~~ | Oy oy! |sin 6 rcos 0 
Or (ey) 
=rcos*6+rsin?é=r (cos?6 + sin’6) = Pp, 
Example: 
Exercise: 
Problem: 


Finding the Jacobian 
Find the Jacobian of the transformation given in [link]. 
Solution: 


The transformation in the example is T' (u,v) = (wu? — v?, wv) where x = u” — v? and y = uv. Thus the 
Jacobian is 


Equation: 
Ox Ox 
Olz,y) |e ov |_ |2u 
J (u,v) = = = Qu? + 2v?. 
( ? ) O(u, v) — & —2v 
Note: 
Exercise: 


Problem: Find the Jacobian of the transformation given in the previous checkpoint: T (wu, v) = (w+ v, 2v). 


Solution: 
Ox Ox 
Oxy) __ Ou ‘Ov a 1 1 _ 
J(u,%) = aay = Oy mln al =" 
Ou Ov 


Hint 


Follow the steps in the previous two examples. 


Change of Variables for Double Integrals 


We have already seen that, under the change of variables T (u,v) = (a, y) where = g(u,v) andy=h(u,v),a 
small region AA in the ry-plane is related to the area formed by the product AwAv in the wv-plane by the 


approximation 
Equation: 


AA = J (u,v)Au, Av. 


Now let’s go back to the definition of double integral for a minute: 
Equation: 


// f(x,y)dA = mila, > > f (wij, Yig AA. 
R 


i=1 j=l 


Referring to [link], observe that we divided the region S in the uv-plane into small subrectangles S;; and we let 
the subrectangles R,; in the zy-plane be the images of $;; under the transformation T (u,v) = (a, y). 


(Xj, Yij) 


The subrectangles $;; in the wv-plane transform into subrectangles R;; in the ry-plane. 


Then the double integral becomes 


Equation: 
|] fea = tim SOY Fleasuy)A = tim OF 
: i=1 j=l i=1 j=l 


Notice this is exactly the double Riemann sum for the integral 


Equation: 
ff fatearroom pea 


(9 (wiz, Vij), (uaz, Vig)) [J (wiz, vig) |AuAv 


au dv. 


Note: 

Change of Variables for Double Integrals 

Let T (u,v) = (x,y) where x = g(u, v) and y = h(u, v) be a one-to-one C' transformation, with a nonzero 
Jacobian on the interior of the region S$ in the wv-plane; it maps S into the region R in the xy-plane. If f is 


continuous on R, then 
= O(z,y) 
i f (2,y)dd = i F(a(uv)h(ue)) | 5% awa 


Equation: 


With this theorem for double integrals, we can change the variables from (2, y) to (w, v) in a double integral 
simply by replacing 

Equation: 

O(z,y) 
O(u, v) 


dA = dx dy = | tudy 


when we use the substitutions z = g(u, v) and y = h(w, v) and then change the limits of integration accordingly. 


This change of variables often makes any computations much simpler. 


Example: 
Exercise: 


Problem: 
Changing Variables from Rectangular to Polar Coordinates 


Consider the integral 


Equation: 
2 V2x—x2 
/ i Va +y? dydz. 
0 0 


Use the change of variables x = r cos # and y = r sin 8, and find the resulting integral. 


Solution: 


First we need to find the region of integration. This region is bounded below by y = 0 and above by 
y = V 2x — x? (see the following figure). 
y 


Changing a region from 
rectangular to polar 
coordinates. 


Squaring and collecting terms, we find that the region is the upper half of the circle x” + y? — 22 = 0, that 
is, y*? + (2 — 1 = 1. In polar coordinates, the circle is r = 2 cos @ so the region of integration in polar 
coordinates is bounded by 0 < r < cos@and0 <@< a 


The Jacobian is J (r,@) = r, as shown in [link]. Since r > 0, we have |J (r, @)| =r. 


The integrand \/ x? + y? changes to r in polar coordinates, so the double iterated integral is 
Equation: 


2 V2a—x? 7/2 2cosd m/2 2cosd 
i i! Ja? + ydy dr = / i r|J (r,0)|dr dé = i / r?dr dé. 
0 0 0 0 0 0 
Note: 
Exercise: 
Problem: 
V1—2? 


1 
Considering the integral i i) (a SF y’)dy dz, use the change of variables x = r cos @ and y = r sin 0, 
0 0 


and find the resulting integral. 


Solution: 


m/2 1 
i [ra do 
0 0 

Hint 


Follow the steps in the previous example. 


Notice in the next example that the region over which we are to integrate may suggest a suitable transformation for 
the integration. This is a common and important situation. 


Example: 
Exercise: 


Problem: 
Changing Variables 


Consider the integral | / (a — y)dy dx, where R is the parallelogram joining the points (1, 2), 
R 


(3, 4), (4,3), and (6, 5) ([link]). Make appropriate changes of variables, and write the resulting integral. 
ne (6,5) 


(3, 4) 


The region of integration for the 
given integral. 


Solution: 


First, we need to understand the region over which we are to integrate. The sides of the parallelogram are 
x-y+1=0,2-y-—-1=0,2-—3y+5=0,andz — 3y+ 9 = 0 ((link)). Another way to look at them 
ise = oj —= ae = — Iles — oi) = ay, aintal ae — sin) = W) 


Clearly the parallelogram is bounded by the lines y= 2+ 1,y=a2—1,y= = (x +5), and y= =(x + 9). 


Notice that if we were to make u = x — yand v = x — 3y, then the limits on the integral would be 
—1<u<land—9<v<-—5B. 


To solve for x and y, we multiply the first equation by 3 and subtract the second equation, 
3u — v = (3a — 3y) — (a — 3y) = 2a. Then we have x = 7 * Moreover, if we simply subtract the 
second equation from the first, we get u — v = (x — y) — (a — 3y) = 2yandy = 45°. 


4 y=x+1 


A parallelogram in the zy-plane that we want to transform by a change in 


variables. 


Thus, we can choose the transformation 
Equation: 


and compute the Jacobian J (u, v). We have 


Equation: 
d(e,y) | 2] 3/2 -1/2) 3 1 1 
De) a a — =—-——+—=-——_, 
A(u,v) |% | |1/2 -1/2 cee 2 
Ou a) 
Therefore, | J (u,v)| = +. Also, the original integrand becomes 
Equation: 
1 1 
£ y= 5 8u v u+vr = >[su u| = —[2u] = u. 
Therefore, by the use of the transformation 7’, the integral changes to 
Equation: 
= il =—§ i ; 
|[@-vdyar - i [1 (ue)udude = / [(G)uaua 
R =) =il =) =il 


which is much simpler to compute. 


Note: 
Exercise: 


Problem: 


4 
ay dz, where R is the trapezoid 
(x — y) 
bounded by the lines  — y = 2,4 — y= 4, x2 = 0, and y = 0. Write the resulting integral. 


Make appropriate changes of variables in the integral | i! 


Solution: 
4 2 
mat — ij AS (1 on 
z= 5(u+u) andy = >(v— 4) and lie au dv. 
=i 29) 
Hint 


Follow the steps in the previous example. 


We are ready to give a problem-solving strategy for change of variables. 


Note: 
Problem-Solving Strategy: Change of Variables 


1. Sketch the region given by the problem in the zy-plane and then write the equations of the curves that form 
the boundary. 

2. Depending on the region or the integrand, choose the transformations x = g(u, v) andy = h (u,v). 

3. Determine the new limits of integration in the wv-plane. 

4. Find the Jacobian J (u, v). 

5. In the integrand, replace the variables to obtain the new integrand. 

6. Replace dy dx or dx dy, whichever occurs, by J (u, v)du dv. 


In the next example, we find a substitution that makes the integrand much simpler to compute. 


Example: 
Exercise: 


Problem: 
Evaluating an Integral 


Using the change of variables u = x — yandv = z + y, evaluate the integral 
Equation: 
R 


where R is the region bounded by the lines z + y = 1 and x + y = 3 and the curves x? — y? = —1 and 
Ne y" = 1 Gee the first region in [link]). 


Solution: 
As before, first find the region R and picture the transformation so it becomes easier to obtain the limits of 


integration after the transformations are made ([link]). 
Yh 


x+y=1 


Transforming the region R into the region S to simplify the computation of an integral. 


Given u = x —yandv=2+ y, we have x = a and y = +;* and hence the transformation to use is 
(0) = (S, =. The lines « + y= 1 anda + y = 3 become v = 1 and v = 3, respectively. The 
curves 2? — y* = 1 and x” — y? = —1 become wv = 1 and uv = —1, respectively. 


Thus we can describe the region S (see the second region [link]) as 
Equation: 


=I 1 
s={ (wit svs3,— sus ct, 
v v 


The Jacobian for this transformation is 


Equation: 
Ox Ox 
, Si) FS 22 oe 


Therefore, by using the transformation T’, the integral changes to 


Equation: 
A 3 1fv 
// (x —y)e™ ¥dA= s/ / ue’ du du 
1 -1/v 
Doing the evaluation, we have 
Equation: 
F 3 I/v F 
3 / / ue“’du dv = — ~ 0.490. 
2, 3e 
1 -1/v 
Note: 
Exercise: 
Problem: 


Using the substitutions z = v and y = /u + v, evaluate the integral | / y sin (y? = «)dA where R is the 


R 
region bounded by the lines y = »/z,x = 2,andy = 0. 


Solution: 
5 (sin 2 — 2) 
Hint 


Sketch a picture and find the limits of integration. 


Change of Variables for Triple Integrals 


Changing variables in triple integrals works in exactly the same way. Cylindrical and spherical coordinate 
substitutions are special cases of this method, which we demonstrate here. 


Suppose that G is a region in www-space and is mapped to D in xyz-space ([link]) by a one-to-one C! 
transformation T (u,v, w) = (x, y, z) where z = g(u,v, w), y= h(u,v, w), andz = k (u,v, w). 


w z 


X = g(u, Vv, w) 
y = Aly, v, w) 
Zz = ku, v, w) 


SS 


Cartesian uvw-space Cartesian xyz-space 


A region G in wuw-space mapped to a region D in xyz-space. 


Then any function F' (x, y, z) defined on D can be thought of as another function H (wu, v, w) that is defined on G: 
Equation: 


F (z,y,z) = F(g(u,v, w), h (u,v, w), k (u,v, w)) = A (u, v, w). 


Now we need to define the Jacobian for three variables. 


Note: 
Definition 
The Jacobian determinant J (uw, v, w) in three variables is defined as follows: 
Equation: 
Ox oy Oz 
Ou Ou Ou 
_— | dx oy Oz 
J(u, Us w) Ov Ov Ov |° 
Ox oy Oz 
Ow dw Ow 


This is also the same as 
Equation: 


Ou Ov Ow 
— | Oy Oy Oy 
J(u, Oo, w) Ou Ou Ow |" 
Oz Oz Oz 
Ou Ov Ow 


The Jacobian can also be simply denoted as ae 


With the transformations and the Jacobian for three variables, we are ready to establish the theorem that describes 
change of variables for triple integrals. 


Note: 

Change of Variables for Triple Integrals 

Let T (u,v, w) = (2, y, z) where z = g(u,v,w),y = h(u,v,w), and z = k(u, v, w), be a one-to-one C! 
transformation, with a nonzero Jacobian, that maps the region G in the uvw-plane into the region D in the 
xyz-plane. As in the two-dimensional case, if F' is continuous on D, then 

Equation: 


[ff Ftemsar = [ff Fotwowy nw seik(us)|srerray tude 


_ |[[# scan) [si eaegcan|febale ai. 
G 


Let us now see how changes in triple integrals for cylindrical and spherical coordinates are affected by this 


Example: 
Exercise: 


Problem: 
Obtaining Formulas in Triple Integrals for Cylindrical and Spherical Coordinates 


Derive the formula in triple integrals for 


a. cylindrical and 
b. spherical coordinates. 


Solution: 


a. For cylindrical coordinates, the transformation is T(r, 0, z) = (a, y, z) from the Cartesian r6z-plane to 
the Cartesian zyz-plane ((link]). Here x = rcos 0, y = rsin 9, and z = z. The Jacobian for the 
transformation is 
Equation: 


Ox Ox Ox 


‘Or 60 Oz 
_ Ax,y,z) __ | dy Oy Oy 
J(r,9,2) = Gaz Or BG CBe 
Oz Oz 
ar 30 Bz 
cos 0 —rsin@ 0 
=lsin0  rcos@ 0| = rcos?6+rsin26 =r (cos?6 4b sin’@) =T. 
0 0 il 


We know that r > 0, so | J (r, 0, z)| = r. Then the triple integral is 


Equation: 
[ff ten Z)dV = Jf #0050, rsin 0, 2)r dra de 
D G 
Cube with sides 
parallel to the 
z coordinate axes Zz 
Z = constant 
D 
xX = rcosé 
y =rsino 
zZ=Z 
r = constant 
@ = constant 
x 
Cartesian réz-space Cartesian xyz-space 


The transformation from rectangular coordinates to cylindrical coordinates can be treated as a 
change of variables from region G in r0z-space to region D in xyz-space. 


oO 


. For spherical coordinates, the transformation is T (p, 6, y) = (2, y, z) from the Cartesian pOy-plane to 
the Cartesian xyz-plane ((link]). Here x = psin y cos 0, y = psin y sin 0, and z = pcos wy. The 
Jacobian for the transformation is 


Equation: 
cee 
a( ’ op 08 dp sinycos9  —psinysinOd pcosypcosé 
Lez : : ; : 
J, 2) = a(p 5 = 7 oy ze =|singsinO —psinycos@ pcosysin§@|. 
a Oz Oz Oz cos 8 0 —psiny 
Op 00 Op 


Expanding the determinant with respect to the third row: 
Equation: 


—psinysinO pcosycosé singycos@ -—psinysin#@ 


= COS 


psin y sin 0 pcos y sin 8 snysin@O psinwcosé 
= cos Y (—p’sin y cos y sin? — p’sin y cos y cos”6) 
—psiny (0 sin cos?6 + psin?y sin”@) 
= —p’sin y cos? (sin70 = cos’6) — p’sin y sin*p (sin?0 + cos’6) 
= —p’sin y cos*y — p*sin y sin? 


=p sini (cosy = sin’) = =p sino. 


Since 0 < y < zm, we must have sin y > 0. Thus |J (p, 6, »)| = |—p’sin y| = p’sin . 
Cube with sides 
parallel to the @ = constant 
g coordinate axes z 


p = constant 


X =psingcosé 
y=psingsing 
Z=pcos¢ 


Cartesian pgé-space Cartesian xyz-space 


The transformation from rectangular coordinates to spherical coordinates can be treated as a 
change of variables from region G in p0y-space to region D in xyz-space. 


Then the triple integral becomes 
Equation: 


[]f swaav = fff t(esin ¢ e056, psin g sind, pos ¢)p*sin y dp dp d8. 
D G 


Let’s try another example with a different substitution. 


Example: 
Exercise: 


Problem: 
Evaluating a Triple Integral with a Change of Variables 


Evaluate the triple integral 


Equation: 


3 4 (y/2)+1 


(« ar =) de dy dz 
0 (0 y/2 


in xyz-space by using the transformation 
Equation: 


u = (2a — y)/2,v = y/2, andw = z/3. 


Then integrate over an appropriate region in www-space. 
Solution: 


As before, some kind of sketch of the region G in xyz-space over which we have to perform the integration 
can help identify the region D in uvw-space ([link]). Clearly G in xyz-space is bounded by the planes 
z=y/2,c¢ = (y/2)+1,y=0, y = 4, z = 0, and z = 4. We also know that we have to use 

u = (2x — y)/2,v = y/2, and w = z/3 for the transformations. We need to solve for x, y, and z. Here we 
find that 7 = w-- vu, yi— 2, and z — Sw. 


Using elementary algebra, we can find the corresponding surfaces for the region G and the limits of 
integration in wuw-space. It is convenient to list these equations in a table. 


Equations in xyz for the Corresponding equations in ww for Limits for the 
region D the region G integration in wow 
cy 2 u+v=2u/2=v u=0 

i= v2 U0 — (20/2) 4 =o, | Ul 

y=0 2vu=0 v=0 

y=A4 2v=4 U— 2 

z=0 3w = 0 w=0 


p= so 3w = 3 ay = Il 


Rear plane: 
x= Z,ory = 2x 


Front plane: 
=2+1ory=2x-2 


The region G' in wuw-space is transformed to region D in ryz-space. 


Now we can calculate the Jacobian for the transformation: 


Equation: 
Ox Ox Ox 
Ou Ov Ow 1 1 0 
F(ujv,u)=|% % &=|9 2 of=6 
oz Oz Oz 0 oO 3 
du dv dw 


The function to be integrated becomes 
Equation: 


z 3w 
fayz=ats St a Aare =ut+v+w. 


We are now ready to put everything together and complete the problem. 
Equation: 


3 4 
Zz 
ah / (w+ =)de dy dz 
0 
Tia eet 1p 
=f J foorrrmutenmiddede= ff [corre etd 


oO © 
il 


0 0 
on a Log a 
=i) (utv+ujdudedu=6f [|S + ou+wu| du dw 


0 0 0 0 0 0 
2 


i 9 1 
-6f (5 +0 w)dodw = 6 f [5 +S tw] dw 


0 0 0 0 


= [e+ 20ydu= fous u* pars 
0 


Note: 
Exercise: 


Problem: Let D be the region in xyz-space defined by 1 < x < 2,0< ay<2,and0<z< 1. 


Evaluate /I/ (xy ae 3zyz) dx dy dz by using the transformation u = xz, v = zy, and w = 3z. 


Solution: 


Boy 9 
If{G+ du dv dw = 2+1n8 
0 1 


0 
Hint 


Make a table for each surface of the regions and decide on the limits, as shown in the example. 


Key Concepts 


e A transformation T’ is a function that transforms a region G in one plane (space) into a region # in another 
plane (space) by a change of variables. 

¢ A transformation T : G — R defined as T (u, v) = (a, y) (or T (u,v, w) = (2, y, z)) is said to be a one-to- 
one transformation if no two points map to the same image point. 


¢ If f is continuous on R, then |[seuaa = |[ te (u,v), h(u, v)) | cay 
R Ss 


= 


= 
eS 


e If Fis continuous on R, then 
Equation: 


fF ge = [ff Po (u,v, 1), h (u,v, w), (us v5 w)) | ass 


= JI H (u,v, w) |J (u, v, w)|du du dw. 
G 


In the following exercises, the function T : S — R,T (u,v) = (a, y) on the region 
S ={(u,v)|0 <u <1,0 < v < 1} bounded by the unit square is given, where R C R? is the image of S under 
so 


a. Justify that the function T' is a C! transformation. 
b. Find the images of the vertices of the unit square S through the function T. 
c. Determine the image R of the unit square S and graph it. 


Exercise: 


Problem: x = 2u, y = 3v 


Exercise: 
Problem: z = >,y = 


Solution: 


a. T (u,v) = (g (u,v), h(u,v)), 2 = g(u,v) = 5 andy = h(u,v) = 3. The functions g and h are 
continuous and differentiable, and the partial derivatives g, (u,v) = +, 


- : cas = Oph ha = sae (u,v) = $ are continuous on S; b. T' (0,0) = (0,0), 


= (5, 0),T = (0, +) and T' (1,1) = (5 ; +); c. R is the rectangle of vertices 

a = (i 5 30); (c ae ar 1) in the zy-plane; the following figure. 

y 
0.4 
0.3 
0.2 
0.1 

0 0.1 0.2 0.3 0.4 0.5* 

Exercise: 


Problem: z = u—v,y=u+vu 


Exercise: 


Problem: z = 2u — v,y = u+ 2v 


Solution: 


a. T (u,v) = (g (u,v), h(u,v)), 2 = g(u,v) = 2u—v, andy =h(u, aes The functions g and h 
are continuous and differentiable, and the partial derivatives a ,v) = 2, gv ee v) = -1, ha (u,v) = 1, 
and hy (u,v) = 2 are continuous on S; b. T' (0,0) = (0,0), T (1,0) = (2,1), T (0,1) = (—1,2), and 

T (1,1) = (1,3); c. B is the parallelogram of vertices (0, 0), (2,1), (1,3), and (—1, 2) in the zy-plane; see 
the following figure. 


y 
3 


—2 -1 0 1 2* 
Exercise: 


Problem: x = wu”, y = v” 

Exercise: 
Problem: x = u®, y = v® 
Solution: 
a. T (u,v) = (g(u,v),h(u,v)), 2 = g(u,v) =u, and y = h(u,v) = v®. The functions g and h are 
continuous and differentiable, and the partial derivatives g,, (u,v) = 3u?, g, (u,v) = 0, h, (u,v) = 0, and 
hy (u,v) = 3v? are continuous on S; b. T (0,0) = (0,0), T (1,0) = (1,0), T (0, 1) = (0,1), and 


T (1,1) = (1,1); c. B is the unit square in the zy-plane; see the figure in the answer to the previous 
exercise. 


In the following exercises, determine whether the transformations T’: S — R are one-to-one or not. 
Exercise: 


Problem: z = u”, y = v, where S is the rectangle of vertices (—1, 0), (1,0), (1,1), and (—1, 1). 


Exercise: 
Problem: z = u‘, y = u + v, where S is the triangle of vertices (—2, 0), (2,0), and (0, 2). 


Solution: 


T is not one-to-one: two points of S have the same image. Indeed, T (—2,0) = T (2,0) = (16, 4). 


Exercise: 


Problem: « = 2u, y = 3v, where S is the square of vertices (—1, 1), (—1, —1), (1, —1), and (1,1). 


Exercise: 


Problem: T (u,v) = (2u — v, wu), where S is the triangle of vertices (—1, 1), (—1, —1), and (1, —1). 
Solution: 


T is one-to-one: We argue by contradiction. T (ui, v1) = T (uz, v2) implies 2u, — v1 = 2u2g — v2 and 
Uy = U2. Thus, uy = ug and vy = vo. 


Exercise: 


Problem: « = u+v+w,y=u-+v,z = w, where S = R=R?. 
Exercise: 

Problem: x = u? +v+w,y =u? + v, z= w, where S = R= R?. 

Solution: 


T is not one-to-one: T (1, v, w) = (—1,v, w) 


In the following exercises, the transformations T’: S — # are one-to-one. Find their related inverse 
transformations T’' : R — S. 
Exercise: 


Problem: x = 4u, y = 5v, where S = R= R?. 
Exercise: 

Problem: x = u + 2v,y = —u+ v, where S = R=R?. 

Solution: 

ah 


Exercise: 


Problem: x = e”"*”, y = e””, where S = R? and R = {(z, y)|a > 0, y > 0} 
Exercise: 

Problem: z = Inu, y = In (uv), where S = {(u, v)|u > 0,v > 0} and R= R?. 

Solution: 


u=e*,vu=e 7 


Exercise: 


Problem: « = u+v-+w,y = 3v,z = 2w, where S = R=R’. 


Exercise: 


Problem: « = u+v,y = v+w,z=u-+uw, where S = R= R’. 


Solution: 


In the following exercises, the transformation T : S + R,T (u,v) = (a, y) and the region R C R? are given. 
Find the region S C R?. 
Exercise: 


Problem: z = au, y = bv, R = {(z,y)|x? + y? < a7b"}, where a,b > 0 


Exercise: 


Problem: x = au, y = bv, R {(z,y)|2 + 4 <1}, where a,b > 0 


Solution: 
S= {(u, v) |u? +ur< 1} 
Exercise: 


U. 


*y=t,2=4, R= {(z,y)|x? +y? +27 < 1}, wherea,b,c >0 


Problem: x 


Exercise: 


Problem: x = au, y = bv,z = cw, R {(c,y)|& — & - & <1,2>0}, wherea,d,c > 0 


Solution: 


R= {(u,v, w)|u? — v? — w? <lw> Oo} 


In the following exercises, find the Jacobian J of the transformation. 
Exercise: 


Problem: « = u+ 2v,y = —u+v 


Exercise: 


3 
Problem: z = + ,y= 


Solution: 

3 

2 
Exercise: 

Problem: z = e yore 
Exercise: 

Problem: z = ue’, y = e~ 

Solution: 


1 


Exercise: 


Problem: z = ucos(e”), y = usin (e”) 
Exercise: 

Problem: z = vsin (u’), y = VCOS (u”) 

Solution: 

2uv 


Exercise: 


Problem: z = ucoshv, y = usinhv, z = w 


Exercise: 
Problem: xz = v cosh (=). y = vsinh (—), z=ut+u’? 


Solution: 


UV 


uz 


Exercise: 


Problem: zc =u+v,y=v+w,z=u 


Exercise: 


Problem: c = u—v,y=u+v,z=ut+tv+w 
Solution: 


2 


Exercise: 


Problem: The triangular region R with the vertices (0,0), (1,1), and (1, 2) is shown in the following figure. 


0 0.2 0.4 0.6 0.8 1* 


a. Find a transformation T : S > R, T (u,v) = (x,y) = (au + bu, cu + dv), where a, b,c, and d are real 
numbers with ad — bc # 0 such that T~! (0,0) = (0,0), 7+ (1, 1) = (1,0), and T? (1, 2) = (0,1). 


b. Use the transformation T to find the area A (R) of the region R. 
Exercise: 
Problem: The triangular region R with the vertices (0,0), (2,0), and (1,3) is shown in the following figure. 


y 
3 


oO 


1 2* 
a. Find a transformation T : S + R, T (u,v) = (x,y) = (au + bv, cu + dv), where a, b,c and d are real 
numbers with ad — be # 0 such that T’~! (0,0) = (0,0), T’~! (2,0) = (1,0), and T’~! (1,3) = (0,1). 
b. Use the transformation T to find the area A (R) of the region R. 


Solution: 


a. T (u,v) = (2u + v, 3v); b. = area of R is 


7 y)/3 5 1 1-wu 
A= / / a= ian (2,9) avdu= ff sdvdu—3. 
0 (u, v) 
0 y¥/3 0 0 


In the following exercises, use the transformation u = y — x, v = y, to evaluate the integrals on the parallelogram 
R of vertices (0,0), (1,0), (2,1), and (1, 1) shown in the following figure. 


y 
1 
0.8 
0.6 


0.4 


0.2 


0 0.5 1 15 


Exercise: 


Problem: // (y—a)dA 
R 


Exercise: 
Problem: | fl (y? _ zy) dA 
R 


Solution: 


wl 


In the following exercises, use the transformation y — z = u, x + y = v to evaluate the integrals on the square R 
determined by the lines y = 2, y = —4% + 2,y = x + 2, and y = —z shown in the following figure. 


y 
2 


=i —0.5 0 0.5 
Exercise: 
Problem: | / en dA 
R 
Exercise: 
Problem: | / sin (a — y)dA 
R 
Solution: 
—1+cos2 


In the following exercises, use the transformation z = u, 5y = v to evaluate the integrals on the region R bounded 
by the ellipse x” + 25y” = 1 shown in the following figure. 


Exercise: 


Problem: // v/a? + 25y2 dA 
R 


Exercise: 
2 
Problem: | / (a? + 25y”) dA 
R 


Solution: 
7 


15 


In the following exercises, use the transformation u = 2 + y, v = x — y to evaluate the integrals on the 
trapezoidal region R determined by the points (1,0), (2,0), (0,2), and (0, 1) shown in the following figure. 


y 
2 


0 


Exercise: 


Problem: | / (a? — Qey+ y”) et 4dA 
R 


Exercise: 


Problem: // (x? + 3a7y + 3ay? + y°) dA 
R 


Solution: 


31 


5 
Exercise: 

Problem: 

The circular annulus sector R bounded by the circles 42? + 4y* = 1 and 9x? + 9y? = 64, the line 


x = yV3, and the y-axis is shown in the following figure. Find a transformation T' from a rectangular region 
S in the r6-plane to the region R in the zy-plane. Graph S. 


y 
3 


2.5 


0 05 115 2 25 3% 


Exercise: 
Problem: 
The solid R bounded by the circular cylinder x? + y* = 9 and the planes z = 0,z = 1, x = 0, andy = Ois 


shown in the following figure. Find a transformation T from a cylindrical box S in r@z-space to the solid R 
in xyz-space. 


Solution: 


T (r,0,z) = (r-cos6,rsin 6, z);S = [0,3] x [0,2] x [0,1] inthe réz-space 


Exercise: 


Problem: 


= 1 
2 2 _ 
Show that |[s (Vz - aa = avi | F (pp dp, where f is a continuous function on (0, 1] and 
R 0 


R is the region bounded by the ellipse 5x? + 3y? = 15. 
Exercise: 


Problem: 


1 
Show that I f (v 16a? +- 4y? + 2)dv =7F (E (p)p7dp, where f is a continuous function on (0, 1] 
R 0 


and R is the region bounded by the ellipsoid 16a? + 4y? + z? = 1. 
Exercise: 
Problem: 
[T] Find the area of the region bounded by the curves zy = 1, ry = 3, y = 2a, and y = 3z by using the 


transformation u = ry and v = a Use a computer algebra system (CAS) to graph the boundary curves of 
the region R. 


Exercise: 


Problem: 


[T] Find the area of the region bounded by the curves x?y = 2, x?y = 3, y = 2, and y = 2z by using the 
transformation u = x7y and v = ae Use a CAS to graph the boundary curves of the region R. 


Solution: 


The area of R is 10 — 4/6; the boundary curves of R are graphed in the following figure. 


0 "08 0.9 1.0 11 12 13 1.4* 


Exercise: 


Problem: 


z+1 


1 2 
Evaluate the triple integral / / / (y + 1)dx dy dz by using the transformation u = x — z, 
01 


z 


v = $y, and w= 4. 


Exercise: 


Problem: 
2 6 32z+2 

Evaluate the triple integral / i / (5 — 4y)dx dz dy by using the transformation 
0 4 32 


u= 2 — 3z,v = 4y,andw = z. 
Solution: 


8 
Exercise: 
Problem: 
A transformation T : R? + R?,T (u,v) = (2, y) of the form x = au + bu, y = cu + dv, where 


a, b,c, and d are real numbers, is called linear. Show that a linear transformation for which ad — bc 4 0 
maps parallelograms to parallelograms. 


Exercise: 
Problem: 
The transformation Ty : R? > R?, Ty (u, v) = (2, y), where x = ucos 6 — vsin 0, y = usin + vcos8, 


is called a rotation of angle 6. Show that the inverse transformation of Tg satisfies Ty t= T_9, where T_9 is 
the rotation of angle —0. 


Exercise: 
Problem: 
[T] Find the region S in the uv-plane whose image through a rotation of angle 7 is the region R enclosed by 
the ellipse x? + 4y? = 1. Use a CAS to answer the following questions. 
a. Graph the region S. 
b. Evaluate the integral | / e *“’du dv. Round your answer to two decimal places. 


Exercise: 
Problem: 
[T] The transformations T; : R? > R?,i = 1,..., 4, defined by Ty (u,v) = (u, —v), 


T2 (u,v) = (—u, v), T3 (u, v) = (—u, —v), and Ty (u, v) = (v, w) are called reflections about the 
x-axis, y-axis, origin, and the line y = z, respectively. 


a. Find the image of the region S = {(u, v)|u? + v? — 2u — 4u + 1 < 0} in the zy-plane through the 
transformation T; 0 T) o T3 0 T4. 
b. Use a CAS to graph R. 


c. Evaluate the integral sin(u”)du dv by using a CAS. Round your answer to two decimal places. 


Solution: 


a R= { (2, y)|y? te? —Qy—4da+1< 0}; b. R is graphed in the following figure; 


c. 3.16 
Exercise: 
Problem: 
[T] The transformation T;,1,1 : R? > R?,Ti1,1 (u,v, w) = (2, y, z) of the form z = ku, y= v,z= wu, 
where & # 1 is a positive real number, is called a stretch if k > 1 and a compression if 0 < k < 1 inthe 


x-direction. Use a CAS to evaluate the integral i i / e (427+ 99° +252") Gp dy dz on the solid 


Ss 
S = {(z,y, z)|4a? + 9y? + 252? < 1} by considering the compression T2,3,5 (u,v, w) = (a, y, z) defined 
by x = >,y = 3, and z= =. Round your answer to four decimal places. 


Exercise: 


Problem: 


[T] The transformation Tyo : R? > R’, Tao (u,v) = (u+ av, v), where a # 0 is areal number, is called a 
shear in the z-direction. The transformation, Tyo : R? > R?,T.» (u,v) = (u, bu + v), where b 4 Oisa 
real number, is called a shear in the y-direction. 


a. Find transformations To o T3 9. 


b. Find the image # of the trapezoidal region S bounded by u = 0,v = 0,v = 1, andv = 2 — wu through 
the transformation To 2 0 T’3 9. 


c. Use a CAS to graph the image R in the zy-plane. 
d. Find the area of the region R by using the area of region S. 
Solution: 


a. To,2 0 T3,9 (u, v) = (u+ 3v, 2u + 7v); b. The image S is the quadrilateral of vertices 
(0, 0), (3,7), (2,4), and (4, 9); c. S is graphed in the following figure; 


3 
d.3 


Exercise: 


Problem: 


Use the transformation, x = au, y = av, z = cw and spherical coordinates to show that the volume of a 
arty? 2 «. Ana?c 
ae ta als 3". 


region bounded by the spheroid 
Exercise: 


Problem: 


Find the volume of a football whose shape is a spheroid ae + £ = 1 whose length from tip to tip is 11 


inches and circumference at the center is 22 inches. Round your answer to two decimal places. 


Solution: 

202 ~ 282.45 in® 
Exercise: 

Problem: 


[T] Lamé ovals (or superellipses) are plane curves of equations (2) _ (4) " — 1, where a, b, and n are 
positive real numbers. 


a. Use a CAS to graph the regions R bounded by Lamé ovals fora = 1,b = 2,n = 4andn = 6, 
respectively. 

b. Find the transformations that map the region R bounded by the Lamé oval «4 + y* = 1, also called a 
squircle and graphed in the following figure, into the unit disk. 


c. Use a CAS to find an approximation of the area A(R) of the region R bounded by x* + y* = 1. Round 
your answer to two decimal places. 


Exercise: 


Problem: 


[T] Lamé ovals have been consistently used by designers and architects. For instance, Gerald Robinson, a 
Canadian architect, has designed a parking garage in a shopping center in Peterborough, Ontario, in the shape 
of a superellipse of the equation (2) we (¥) " = 1 with -= $ and n = e. Use a CAS to find an 
approximation of the area of the parking garage in the case a = 900 yards, b = 700 yards, and n = 2.72 
yards. 


Solution: 


A(R) ~ 83,999.2 


Chapter Review Exercises 


True or False? Justify your answer with a proof or a counterexample. 


Exercise: 
b d d b 
Problem: | [ f(c,y)dydx = [| #(0,y)ayae 
a Cc c a 
Exercise: 


Problem: Fubini’s theorem can be extended to three dimensions, as long as f is continuous in all variables. 


Solution: 


True. 


Exercise: 


2rn 1 1 


Problem: The integral / / i, dz dr dé represents the volume of a right cone. 
0 oO fr 


Exercise: 


Problem: The Jacobian of the transformation for z = u? — 2v, y = 3u — 2uv is given by —4u? + 6u + 4v. 
Solution: 


False. 


Evaluate the following integrals. 
Exercise: 


Problem: // (5a°y? — y*)dA, R= {(2,y)|0< 2 <2,1<y<4} 
R 


Exercise: 


y 
Problem: —“"__ dA, D= {(a2,y)]0<a<1l,-r<y<a 
// 3e2 41 {(z,y)| y Ni 


Solution: 


0 


Exercise: 


Problem: | / sin (a? + y’)dA where D is a disk of radius 2 centered at the origin 
D 


Exercise: 


1 1 
Problem: / / aye® dx dy 
0 y 


Solution: 


1 


4 
Exercise: 


l 2 2=zZ 


Problem: / i i. 6dy dx dz 


aD oD 


Exercise: 


Problem: | ]f suav, where r= {(z,y,2)|0<@<10<y<2,0<z< 9—y} 
R 


Solution: 


1.475 


Exercise: 


2 
Problem: / 
0 


Exercise: 


I m/2 3 

Problem: / / [ esin(o)dp ap a0 
0 0 1 

Solution: 


52 


3 0 


Exercise: 


1 View V1-o 
Problem: / : i dz dy dx 
0 Vine? —./1-2?-y 


For the following problems, find the specified area or volume. 
Exercise: 


Problem: The area of region enclosed by one petal of r = cos (48). 


Solution: 


ome 
16 


Exercise: 


Problem: The volume of the solid that lies between the paraboloid z = 2”? + 2y? and the plane z = 8. 


Exercise: 


Problem: The volume of the solid bounded by the cylinder x? + y” = 16 and from z = 1toz +2 =2. 


Solution: 


93.291 
Exercise: 


Problem: 


The volume of the intersection between two spheres of radius 1, the top whose center is (0,0, 0.25) and the 
bottom, which is centered at (0, 0, 0). 


For the following problems, find the center of mass of the region. 
Exercise: 


Problem: p(x, y) = zy on the circle with radius 1 in the first quadrant only. 


Solution: 


(a5 a5) 


Exercise: 


Problem: p(z, y) = (y + 1),/z in the region bounded by y = e*, y = 0, andz = 1. 


Exercise: 


Problem: p(x, y, z) = z on the inverted cone with radius 2 and height 2. 


Solution: 


(0,0, ) 
Exercise: 


Problem: 


The volume an ice cream cone that is given by the solid above z = ,/(x? + y”) and below 
2+arty? =z. 


The following problems examine Mount Holly in the state of Michigan. Mount Holly is a landfill that was 
converted into a ski resort. The shape of Mount Holly can be approximated by a right circular cone of height 1100 
ft and radius 6000 ft. 

Exercise: 


Problem: 


If the compacted trash used to build Mount Holly on average has a density 400 lb/ ft?, find the amount of 
work required to build the mountain. 


Solution: 


1.4527 x 101° ft-lb 
Exercise: 
Problem: 
In reality, it is very likely that the trash at the bottom of Mount Holly has become more compacted with all 
the weight of the above trash. Consider a density function with respect to height: the density at the top of the 


mountain is still density 400 lb/ ft and the density increases. Every 100 feet deeper, the density doubles. 
What is the total weight of Mount Holly? 


The following problems consider the temperature and density of Earth’s layers. 
Exercise: 


Problem: 
[T] The temperature of Earth’s layers is exhibited in the table below. Use your calculator to fit a polynomial 


of degree 3 to the temperature along the radius of the Earth. Then find the average temperature of Earth. 
(Hint: begin at 0 in the inner core and increase outward toward the surface) 


Layer Depth from center (km) Temperature “C 


Rocky Crust 0 to 40 0 
Upper Mantle 40 to 150 870 
Mantle 400 to 650 870 
Inner Mantel 650 to 2700 870 
Molten Outer Core 2890 to 5150 4300 
Inner Core 5150 to 6378 7200 


Source: http://www.enchantedlearning.com/subjects/astronomy/planets/earth/Inside.shtml 


Solution: 


y = —1.238 x 10-‘x* + 0.001196x? — 3.666z + 7208; average temperature approximately 2800°C’ 
Exercise: 
Problem: 


[T] The density of Earth’s layers is displayed in the table below. Using your calculator or a computer 
program, find the best-fit quadratic equation to the density. Using this equation, find the total mass of Earth. 


Layer Depth from center (km) Density (g/cm3) 
Inner Core 0 12.95 

Outer Core 1228 11.05 

Mantle 3488 5.00 

Upper Mantle 6338 3.90 

Crust 6378 2.55 


Source: http://hyperphysics. phy-astr.gsu.edu/hbase/geophys/earthstruct.html 


The following problems concern the Theorem of Pappus (see Moments and Centers of Mass for a refresher), a 
method for calculating volume using centroids. Assuming a region R, when you revolve around the x-axis the 
volume is given by V, = 27Ay, and when you revolve around the y-axis the volume is given by V, = 27 Az, 
where A is the area of R. Consider the region bounded by x” + y? = 1 and abovey = z+ 1. 

Exercise: 


Problem: Find the volume when you revolve the region around the z-axis. 


Solution: 


COE) 


Exercise: 


Problem: Find the volume when you revolve the region around the y-axis. 


Glossary 
Jacobian 
the Jacobian J (u, v) in two variables isa 2 x 2 determinant: 
Equation: 
Or Oy 
_ | Ou Ou |, 
J (u,v) = i oy 
Ov Ov 


the Jacobian J (wu, v, w) in three variables isa 3 x 3 determinant: 


Equation: 
Ox oy Oz 
Ou Ou Ou 
—. || dz Oy Oz 
J(u, U5 w) Ov Ov Ov 
Ox oy Oz 


Ow Ow Ow 


one-to-one transformation 
a transformation T : G — R defined as T (u, v) = (a, y) is said to be one-to-one if no two points map to the 
same image point 


planar transformation 
a function T that transforms a region G in one plane into a region R in another plane by a change of variables 


transformation 
a function that transforms a region G in one plane into a region R in another plane by a change of variables 


